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On the Design of Nozzles —— 


E OBTAINED an exact particular solution of the non- 

linear equation, which is the principal portion of the 
exact equation of three-dimensional gas flow over a large 
range of Mach numbers (0 < M < 1.7). The result is ap- 
plicable to the design of nozzles. 

1 Let the gas move essentially along the z axis. Let us 
imagine, for example, sufficiently long nozzles, which become 
an elementary small jet in the limit, and assume that we want 
to find more rigorous solutions over a larger range of Mach 
number variation than obtainable under the assumption of a 
small jet. Itis clear that for such a flow it is possible to dis- 
card all the terms with mixed derivatives gy, Gre and ys 
from the exact equation of motion of gas 


u? v? w? 
2uv 2uw 2vw 
Pry — — =0 [1.1] 


by virtue of the smallness of v = g, and w = ¢,, and it is 
also possible to neglect the quantities v? and w* in the co- 
efficients of 9.2, gy, and ¢,,. We thus arrive at the equation 


F(u)gzz + + Ga = 0 


or 
F [ex(1 + + =0 
G = u = cx(1 + [1.2] 
where 
cx? — u? k+1 
Gz F(u) ce? — (: [1.3] 


and cx is the critical velocity of sound. The steeper the walls 
of the nozzle, the less exact is Equation [1.2]. However, 
even here Equation [1.2] is the principal part of the exact 
Equation [1.1] in the vicinity of the velocity of sound and in 
the region of small subsonic velocities. Apparently, the best 
results can be obtained with Equation [1.2] in the design of 
an axially symmetrical nozzle, for in this case definite veloci- 
ties are reached at the minimum expansion. The solution 
can be obtained, however, only for the equation 


—[(K + 1)u’ + Ku’ + + Gu’ = 0 [1.4] 


In the figure, curve 1 represents the exact function F, line 2 
represents —(xk + 1)u’, which corresponds to the Tricomi 
equation, and curve 3 represents the function —[(« + 1)u’ 
+ xu’?]. In the vicinity of u = cy (i.e., u’ = 0), Equation 
[1.4], like Equation [1.2], is equivalent to the Tricomi equa- 
tion. 

2 As in (1),! we seek for Equation [1.4] a particular solu- 
tion of the form 

= + foly, + fily, 2) [2.1] 

Substituting [2.1] into [1.4] and setting equal to zero the 

Translated from Jzvestiia Akademii Nauk SSSR, Otdelenie 
Tekhnicheskikh Nauk, Mekhanika i Mashinostroenie (Bull. USSR 
Acad. Sci., Div. Tech. Sci., Mechanics and Machine Building), 
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coefficients of the powers of x, we obtain the system of equa- 
tions 


oy? + dz? = 
dy? + dz? 8kfi7fe 4(k + 1)fi? 0 [2.2] 
2fil(k + L)fe + = 0 


dy? dz? 


For axially symmetrical flow, the system [2.2] becomes a 
system of ordinary differential equations 


afi(r) 1 dfilr) 


— &xf3 = 2.3 
dr? r dr 2-3} 
1 — 8kfi(r)f(r) — 4(k + LfX(r) = 0 [2.4] 
dr? r dr 
FIAT) + + = 2.5] 
dr? r dr 


where r = y? + z*. Equation [2.3] is satisfied by 
1 2k r. 

Let us consider, however, a solution that is continuous on 
the symmetry axis. In general form, such a solution of 
Equation [2.3] is 

fi = Ap(Ar) [2.6] 
where A is an arbitrary constant and p(r) is a particular solu- 
tion. Let, for the sake of being definite, p(0) = 1. Obviously 
p’(0) = 0. The function f;(r) is simultaneously the particu- 
lar solution of the homogeneous equation corresponding to 
[2.4]. From this we have 


Hn = {B+ (frase) ae 
0 0 


where B is an arbitrary constant. For [2.5] we have 


far) = 2 f e+ + rar [2.8] 
0 0 
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The values for p(r) and p’(r) are best obtained by numerical 
interpolation of the equation with respect to the function 


kK(r) = fi-(r) [2.9] 
The function p(r) is approximated quite closely by the func- 


tion sec (2.2335 r). We obtain the wall of the nozzle by cal- 
culating the streamline 


{2.10} 


The sonic line passes through the origin for B = 0. 

Calculations for B = 0 and A = 1 show that short 
transonic and subsonic sections have nozzles with ordinates 
y>02atz=0. 

On the symmetry axis we have in this case \ = 1 + 2z. 
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1 Yur’ev,I. M., “On Nozzle Design,” Prikladnaya Mat 
(Applied Math. and Mechanics), vol. 19, no. 1, 1955. 
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Reviewer’s Comment 


It is the opinion of the reviewer that the inclusion of the 
second-order term in Equation [1.4] is not necessary in 
considering nozzle flows in the throat region. However, in 
treating flow problems containing stagnation points, it would 
be necessary to use an equation such as Equation [1.4] 
for proper representation of the function F. Tomotika and 
Tamada (1) also have treated the problem, including a 


second-order term, for similar reasons, but the method given 
by Yur’ev here lends itself to easy numerical calculations. 

G. V. R. Rao 

Rocketdyne Division 

North American Aviation, Inc 


1 Tomotika, S. and Tamada, K., ‘‘Studies on Two-Dimensional Transonic 
Flows of Compressible Fluid—Part III,’”’ Quart. Appl. Math., vol. IX, no. 2. 
July 1951, pp. 129-147. 


Calculation of the Heating of Shells 
at High Speeds I. N. SOKOLOVA 


The shell of an aircraft may become highly heated when flying at high speed. Thin metallic 
shells, thanks to the high heat conduction of the metal, rapidly heated throughout their thickness 
to an equal temperature. In nonmetallic shells, even thin ones, the heating is much slower than 
in metallic ones, and the temperature of the internal surface may be considerably different from that 
of the outer surface. The solution of the problem of heating such shells under nonstationary con- 
ditions involves considerable difficulty, since the connection between the heat flux inside the wall 
and the surface temperature must be determined from the heat conduction equation, an exact 
solution of which is difficult to obtain for arbitrary boundary conditions. An attempt is given to 
obtain an approximate solution of such a problem. 


1 Statement of the Problem 


E CONSIDER the heating of a shell placed in a stream 

of gas of variable velocity, pressure and density. The 
inner surface of the shell is assumed to be thermally insulated. 
It is assumed that the heat flow along the surface can be neg- 
lected compared with the transverse flow, a condition that 
is satisfied when 


6*(d In kn/dr) 1 


where 
6* = the thickness of the shell 
kn = heat transfer coefficient 
x = coordinate along the shell 


T ‘ranslated from [zvestiia Akademii Nauk SSSR, Otdelenie Tekh- 
nicheskikh Nauk, Energetika i Avtomatica (Bull. USSR Acad. Sci., 
Div. Tech. Sci., Power and Automation), no. 3, 1959, pp. 90-94. 
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In this case the problem becomes one-dimensional. We 
choose the origin on the external surface of the shell and direct 
the only spatial axis y from the outer into the inner surface. 
The temperature distribution over the thickness of the shell 
is determined by the heat conduction equation 


27 
[1.1] 
ot oy? 
where 
T = temperature 
t = time 
k* = coefficient of heat conduction of the material 


On the outer surface we have the condition of heat balance 


wep( T° — + — 
— To’) = [1.2] 
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where 


Po, Uw, To = density, velocity and temperature in the in- 
coming stream 


Cp = specific heat at constant pressure 

= adiabatic wall temperature 

T; = temperature of the outer surface 

A* = coefficient of heat conduction of the material of 
the plate 

o Stefan-Boltzmann coefficient 


emissivity coefficient 


On the inner surface the condition of heat insulation is 
satisfied 


(OT /dy)y-5* = 0 [1.3] 
The initial condition is 


at t=0 [1.4] 


The flow parameters M.., p~, p~ and are arbitrary func- 
tions of time. The heat transfer coefficient ka depends on the 
numbers M,, and R.., and on the structure of the boundary 


layer. 


2 The Method of Parabolas 


Since an exact solution of the equation of heat conduction, 
satisfying the foregoing boundary conditions, is very difficult 
to find, we shall try an approximate solution. The approxi- 
mation will consist of replacing the actual temperature pro- 
file along the thickness of the wall by a parabolic profile 


T = + ((T1 — T2)/6**\(6* — y)? [2.1] 


(Tis the temperature of the internal surface), and of solving 
the equation of heat conduction in the mean. We integrate 
the equation of heat conduction over the thickness of the wall, 
substituting Equation [2.1] for the temperature, and obtain 
a tsa [2.2] 
In the heat balance equation, for the sake of simplicity, we 
linearize the radiated heat flux 


Tit — Tot = — To) (2.3) 


which makes for a somewhat excessive value for the sought 
temperature. After this simplification the equation of heat 
balance becomes 


knpoU — T;) — (2A*/6*)(T; — T2) = 
4e0T..7, — To) [2.4] 


We introduce 


A = 3k*/6* 
V = kapoUwcp + 24*/6* + 4eoT 
M = —2d*/6* 


Q = 4eoT + knpoU 


With this notation, Equations [2.2 and 2.3] become 


d 1 
MT.+ NT, = Q a (7. + = A(T; — T2) [2.6] 


Introducing 
z2=T7.+(1/2)7, [2.7] 


the system [2.6] can be reduced to a single equation 


dz M+N (3/2)Q 


(1/2)M —N* (1/2)M—WN 


the solution of which is found in quadratures 


| 3 Qat | 
z = F(t) | constant — =~ A 
2 F(t)((1/2)M — N] 


(r ( 4 w) [2.8] 


After determining z the quantities 7; and T; are calculate: 
from 


_ _Mz-Q (1/2)Q — Nz 
1 (1/2)M — N (1/2)M — N 


7 


The constant in the expression for z is determined at a cer- 
tain value ¢ = t, different from 0. Near the start of motion, 
ie., before the instant t; is reached, it is necessary to solve the 
problem by a somewhat different method. The point is that 
the profiles of the temperature at the initial period of motion 
are poorly approximated by a parabola extending over the 
entire thickness of the shell, since actually the wall is heate 
in a short time to a smaller depth than that given by a par:- 
bolic profile. Therefore in the first instant of motion the 
temperature must be determined by a more accurate methoi. 


3 Method of Incomplete Parabolas 


We shall approximate the temperature profile by means of 
a parabola not over the entire thickness of the wall 6*, but 
only over a thickness 6(¢), on the inner boundary of whicli 
y = 6(t) the condition 07 /dy = 0 is satisfied. 

Thus, we assume that the heat wave is propagated not in- 
stantaneously, but at a corresponding thickness 6(¢) at each 
instant. Consequently, the temperature profile will have the 
form 


T = T) + (Ti — 16(t) — [3.1] 
As before, the equation of heat conduction will be solved in 
the mean. Integrating the equation of heat conduction from 
0 to 6(é) and inserting the expression for 7’ and for (O7'/Oy),—0, 
we get 
d 6k* 
— [6(t)(7;, — T. = — (7, — T, 3.2 
di 1 0) [3.2] 
subject to initial conditions [1.4]. 
In this equation the unknowns are 6(t) and T;.. They are 
determined by simultaneous solution of this equation and 
the equation of heat balance. We introduce the notation 


— To) = S 3.3] 

Then the system of Equations {2.2 and 2.4] assumes the 
form 

dS 6k* 


— =-—§ 


g(t)? — [knpoU wep + — 


2\*S =0 [3.4] 
where 
= knpoU T — To) + 4€0T — To) [8.5] 


From a solution of this system we determine S and 6(t), and 
T; is determined from Equation [3.3]. The solution is carried 
out by this method up to the instant é, at which 6(¢) becomes 
equal to the thickness of the wall 6*. 

The constant in the form [2.8] is determined at t = t; from 
the value 


constant = + (1/2)7,(t) 13.6] 


4 Estimate of the Accuracy of the Approximate 
Solution 


To estimate the accuracy of the foregoing approximate 
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solution for the heat conduction equation, let us compare the 
exact and approximate solutions obtained for particular 
specified boundary conditions. 

a) The temperature of the external surface is specified in the 
form 


T; =at+b [4.1] 


The exact solution of the heat conduction equation, subject 
to such a choice of 7), has the form 


1606** 1 
9 3 
0 (2m + 1) 
—k*(2m + (2m + 


The approximate solution for the same problem has the form 


5 
T = + (Sake — for t< [43] 


exp |-# 


3k* — y)? 6*? 
exp (- (t — for t> [4.4] 


Figs. 1 and 2 show a comparison of the temperature profiles, 
calculated from the exact (curve 1) and approximate (curve 2, 
approximation by one parabola) equations for plates 30-mm 
thick, at values of k* equal to 0.12 cm? per sec and 0.008 cm? 
per sec, and for the different values of ¢ as indicated on these 
curves. The maximum discrepancy between the exact and 
approximate solutions amounts to approximately 7 per cent. 

b) The temperature gradient is specified 


or 
Oy | y= 


T=T7T,- 


f(t) [4.5] 


In this case the exact solution of the heat conduction equation 
has the form 


T(y, t) = To 
ApRIL 1960 
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For the sake of simplification, we put f(t) = constant = B. 
Then the exact temperature 7’ = T'(y, t), after improving the 
convergence, becomes 


y k*Bt 
r= (1-32) - 


1 
exw (- 1) cos [4.7] 


The corresponding approximate solution will be 


B = 
T(y, t) = To — ov — y]? 


6*2 
fort < [4.8] 
k*B B 
T(y, t) = To (t — t;) age — y)? 
6*2 
for t > [4.9] 


Figs. 3 and 4 show a comparison of the exact and approxi- 
mate solutions for the foregoing two plates: Curves 1 are for 
the exact solution, curves 2 are for approximation by one 
parabola. The error in this case does not exceed 4 per cent. 

From a comparison of the exact and approximate solutions 
for particular types of boundary conditions, it can be as- 
sumed that the accuracy of the approximate method is of the 
same order of magnitude even for arbitrary boundary condi- 
tions. 

Greater accuracy can be attained by approximating the 
temperature profile by two parabolas, each satisfying, in the 
mean, the equation of heat conductions at half the thickness 
6(t). We choose the parabolas in the form 


T® =a +bytay? for y < (1/2)6(t) [4.10] 


T® =a, +b (v 0) +e (y - for 


5 <y <6) [4.11] 


The coefficients of the parabolas are determined from the 
boundary conditions and from the condition of the joining 
of the parabolas—the equality of the temperatures and their 
derivatives at y = (1/2)6(t). 

Let us use the previous examples to follow the variation in 
the accuracy in the case when the temperature profile is ap- 
proximated by two parabolas. 

Casea) T, = at 

By solving the equations of heat conductions in the mean, 
we obtain at ¢ > 6*?/8k* 

= (1/8)at &(t) = [4.12] 
where 9 is value of the temperature on the line where the 
parabolas join. 


The temperature profiles are determined from the expres- 
sions 


3 5 y* 


r= 
= 
2 V 8k*t 


The unknown 0 and T; at t > 6*?/8k* are found by solving 
the heat conduction equation in the mean; we shall not give 
these expressions here because of their complexity. The ex- 
pressions for the temperature profiles will be 


T = T, + (80 — 47, — 5 
y? 


(47; + 87: — 120) 


[4.14] 


y — (1/2)6* 


T® = 0+ — 40) 6* 
[y — (1/2)6*)? 
(490 — 47.) 
Case b) 
ot 
4.15 
[4.15] 


The approximate solution at t < 6**/9.6k* vields in this 
case 


= V9.6kt © = —B&(t)/16.5 [4.16] 


The temperature profiles have the form 
= —Bé(t) | 0.3712 — =~ 4+ 0.7576 —— 
ol 


T® = —0.2424B6(t)[1 — y/8(t)]? 


When t > 6*?/9.6 


TO = (20 —T.—- + By + 


1 B 
4T, — 40) —— — | 7? 
— (1/2)a(¢ 
= + (47; — — 
(y — (1/2)a(t))? 
49 — 4T,) 
6(t)? 
[4.18] 


6 and 7.2 are determined when ¢ > 6*2/9.6 from the solu- 
tion of the equations of heat conduction in the mean. 

The temperatures, calculated from these formulas, are also 
shown in Figs. 1, 2, 3 and 4 (curves 3). The accuracy is con- 
siderably increased when approximation with two parabolas 
is used. 


Original submitted Aug. 18, 1958 


Reviewer’s Comment 


This paper represents an interesting application of an 
approximate technique of solving partial differential equa- 
tions known as the integral method. In this type of approach 
the equation is satisfied in the mean, the mean being taken 
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over a certain thickness variously called the boundary layer 
thickness or the thermal thickness, depending on the physical 
problem considered. 

Von Karman (1) and Pohlhausen (2) originally developed 
the method. It was extended by Gruschwitz (8), and Wieg- 
hardt (4) applied it to the solution of the energy equation. 
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Most recently, it has been utilized in a series of papers by 
Goodman (5) in treating Stefan’s problem (heat transfer at a 
moving boundary). 

In the present paper, Sokolova has treated heat transfer 
with either a stationary or moving inner boundary. He finds 
that the use of two polynomials leads to greater accuracy than 
the use of a single polynomial (parabola). It would have 
also been interesting for the author to compare his two-para- 
bola approach with the use of an exponential curve which also 
satisfies the boundary conditions. 

Another extension of the present work, which would be of 
great interest, would be the specification of the environmental 
conditions pe and 7’.., and the flight speed U.. as a function of 


time, rather than the boundary conditions [4.1] or [4.5] se- 
lected by the author. 

—Srincuaire M. Scana 

General Electric Co. 


1 von Karman, T., ‘‘Uber laminare und turbulente Reibung,”’ Z. Ange- 
wandte Math. und Mech., vol. 1, 1921, translated as NACA TM 1092. 

2? Pohlhausen, J., ““‘Zur naherungsweisen Integration der Differential 
gleichungen der laminaren Grenzschicht,” Z. Angewandte Math. und Mech., 
vol. 1, 1921, p. 252. 

3 Gruschwitz, E., ‘Die turbulente Reibungsschicht in ebener Strémung 
die Druckabfall und Druckansteig,”’ Ing. Arch., vol. 2, 1931, p. 321. 

4 Wieghardt, K., ‘‘Uber einen Energiesatz zur Berechnung laminaren 
Grenzschichten,” Ing. Arch., vol. 16, 1948, p. 231. 

5 Goodman, T. R., ‘“‘The Heat Balance Integral and its Application to 
Problems Involving a Change of Phase,” Trans. ASME, 1958, p. 335. 


Analysis of Operation of Simplest 
Ramjet Combustion Chamber Under 


Flying Conditions 


HE PARAMETERS at the inlet of a ramjet combustion 

chamber are subject to considerable variation in flight. 
The variations will be particularly strong for launching-type 
engines. The pressure and the temperature at the chamber 
inlet depend on the altitude and velocity of flight. In ad- 
dition, under different control methods, the rate of flow into 
the chamber and the composition of the mixture do not re- 
main constant. All this naturally influences the operation of 
the chamber and the heat release ratio. 

In this paper we analyze the influence of flight conditions 
on the required dimensions of the combustion chamber. 
Such an analysis is made possible by the existence of equa- 
tions that relate the principal combustion characteristics with 
the pressure, temperature, mixture composition and flow tur- 
bulence. These relations have been investigated essentially 
for combustion in a turbulent homogeneous mixture, in which 
the evaporation of the fuel and the mixing of the fuel vapor 
with the air are fully completed by the instant combustion 
begins. Such conditions occur at high temperatures at the 
chamber inlet, i.e., at high flight velocities (17 > 3). How- 
ever, even at lower temperatures (lower flight velocities), 
when the fuel does not have time to evaporate completely, it 
is also possible to use as a first approximation the equations 
for a homogeneous mixture, since they are quite close to the 
relations that hold for an inhomogeneous mixture. 

Let us first establish the variation of the speed of flame 
propagation in a turbulent stream wu, and of the combustion 
time ¢, as functions of the pressure and temperature. Al- 
though the dependence of the speed of flame propagation and 
of the combustion time on the pulsational speed w’ and on the 
normal speed u, have by now been sufficiently wel! investi- 


Translated from Jzvestiia Vysshych Uchebnykh Zavadenii 
MVO, Aviationnaya Tekhnika (Bull. Higher Institute of Learning, 
Aviation Engng.), no. 2, 1959, pp. 122-133. 
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A. V. Talantov 


Kazan’ Aviation Institute 


gated, the influence of pressure and temperature has not been 
so well studied. We shall attempt to establish the influence 
of these factors, using the customary relations for the speed wu, 
and the time t. 

In accordance with the theory of “‘surface’’ combustion and 
with experiment (1,2), the speed of flame propagation in a 
turbulently flowing homogeneous mixture can be written 


= {1] 
where 


AwiaA 
a = 0.6 to 0.8 


There are no indications in the literature concerning the de- 
pendence of the pulsational speed on the stream temperature 
(4). The normal speed is proportional to the 1.8 power of 
the temperature (3). In the range of moderate ratios of pulsa- 
tional to normal speed, which is the range of practical interest 
(but at w’/u, > 1), we can assume a © 0.8. Then 


~ T%4 [2] 


This result is in satisfactory agreement with experiment (1,4). 
The pressure dependence of u, computed from Equation 
[1] is less than that found from experiment (4,5). 
In accordance with (3 and 4), we can assume 


Un w! ~ 
Then the pressure dependence on the speed of flame propaga- 
tion in turbulent flow becomes 


1 Numbers in parentheses indicate References at end of paper. 
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fr 
: 
O 
ul 
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On the other hand, experiment yields (4,5) 
ps 


One cause for this disagreement is the damping of the tur- 
bulence behind the grids (4). As the pressure is reduced and 
the speed u, is decreased, the flame becomes longer. Conse- 
quently the average value of the pulsational speed, which de- 
termines u;, decreases, and this causes an increase in the ex- 
ponent for p. In practice we frequently produce forced turbu- 
lization of the stream, and we shall therefore assume in our 
calculations the experimental value 


[3] 


We estimate the combustion time from the following rela- 
tion, which is theoretically derived (2) 


[4] 
w Un 
Here h is the scale of turbulence. 


This relation is in satisfactory agreement with the experi- 
mental data on the effects of lb, w’ and up, at B = 3. 


t/t, 


05 


0 T/Te 


Fig. 1 
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If it is assumed that the pulsational speed and the scale of 
turbulence are independent of the temperature (4), then the 
dependence of the combustion time on the temperature, with 
allowance for the influence of the temperature on the normal 


speed, can be written 
yf 1.8 
In E (7) ] 


(+2) 
Inf1l+— 
Ung 

We shall compare this theoretical equation with the only 
experimental data known on this subject, namely those of 
Doroshenko and Nikitskii (4). 

The points on Fig. 1 show the dimensionless values of the 
combustion time vs. the temperature. The curves represent 
the values calculated from Equation [5]. The solid line is for 
a pressure of 200 mm mercury, the dotted one for 600 mm 
mercury. 

Thus, theory and experiment show that the combustion 
time decreases with increasing temperature. It is possible 
to use Equation [5] in approximate calculations. 

In examining the pressure dependence of the combustion 
time it is necessary to take into account not only the influence 
of the pressure on the normal and pulsational speeds, but also 
on the turbulence scale (4) 


t = [5] 


L~ 


Taking all this into account, we can write for the combus- 


tion time 
0.5 
ann [t+ 
Un/Jo 


The ratio of the combustion time to the initial time, calcu- 
lated from Equation [6], as well as the results of the experi- 
ment of Doroshenko and Nikitskii (4), are plotted against the 
pressure ratio in Fig. 2. As in the case of the temperature de- 
pendence of the time, the curves differ little at various ratios 
of w’/u, The solid line on the curve corresponds to a tem- 
perature of 373 deg abs, the dotted line is for 573 deg abs. 

It is seen that both theory and experiment, which are in 
satisfactory agreement, evidence a rather intense increase in 
combustion time with decreasing pressure. We shall use 
Equation [6] in further calculations. 

When the pressure and temperature change simultaneously 
with the mixture composition and the flow rate, the equation 
for the combustion time becomes 


Po w 
w’ 
Un Jo 


Here, as before, the subscript zero pertains to initial values, 
for which the initial time & was determined, and un, and 


Un, are the normal speeds at the current and initial values of 


a. This equation is correct only when the effect of each 
factor is independent of the remaining parameters. 

The foregoing equations make it possible to calculate the 
length that the combustion chamber must have in order for 
the process to be completed in it. 

To calculate the first component X, of the chamber 
length, from the edge of the stabilizers to the place where the 
flames join, we can use a system of equations similar to that 
given in (6). To simplify the calculations we shall assume, in 
addition to the other assumptions made in that paper, that 
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the flame is infinitesimally thin. The system of equations is 
then written in the following form: 
Equation of conservation of energy for a fresh mixture 


k-1 

= 1 — —— — 1) [8] 
2 
where 

Tm, Um = dimensionless temperature and speed of fresh mix- 
ture, referred to the corresponding parameters at 

the chamber inlets 

My = Mach number at the chamber inlet 


Adiabatic equation for the fresh mixture 
Tn = [9] 


where 7 is the pressure in the section under consideration 
divided by the pressure at the chamber inlet. 

Equation of conservation of energy for the flow of combus- 
tion products 


k-1 
where 
Tp, Up = dimensionless temperature and speed of combustion 
products 
A, = 1+ q/epTo 
q = heat supplied per kg of air 


Equation of conservation of mass 
— tm (1 — 9) + — = 1 [11] 
Tin Tp 


where 7 is the dimensionless ordinate of the flame front, re- 
ferred to half the gap between stabilizers. 
Momentum equation 


1 


Equation for the fraction of heat liberated 
= [13] 


Equation of conservation of mass of fresh mixture for two 
neighboring sections, 7 — 1 and 7, in converted form 


point of joining of the flame fronts to the completion of com- 
bustion, can be found by integrating the equation of motion 
of the mixture in the zone 


dX = w,-dt 


Here w, is the rate of displacement of the mixture in the zone. 

This rate is usually variable and a function of the heat 
liberation and—the burnup (r) in the zone. An exact calcu- 
lation of X is possible by means of a procedure detailed in (6). 
Taking into account the fact that in this case the main varia- 
tion of the parameters is in the axial direction, and neglecting 
the change in the parameters in the transverse direction, we 
can determine X approximately by considering the problem 
as one-dimensional. The rate of flow in this zone can be 
written in this case as 


+ — 


= 


where 6 = T/T» is the ratio of the combustion temperature 
to the temperature at the entrance to the zone. This equation 
for the rate takes into account only the density changes due 
to combustion; the compressibility is disregarded and the 
formulation of the problem is hydrodynamic. 

The distance X becomes 


X= wo [1 + r(@ — 1)]dt 


Here r is a function of time. The form of this function is 


given in (2) as 


un\*} 1 
r=3 


(1 + w’/un)? 


where t, = l/w’ is the characteristic time. 
integrating we get 


1 | [16] 


Inserting r and 


lin 
X = wol + — 1) (1 X 
w 


AG + = ;) x 
3 1 + w’/Un 


t 1 1 1 
)] 


—(4/Tm)i(1 — Ni 


AX = — Ni-1)Umi—-1 
+ (4 /Tm)i] the 


Here 


fi, = w,/wo = dimensionless speed of flame propagation, re- 
ferred to the speed at the chamber inlet 

dimensionless distance along the axis between neighbor- 
ing sections 7 — 1 and ?, referred to half the gap be- 
tween stabilizers 


> 


It is also possible to use the approximate relation 
Tp/Tm = Xp [15] 


These equations are sufficient to determine the unknown 
values of 7, Um, Up, Tm, Tp, 0, rand AX, one of which serves as 
the independent variable. 

The calculation is best carried out in the following se- 
quence. We first find the dependence of the dimensionless 
ordinate 7 and of the other parameters on the combustion 
coefficient r. Then we use Equation [14] to calculate the 
value of AX, which when summed yields X,;. 

The second component X, of the chamber length, from the 
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The total length of the combustion zone can be obtained 
from this expression by putting ¢ = t, 


X; = Bu + (6-1) x 


w’ 3 3[1 + 3 

1 w’ 1 

1 + (w’/un) +2) + +( 


1 w 1 
[1 + (w’/un)]? [1 + (w’/un)] i}} 


Experiment yields that B is equal to 3. Knowledge of the 
path X covered by the hot mixture (Eq. [17]) and the 
knowledge of the burnup r (Eq. [16]), as functions of the 


time make it also possible to calculate the heat liberated along 
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the length, r = f(X). The burnup is found to be quite slug- 
gish in the final portion, where a change in the coefficient 
of heat liberation r from 0.95 to 1 requires approximately the 
same distance as a change of r from 0 to 0.95. In real cham- 
bers it is therefore useless to strive for complete heat liberation 
in the final portion at the expense of considerably lengthening 
the chamber, for in the final analysis, the losses produced 
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thereby will offset any gain from the increased heat liberation. 

We used the foregoing procedure and the derived equations 
to calculate the combustion process in ramjet engine cham- 
bers. 

One of the objects of calculation was a subsonic ramjet en- 
gine with nonadjustable transition sections, operating at alti- 
tudes from 0 to 8000 m in the range of M = 0.4 to 0.8. The 
characteristic distance | between stabilizers was taken to be 40 
mm, the degree of turbulence « = 0.1, the turbulence scale 
lb = 0.50em. The fuel was gasoline and the mixture composi- 
tion stoichiometric. The parameters of the mixture entering 
into the combustion chamber were calculated by the proce- 
dure customarily used in engine theory. 

Fig. 3 shows the required Jengths X ; and X, as functions of 
the rate of flow (Mach number) for 0, 5000 and 8000 m. 

An examination of the results obtained makes it possible to 
draw the following conclusions. 

The principal portion of the chamber, for a chosen size o/ 
gap between stabilizers, is the combustion zone. The actual 
distance from the section through the stabilizers to the point 
of joining of the flame fronts is less than 10 per cent of the 
total length of the flame. Consequently, in those cases when 
no great accuracy is required, this quantity can be neglected 
to simplify the calculations. It is appropriate to emphasize 
the need for a thorough investigation of the combustion tin: 
in turbulent stream, in contrast with the heretofore preferred 
investigation of the rate of flame propagation. 

In accordance with the preceding, the dependence of the 
total length of the combustion’ chamber on the altitude and on 
the flying speed are determined essentially by the combustion 
time. As the flight speed is increased, the temperatures and 
the pressure at the inlet to the chamber increase, and this 
should lead to a reduction in the required length. This does 
not happen, however, owing to the increase in speed of flow 
through the engine, including the chamber inlet, with increas- 
ing Mach number. 

The required chamber dimensions are greatly influenced by 
the change in flying altitude. As this altitude increases the 
combustion conditions are substantially worse; the length of 
the chamber for complete combustion increases. 

It is easy to determine the coefficient of completeness of 
heat liberation in the chamber if its actual length is specified 
and the variation of burnup along the length of the chambe: 
is known. Fig. 4 shows the dependence of completeness of 
heat liberation 2 for the engine under consideration, for a 
chamber 1.3m long. The resultant relations are in satisfac- 
tory agreement with the actual variation of 7, with altitude 
and with Mach number, 

We also used the procedure developed here to calculate the 
dimensions of a combustion chamber for engines operating at 
higher speeds, M = 1 to 4, and at altitudes from 0 to 20,000 
m. 

It must be borne in mind that these results were obtained 
not for one particular engine operating under different con- 
ditions, but for a series of engines, each designed for its own 
operating mode. 

The coefficient of speed at the chamber inlet was assumed 
to be 0.2 for all modes with the exception of M = 1, where 
A = 0.15. The coefficient of pressure recovery o in the 
system of steps at the inlet to the engine was assumed to have 
values of 1, 0.9, 0.7 and 0.5 for M = 1, 2, 3 and 4, respec- 
tively. The composition of the mixture was assumed 
stoichiometric. The turbulence scales were 0.5 and 1 cm; 
the degree of turbulence ¢€ was 0.1. 

The results of calculations are shown in Figs. 5-8. One of 
the interesting results of the calculation is that in the region 
of M = 4, owing to the high temperatures and large pressures 
at the chamber inlet, the rate of flame propagation becomes 
equa! to the flow rate and even exceeds it. These results are 
reliable only to the extent that the foregoing relations, ob- 
tained for small pressures and moderate temperatures, can be 
extended to include large pressures and high temperatures. 
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A direct experimental verification of these relations in this 
region has not yet been made. It is possible that the velocity 
u, and the combustion time ¢ wil] not obey the foregoing rela- 
tions because of processes that may occur in the mixture of 
high temperatures prior to ignition. Equally valid is the 
assumption that owing to the high temperatures the principal 
process whereby heat is supplied to the chamber under these 
conditions (1/> 4) is self-ignition, rather than flame propaga- 
tion. However, with all this, one can expect the parameters 
of the combustion chamber to be of the same order of magni- 
tude, and in any case not worse than those obtained in the 
foregoing calculations. 

At M = 1 or 2, the size of the initial section X, is relatively 
small. Thus, for supersonic engines, as for subsonic ones, it 
is possible to neglect in approximate calculations the length 
of the initial section (before joining of the flame fronts). 

The value of the combustion time and the corresponding 
required chamber length for a supersonic engine depend very 
strongly on the flying altitude and speed. The length of the 
chamber decreases with increasing Mach number. The effect 
of the flying speed is explained by the increased temperature 
and pressure at the inlet to the chamber, which lead to a re- 
duction in the combustion time. A certain increase in the 
speed of the stream entering the chamber at constant X, 
contributing to a lengthening of the zone, does not exert a 
decisive influence. However, the dependence of the required 
Jength of the combustion chamber on the flying speed could 
actually be somewhat weaker. The point is that, at a con- 
stant amount of heat delivered per kilogram of mixture, the 
relative heat delivery decreases with increasing flying speed 
and consequently with increasing stagnation temperature. By 
virtue of this, the maximum possible speed coefficient at the 
chamber inlet increases, and, if the maximum is reached, this 
leads to a lengthening of the zone (the other parameters re- 
maining constant) and to a slowing down in the rate of de- 
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crease of chamber length with increasing flying speed. In 
the subsonic region, an increase in flying speed leads to a slight 
increase in the stagnation temperature and pressure; the in- 
crease in the speed of the stream entering the chamber exerts 
the decisive influence, and the chamber length increases with 
increasing flying speed. 

As the flying altitude is increased, the required chamber 
dimensions increase because of the drop in temperature and 
pressure, 

From an examination of the curves of Figs. 5 and 6 it fol- 
lows that to design a ramjet engine that operates satisfactorily 
in the subsonic region (M < 1) at altitudes greater than 10 
km, serious difficulties are encountered from the point of view 
of effecting combustion with ordinary fuels. The required 
length of the combustion chamber under these conditions will 
be measured in meters, which will lead to large energy losses 
and to an excessive increase in the structure weight. 

It must be borne in mind that the foregoing required cham- 
ber lengths are calculated for complete heat liberation, in- 
cluding the final stagnant section of the burnup. The tur- 
bulence seale in these calculations, equal to 0.5 or 1 cm, is not 
the minimum for real chambers. Consequently, the length ob- 
tained must be considered only as the upper limit. 

It is interesting to determine the change in the necessary 
longitudinal dimensions of a combustion chamber for engines 
having an identical pressure at the chamber inlet. The point 
is that for identical specific rated wing load in horizontal 
flight, the pressure in a ramjet combustion chamber will have 
approximately the same value for different altitudes and flying 
speeds. Therefore, for real engines, intended to fly under 
different conditions, the pressure in the chamber still changes 
insignificantly. The dotted curve of Fig. 5 is drawn through 
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the points for chamber pressures of 1, 2 and 5 atm. It is seen 
that the necessary chamber dimensions change little at con- 
stant pressure. As the altitude is increased and as the flying 
speed is increased, the condition that the pressure be main- 
tained constant necessitates some reduction in the dimensions 
of the chamber. This is explained by the increase in tem- 
perature at the combustion chamber inlet, leading to a reduc- 
tion in the combustion time, in spite of the increase in flow 
speed w. If one assumes w ~ VT (at constant A), and if the 
dependence of time on temperature is taken in accordance 
with Equation [5], then the length of the combustion zone, 
proportional to the product of the velocity and the time, will 
obey the following law 
1.3 
X,~ wt ~ inf + (7) ] 


Un 


It is easy to verify that as the temperature increases Xz 
decreases. The dependence of the total heat release rate q 
and of the space heating rate per atmosphere gp on altitude 
and flying speed are shown in Figs.7 and 8. What is striking 
is the exceedingly strong dependence of g and q, on the flying 
conditions, their decrease with increasing altitude, and their 
increase with increasing flight speed. The dependence of the 
total heat release rate on the flying speed may be even 
stronger, if the increase of X at the chamber inlet, noted 
previously, takes place. An increase in the flow rate at the 
chamber inlet leads to a proportional increase in the heat 
delivered to the chamber; lengthening of the zone is not de- 
cisive in this case. The total heat release rate q under cer- 
tain conditions reaches hundreds of millions and billions of cal/ 
m*-hr, which becomes comparable to the heat release rate 
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in liquid fuel engines. The space heating rate per atmosphere 
reaches hundreds of millions of calories. 

What is interesting is not only the values of the mean 
space heating rate but also the space heating rate per atmo- 
sphere change with the flying conditions. This indicates 
that the method used to design chambers for an average value 
of q is defective, even if one takes for g its value per at- 
mosphere. 

The analysis performed makes possible several recom- 
mendations toward improving the parameters of the process 
in the chamber under fixed flight conditions. 

In accordance with the foregoing, we propose that the fac- 
tor determining the chamber length is the combustion time, 
which obeys Equation [4]. The value of t is affected by the 
normal and pulsational speeds and by the turbulence scale. 
For chambers operating with stoichiometric mixtures of ordi- 
nary hydrocarbon fuels, such as gasoline, kerosene and others, 
there is no possibility of forced operation by increasing the 
normal speed. At the same time, one cannot exclude the 
use of new fuels, which yield substantially greater values of 
normal speed, or the use of additives to the fuel so as to in- 
crease the normal speed. 

In the case when the speed coefficient \ in the chamber is 
quite different from the maximum possible one, an increase 
in this coefficient makes it possible to increase the calorific 
value somewhat. A certain possibility of forcing is contained 
also in an increase of the pulsational speed by insertion of 
specially selected turbulizing grids. 

One opportunity, so far not used, for forcing the chamber is 
to reduce the turbulence scale. This can be done by inserting 
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It is desirable to install the grids 
as close as possible to the combustion zone, so as to avoid 
damping of the turbulence with increasing distance from the 
grids, i.e., an increase in the scale and a decrease in the pulsa- 
tional speed. 


suitable grids of thin rods. 
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Reviewer’s Comment 


The results of Talantov’s calculations are qualitatively cor- 
rect, although the degrees of some of the effects noted in the 
figures are questionable with respect to real engines. For 
example, in the flight speed range between Mach 3 and Mach 
4, the engine will be getting into the spontaneous ignition 
regime, particularly at the lower altitudes where stagnation 
pressures and temperatures are higher for a given flight speed. 

For altitudes above the isothermal range, i.e., above 26,000; 
M, temperatures again rise.) In such engines, it is a poor 
assumption that the comLustion chamber will receive a homo- 
geneous mixture, or even a good distribution of liquid fuel at 
the combustor inlet, because the fuel injectors should be placed 
close to the flameholders to prevent premature burning. 
Significant lengths, relative to the X,’s shown, will be required 
simply for fuel distribution and mixing of burned and partially 
burned parcels of gas in order to obtain reasonable engine 
efficiencies. Hence the engine lengths shown for the Mach 3 
to 4 range in Fig. 5 and 6 are inordinately small, and the 
statement that ‘“...the (combustion chamber) length ob- 
tained must be considered only as the upper limit” is incorrect 
for these high speed engines. 

With respect to the theoretical bases for the analysis, there 
are questions with regard to the assumptions employed for 
both the flame spreading length and the subsequent burnup 
length. The effects of temperature and pressure on turbulent 
flame velocity used by Talantov are presumably based on 
experiments with open or burner-type turbulent flames. Al- 
though the reviewer did not have access to the references 
cited by Talantov, the relationships for the effect of pressure 
on U, quoted are similar to those presented by Goldenberg 
and Pelevin (1) and Khramtsov (2), both of whom employed 
open flames. On the other hand, Scurlock (3) found little 
effect of approach stream turbulence on spreading rate of 
confined turbulent flames such as we must consider in a ram- 
jet combustor. Scurlock and others have shown that a large 
part of the turbulence in confined flames is due to shear be- 
tween burned and unburned gas flows in two regions: Iin- 
mediately behind the flameholder the unburned mixture moves 
faster than the recirculating burned gas, and downstream of 
this piloting recirculation zone the situation is reversed as the 
burned gas accelerates. Further, Markstein (4) points out 
that for fames in ducts the “noise” mode of generation of 
random disturbances (due to pressure oscillations) may have a 
larger effect than shear flow disturbances. Thus for ram- 
jet flames there is little connection between the local turbu- 
lence scale and intensity in either burning region and the 
turbulence parameters of the combustor inlet flow. 

Another point with regard to flame spreading is that ex- 
periments have shown that the flame spreading begins toward 
the rear end of the recirculation zone, the length of which 
could be of the same order as the calculated flame spreading 
length when practical sized flameholders are used. This 
would be important for the higher flight speeds; the tempera- 
ture ratio A, which influences the shear generated turbulence 
mentioned in the foregoing is smaller. 

For the second, or burnup region treated by Talantov, this 
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reviewer wonders what Talantov has used for normal burning 
velocity, U,, in Equations (16-18). Since all of the material 
reaching this region has already traversed the turbulent flame 
front, it is both hotter than 7) and diluted with reaction prod- 
ucts. For such material, U,, is probably 2 to 6 times that for 
the unburned mixture at room conditions (5). The same 
reservations regarding suitable values of turbulence charac- 
teristics for this region apply as noted in the foregoing for the 
flame spreading region. 

Spalding (6) has proposed that, for the high flow velocities 
of practical interest in ramjet design, the length required for 
flame spreading in a duct may be determined largely by en- 
trainment of unburned gas into the hot gas jet which com- 
mences behind the recirculation zone. This approach leads 
to the conclusions that an axial length of 3 to 5 duct widths 
(or flameholder spacing widths) must be allowed for flame 
spreading and that little can be done to alter this required 
length; for example, the analysis suggests that preheating 
(e.g., from higher flight speed) leads to greater lengths because 
of the reduction in i. 

Thus the reviewer believes that the results of this paper 
should be considered only for the purpose of predicting qualita- 
tive trends for ramjet length requirements, which have already 
been estabiished by development experience. Probably of 
more value are papers such as Spalding’s (6) which stimulates 
thought on new theoretical approaches to the conventional 
ramjet problem, and papers challenging the experimentalist 
or designer to conceive new higher output designs for ram- 
jets such as the one by Avery and Hart (7), which dealt with 
maximum attainable heat release rates in hypothetical com- 
bustors with instantaneous mixing of unburned and burning 
gas. 
In spite of the foregoing comments, however, it is encourag- 
ing to find papers such as Talantov’s which attempt to make 
practical use of some of the large mass of data from bunsen 
burner flame studies. 

—Gorpon L. DuGGer 
Applied Physics Laboratory 
The Johns Hopkins University 
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Puncture Problem at Cosmic Velocities a A. LAVRENT'EY 


T IS known from the theory of cumulative charges that 
the mechanism of puncturing of metallic plates by cylin- 
ders or by balls differs greatly at 3-10 km/sec from the 
mechanism at velocities up to 1000 m/sec. Two stages take 
place at high velocities: a—The ball or cylinder, penetrating 
into the obstacle, flows over the surface of the punched-out 
crater; b—inertial expansion of the cavity takes place after 
the “bullet has been annihilated.” The first part can be 
calculated with sufficient accuracy by using an ideal incom- 
pressible liquid as a model, but the calculation of the second 
part is more difficult, although the results available indicate 
that the principal difficulties can be overcome. The com- 
putational schemes yield good agreement with experiment. 

Much less investigated is the puncturing problem at 50- 
100 km/sec; the experimental difficulties at such velocities 
require particular caution as regards the main hypothesis 
made in the theory of each phenomenon. 

As far as I know, the literature contains only one paper 
devoted to this problem—that of K. P. Stanyukovich on the 
formation of the moon’s craters and on the determination of 
the momentum imparted by a falling meteorite. On the 
basis of his computation, Stanyukovich advances the hypoth- 
esis that as a ball drops, its kinetic energy is converted into 
potential energy of the gas into which it is converted by the 
impact; assuming further a few less clearly defined hypothe- 
ses, he arrives at the following conclusions: 

1 The momentum acquired by the struck body is propor- 
tional to the energy of the incident body. 

2 The volume of the crater is also proportional to the 
energy. 

In the present note I propose an incompressible-medium 
model for which the calculations can be carried out to con- 
clusion. My deductions do not agree with those of Stanyu- 
kovich. 


1 One-Dimensional Case 


We start with a consideration of the one-dimensional case, 
when the model is particularly simple. We consider impact 
between a plate of thickness a, flying at a velocity vo, against 
the end of a cylinder (rod) of length 1. The thickness a will 
be assumed much smaller than 1. The problem consists of 
determining the momentum which the rod acquires as a 
result of the impact. 

We assume the striking plate to be incompressible and 
absolutely hard. The rod is considered as a limiting case of 
a set of absolutely hard infinitesimally thin platelets, placed 
infinitely close to each other. During an inelastic impact 
between the striking plate and the first platelet of the rod, 
the momentum is conserved and kinetic energy will be lost as 
a result of the increase of the mass. This will continue as 
each succeeding platelet goes into motion. Let us now ecalcu- 
late, for the limiting case, the distribution of the kinetic energy 
loss in the system along the rod. 

Thus, let a rod of length z be set into motion and let v be 
the velocity of this piece; as a plate of thickness dz is set 
into motion, the change in velocity dv should satisfy the 
relation 


adv + vdx = 0 


Translated from Iskusstvennye Sputniki Zemli (‘Artificial Earth 
Satellites’), No.3, USSR Acad. Sci. Press, Moscow, 1959, pp. 61-65. 
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from which we obtain after integration, noting that v = vp when 
zr=a 
v = v(a/z) 
Let us calculate now the unknown kinetic energy loss 
E = (1/2)zxv? 
We have 
dE = (1/2)0°dz + xvdv = —(1/2)v%dz = 


2 
—(1/2)u? (<) dz = —E,~ dx 
x 


a 
where 
Ey = —(1/2)ar? 


is the initial energy of the striking plate. Assuming that the 
entire energy lost is converted into heat, we obtain the follow- 
ing distribution of heat along the rod 


dE a 
dx x? 


Let us now denote by 7) the minimum heat density necessary 
to vaporize the rod material. From this we find that the 
impact vaporizes the portion of the rod 0 < x < 2, where the 
unknown 2 is found from the relation 


To = E,(a/2»?) 


As the vapor expands it separates from the rod, while the 
remaining portion of the rod acquires a momentum J. Let 
us calculate J, for the case of greatest interest, when 


aXm<Kl 


and the energy going into vaporization is small compared 
with the initial energy of the striking plate. We can assume 
here that the entire energy of the gas goes into kinetic energy. 
We obtain for the velocity of the gas cloud 


V=aV (a/xo)vo 


Our problem has reduced to a purely gasdynamic problem. 
We confine ourselves now to only very rough calculations in 
two cases. 


= arn? 


1 Wemake an additional assumption that all the scattered 
particles of the gas cloud acquire an identical velocity V, 
and then 


V= V 


= avy? 

obtaining finally for the momentum 
I = aV = Varn = an” 


2 The second extreme case we consider is one in which 
each layer dx of the gas cylinder scatters (along the 2 direc- 
tion) independently of the other elements. Then the gas 
element dz, located at a distance x from the end, yields a 
momentum 


Vdz = 2T dz = (dz/x) 
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From this we obtain for the unknown momentum 


I = SnVdz = V/ 2E,a log V2 avo log —= 
a 


2 Three-Dimensional Case 


Let us now perform a calculation analogous to the preceding 
one for the case when a ball strikes a hemispherical plate of 
radius R. For this computation we must schematize the 
model still further. We assume that the action of the ball 
on the plate reduces to an impact of a spherical layer on the 
hall against a spherical cavity in the plate. All the velocities 
of the points of the striking layer are directed along the radii, 
and the velocity distribution coincides with the corresponding 
distribution in an incompressible liquid: The velocity v at a 
point located at a distance r from the center of the layer is 
determined from the formula 


v = v,(a/r)? 


where a is the radius of the internal cavity of the striker; v9 
the velocity of the points of this cavity. Let ka be the radius 
of the external surface of the striker. The medium is repre- 
sented by a model similar to that considered in the linear 
formulation, namely that the medium is a set of infinitesi- 
mally thin spherical layers of an ideal incompressible liquid, 
located infinitesimally close to each other. As a result of the 
impact of the striker, its mass will start increasing as called 
for by the inelastic-impact scheme; as the mass increases, 
the momentum is conserved, and the kinetic energy of the 
system is converted into heat. 

Let us perform the calculation. We determine first the 
radial momentum during the instant when the zone of 
motion in the medium is a hemisphere of radius r; let in this 
ease the radial velocity of the internal surface of the striker 
be w = w(r), then 


= [1] 


Since the medium is incompressible, the velocity v at a point 
located at a distance x from the center will be 


v = w(a/z)? 
From this we obtain for the momentum 
r 
[= = 4rwa%(r — a) 


The condition that the momentum be constant yields 
(r — a)dw + wdr = 0 
w = C/(r — a) 


where the constant C is determined from Equation [1]. 
Thus, finally 


2) 
[2] 


w 
=o)" 


We now find the kinetic energy F under the same conditions 


r a\* a 
E — (1/2) (‘) dx = 27a5w? (: [3] 
a z r 


From this we can obtain the distribution of heat 7’ over the 
body. Then 
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dx 


za 


Fig. 2 


or, taking Equation [2] into account 


r—a a? 
r 


T = 4ra*w 
r—a 


Denoting by E> the initial energy of the striker, we obtain 


2(k — 1)k (k — 1) 
T= ( — (1/2): “| & 
r(r — a)? r r3 


If we now assume that an amount of heat 7) per unit volume 
is needed to vaporize the material, then as the ball strikes 
the body, there is formed in it a gas cavity of radius 

2k(k — 1) 


[4] 
0 


This cavity, defined by expression [4], will characterize the 
lower limit of the size of the funnel. 

Let us now calculate the momentum acquired by the body. 
Again we consider the case when 


ax p<R 


and assume that a small part of the striker energy is con- 
sumed in vaporization of the foregoing mass of funnel. De- 
noting by V the velocity of the gas cloud that escapes from 
the funnel, we obtain 


V = 3% 


(2/3)rp?V? = Ey 
2p? 


From this we obtain for the sought momentum J 


I = (2/3)mp°V = AcE, A= 1) 
3 To 

It is interesting to note that in this scheme, the momentum 
depends greatly on the dimensions of the striking ball. 

In analogy with the two-dimensional case, we can carry 
out the calculations in accordance with the second extreme 
scheme. Let us assume that each element dr of the spherical 
layer of gas scatters independently of the other elements and 
that the momentum imparted is directed along the normal 
to the spherical cavity. 
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The normal velocity V of the gas cloud due to an element dr 
of the spherical layer is determined from the energy balance 


T = 2rr*V? 


If we now denote by a the angle between the symmetry axis 
and a line joining the center of the sphere with the element 
of the gas, then the component J = dI/dr of the momentum 
along the symmetry axis due to the scattering of a ring of 
width rda, a = const, will be 


dJ = 2rr? sin a cos aVda 


ll 


5/2 
rar — 1)" 
2 
Integrating the latter expression with respect to r, along the 
thickness of the gaseous portion of the body, we obtain 
I = — 


where B is a numerical constant. 


Reviewer’s Comment 


Lavrent’ev has made a novel and useful contribution to the 
literature on theoretical prediction of hypervelocity cratering. 
He does not, however, compare his approach with that of 
Opik (1), Grimminger (3), Rostoker (4), or Bjork (5); but 
notes simply that his results are different from those of 
Stanyukovich (2). A thorough review of cratering theories 
[of more ambitious scope than Allison’s paper (6) ], together 
with the correlation with experimental results, now abun- 
dantly available in the U.S. literature, would be most welcome 
at this time. 

In the opinion of this reviewer, Bjork’s treatment (5) of the 
hypervelocity impact process is the most satisfactory, having 
required the fewest assumptions. Bjork neglects the strength 
of the materials; but working with their equations of state 
at high pressures and densities, he applies the hydrodynamic 
equations to yield a complete picture of the positions of all 
particles and shock fronts at each instant in time during the 


cratering process. Correlation of his predictions with results 
of cratering experiments is good. 

—Murray 

General Electric Co. 
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Concerning the Problem of the 
Influence of Internal Degrees of 
Freedom of Particles on the 
Transport Coefficients of a 


Multicomponent Mixture of Gases 


HE HEAT conduction coefficients of shear and bulk 

viscosity of gases can be calculated by the methods of 
kinetic theory. For monatomic gases the kinetic theory gives 
good agreement between the theoretical formulas and the 
experimental data. Subject to a suitably chosen potential of 
interaction between the molecules at ordinary temperatures, 
the theoretical results can also be applied successfully to the 
calculation of the coefficient of shear viscosity of polyatomic 
gases. However, the coefficient of heat conduction contains 


Translated from Fizicheskaya Gazodinamika (‘‘Physical Gas- 
dynamics’), USSR Acad. Sci., pp. 59-69. 
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terms which cannot be obtained from the theory of monatomic 
gases. This was first demonstrated by Eucken. For mona- 
tomic gases the bulk viscosity is connected with the shear 
viscosity by the well-known Stokes equation. This equation 
does not hold for polyatomic gases. An analysis of the in- 
fluence of the internal degrees of freedom of particles on the 
kinetic coefficients was carried out earlier by many investiga- 
tors. Leontovich (1)! solved in general form the Boltzmann 
equation as applied to the propagation of small disturbances 


1 Numbers in parentheses indicate References at end of paper. 
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in a gas at rest for the case w< 1, where w is the frequency 
of propagation of the disturbances in the gas and 7 is the re- 
laxation time for energy exchange between the translational 
and internal degrees of freedom of the particles. The solution 
was carried out by the Hilbert method for a continuous spec- 
trum of the values of the internal energy. Wang Chang and 
Uhlenbeck (2), using the Enskog formalism, obtained the 
transport coefficients in general form, in terms of the effective 
scattering cross sections, for a discrete spectrum of values of 
the internal energy. However, they did not calculate the 
transport coefficients directly. 

Pidduck (3) calculated the kinetic coefficients, using the 
model of “rough” spheres. Ishida (4) solved the Boltzmann 
equation for “solid ovaloids.”’ Jeans (5) used “loaded spheres” 
as his molecule model. However, these papers have one 
common shortcoming. The diameters for elastic scattering, 
which can be strictly calculated on the basis of quantum- 
mechanical collision theory, were replaced by various artificial 
mechanisms for the transfer of energy from the internal de- 
grees of freedom of the particles to the translational degrees 
of freedom. 

In this paper we use the probabilities of various types of 
molecular collisions, obtained quantum-mechaniecally by 
various investigators. The system of generalized kinetic 
equations is solved by the Enskog method (6). 

The system of integro-differential equations for the dis- 
tribution function f; of the components of the gas, with al- 
lowance for the internal degrees of freedom of the particles, 
has the following form 


fv f’; — Sif) > [1] 


where 
fi = f(r, €i, t) 
f'(r, c's, €'%, t) 
Ci, & = velocity and internal energy of ith particle before 
collision 
o's, €'; = same values after collision 
r = radius vector 
t = time 
= |gi,| 
b = collision parameter 
= azimuth angle 

K(é, 

e’, €';) = probability of collision of two molecules with a 


change of the values of €; or €; to €’; or €’; 


For given values of gi;, €:, and €;, K satisfies the normaliza- 
tion condition 


K(€i, €; > €’;) = 1 


The equilibrium distribution functions f; have the follow- 
ing form 


miCi? 
m 
= nj e kT (x ir) [2] 


where 


n; = number of particles of 7th kind per unit volume 

m; = mass of ith particle 

C; velocity of molecule of ith kind relative to local mean 
mass velocity of gas, Co. 


The local mean mass velocity is determined as follows 


ApriL 1960 


pi 
where 
p = gas density, p= nimi 
¢€; = mean velocity of the molecules of the ith kind in the 


laboratory system of coordinates 


In agreement with the Enskog method, we assume that the 
distribution functions f; do not depend explicitly on the 
time, but they are functionals of the macroscopic parameters 
of the medium, which vary with time. In the presence of small 
gradients of the macroscopic parameters, the distribution 
functions f; at each point deviate slightly from the local 
equilibrium functions f°. In the first approximation 


fi = + [3] 


where ®; is a function of C;, e; and of the macroscopic parame- 
ters of the medium. In Equation [3], f; coincides with Equa- 
tion [2] if the following normalization conditions are ob- 
served for the mean number of particles of each kind per unit 
volume, for the mean mass velocity, and for the temperature 


f [4] 
pc) = f [5] 
8 Eee 
dc; [6] 
J, 


= 


fime, = [8] 
= 0 (9] 
» 


In Equations [4 to 9] we have 


Here 


E; = —+86 


Inserting Equation [3] into [1] we obtain 


kn 
dig njli( Pi + [10] 
where 
ni nl + 


In the derivation of the left part of [10] we use the Maxwell 
transport equations, and also the condition that ®; is small. 
In Equation [10] 


Ey kT + eve (= 
i 


dO; = gijbdbdgdc; 
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Here U is the unit tensor, (0/Or)cp is the tensor whose com- 
ponents are the first row and first column, for example, are re- 
spectively 


d 


Uo, Uy and wp are projections of ¢ on the z, y and z axis respec- 
tively; U: (0/dr)co is the biscalar product of the tensors U 
and (0/dr)¢) and 


fe) 
or n n p) or 


Also p is the thermodynamic pressure; p equals nkT. 

The left half of the integral Equation [10] depends linearly 
on the gradients of the macroscopic parameters of the medium, 
and therefore ®; should have the form 

oln T n? > mim; 


— D;;-d; [ll 


A; = €;) 
= + €:) 
= (C3) 


In Equations [11 and 12], B; is a tensor and C;°C; isa tensor 
with a vanishing sum of diagonal components. For example, 
the components of the first row of the tensor C,°C; are U; — 
(1/3)C2, U:V;:, and U;W;. Inserting Equation [11] into 
[10] and equating the right and left side terms at d In T'/dr, 
(d/dr)co, and d; we obtain integral equations from which we 
can determine A;, B;, and Dj;. 

A total solution of the integral equations for A:, B;, and D;; 
can be obtained if detailed information is available on the 
probabilities of the various types of collisions of molecules; 
these probabilities are calculated from quantum-mechanical 
laws. Since the probabilities of the various types of collisions 
for the rotational degrees of freedom have still not been satis- 
factorily worked out, we shall consider here a gas with mole- 
cules that have only vibrational degrees of freedom in addi- 
tion to the translational ones. It is known that as the tem- 
perature is increased, the first to be excited are the rotational 
levels of the molecules, and the vibrational levels follow. 
However, as will be explained below, certain results obtained 
for the vibrational degree of freedom can also be generalized 
to include the case of any internal degree of freedom and, in 
particular, for rotations. In addition, the results will be of 
practical interest if the magnitude of the kinetic coefficient is 
determined essentially by the translational-vibrational 
energy transitions in the collisions. Of the three kinetic co- 
efficients mentioned in the foregoing, the bulk viscosity is 
indeed determined, under certain conditions, by these transi- 
tions. 

We shall consider below the effect of the following type of 
molecule collisions on the kinetic coefficients: 

1 Elastic collisions, occurring without a change in the 
vibrational quantum numbers v of the molecules. Elastic 
collisions of molecules will be described in terms of classical 
mechanics, using a spherically symmetrical interaction po- 
tential. 

2 Elastic collisions, accompanied with a change in the 
numbers v by unity, i.e., collisions that lead to an exchange of 
vibrational energy quanta between molecules. This type of 
collision can occur only between molecules of the same kind. 
The probabilities of such collisions were calculated by 
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Schwartz and Slawsky (7) who used the method of distorted 
waves, and are expressed in the following manner: 


where 
, 
v—>ov =-a 
27?mhw 
= 1 

and 

gn = normal component of the colliding molecules 

a = constant of the interaction potential between the 


molecules in the form Vye~*’ (a can be chosen to fit 
the constants of the Lennard-Jones potential) 


3 Inelastic collisions, accompanied by a change in the 
quantum numbers v of the molecules by unity. The probabili- 
ties of such types of collisions, according to Schwartz and 
Slawsky, are expressed in the following manner 
K(v, vy > v’, = X 


where 
ath ath 
Coefficient of Heat Conduction 
The heat flux vector q is equal to 
J, 


and is expressed in final analysis in terms of A; and D,;._ The 
integral equation for A; has the form 


(2: 1) Cc; = ni nili(Ai + Aj) [16] 


It can be shown that it is possible to solve Equation [16] by 
the Ritz variational method (8). 

As is known, the Ritz method yields rapidly converging 
sequences of functions, and therefore when A, is expanded in 
a series of known functions it is possible to retain merely the 
first terms of the expansion. 

In accordance with the form of the left part of Equation 
[16], we rewrite A; in the following form 

mC;? 
A; ORT + aj kT + a; Je 
where the coefficients a; are determined, if Equation [16] is 
solved by the Ritz method, in the following manner. Insert- 
ing A; into [16], taking the scalar product of the left and right 
halves of [16] successively with (m:C;?/2kT)Ci, (hwi/kT)C; 
and C;, and summing and integrating the resultant equations 
over all possible v; and C;, we obtain three linear equations for 
the coefficients ai“. We can obtain, analogously, (31-3) 
more equations. One of these is linearly dependent on the 
others. The lacking equation for the coefficients a;“ can be 
obtained from the normalization conditions [see Eq. 7 to 9]. 
Two of these [7 and 9], are automatically satisfied. The 
normalization Equation [8] establishes between the ai“ an 
additional dependence, which together with (3-1) equations 
obtained in the foregoing gives a complete system of equa- 
tions for 
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where 


To calculate 


mC? 


f ¢,- (Ac + 


+ A;)de; 


for the elastic collisions, we used the Lennard-Jones potential 
as the interaction potential. The elastic-collision integrals 
have also been tabulated for the Lennard-Jones function and 
are given in (9) in the form of effective temperature de- 
pendent scattering cross sections. These cross sections are 
denoted below by o; and = (o; + For collisions 
of the second type, which give the correction terms in the ex- 
pression for the heat conduction coefficient, we calculate the 
integrals of Equation [17] by using the model of a solid 
sphere with diameter o; and o;;._ The terms that take into 
i.ecount the collisions of the third type in Equation [17] can 
be neglected, if the difference hw; — hw; is not too small. 

The system of 37 equations can be simplified by replacing 
the coefficients a; by 6,“ in the following manner 


5 
a, 6, a, — 
2 k 


As a result of the substitution of a;‘® the system of 37 equa- 
tions for a; and 6; assumes the following form 


8n:Ria + Qua + + 
| 
L;;6;} = 15/4 


1 fah\? k7 3 


{Ma + 2M 56; 2M 


where 
R; = 


Ris = 
= MA76M.M; + 55M2 + 48M;*) — 5MX6M; + 


5M;) 

P;; = —M,(87M;M; + + 30M;*?) + 5Mi(5M; + 
6M;) 

Si; = 2M.(6M; + 5M;) 

Li; = —2M(5M; + 6M;) 

M; = m./mi; 

M; = m;/m;; 

= + m; 


In experiments one usually investigates the heat conduc- 
tion of gases in the stationary state, i.e., for a zero average 
velocity of motion of the molecules of each kind. In order to 
calculate the value of the heat flux corresponding to the ex- 
perimental conditions, we rewrite q in the following form 


I= njC;[(5/2)kT = + f x 
J, 


iT ) Cae, 


Putting C; = 


0, we get 
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where AX = A, + A is the coefficient of heat conduction. 

An corresponds to the coefficient of heat conduction of a 
mixture of monatomic gases and is calculated in the ordinary 
manner (8). 

X» is determined in terms of a; 


[18] 


Here C*,, is the heat capacity of the internal degrees of free- 
dom per molecule. The physical meaning of the terms that 
enter into the expression for \ can be readily interpreted 
using as an example a one-component gas, for which 


3 1 (2 2 ET 
=- —- k+ nDC* 
2m2 [3 \hwo} m [19] 
With 
2m2 


1 h\? kT 
Ny = mDC*, [r+ 2 + — (3) ] 
he m 


Here uy is the shear viscosity of the monatomic gas and D is 
the coefficient of self-diffusion (6). X, is determined only by 
elastic collisions of the first kind. Collisions of the second kind 
do not influence \,. The effect of inelastic collisions on the 
value of A, can be neglected because of the small probability 
of such collisions. 

The value of \», without the correction term, is determined 
only by the elastic collisions of the first kind and corresponds 
to the transfer of heat by self-diffusion of the excited mole- 
cules. The correction term in the expression for \» is due to 
elastic collisions of the second type. This correction is pro- 
duced by the “relay” mechanism of heat transfer, through 
exchange of quanta of internal energy between the molecules 
during the collisions. 

An equation similar to [19], without the correction for Xa, 
was first given by Eucken in the following form 


C*, 
[20] 


Equation [20] differs from [19] without the correction for the 
mechanism of heat transfer by a constant factor in the second 
term 


Actually, Equation [20] was obtained under the assumption 
that the mean free paths of the molecule are the same for the 
transfer of momentum and internal energy. In the derivation 
of Equation [19] it is shown that when the collisions of the 
first and third types are neglected, the internal degrees of 
freedom contributes to \ through the diffusion mechanism 
of heat transfer. The correction term in d, in [19] retains its 
additive meaning when rotations of molecules are taken into 
account, provided the probabilities of inelastic collisions for 
rotations are small. Equation [18], without correction for 
the relay mechanism of heat transfer, represents a generaliza- 
tion of a formula of the Eucken type to include the case of a 
multicomponent mixture of gases, and can be used to ealeu- 
late the coefficient of heat conduction with allowance for any 
arbitrary internal degrees of freedom, provided the probabili- 
ties of the collisions with change of internal state of particle 
are small. 
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Viscosity Coefficients 
For a single component gas, the pressure tensor P equals 
pU+>> f [21] 
The correction term in Equation [21] is expressed in final 


analysis in terms of the tensor B. The integral equation for 
B has the form 


m 2 k \ mC? 


- = n'I(B) [22] 


In accordance with the type of the left half of Equation [22], 
B can be written in the following manner 


m mC? hwv 
= b 0, — 7 


Using the normalization conditions [7 and 9], we can show 
that 


m how E, C*, {mC? 3 
th [ (i =) C: (= 


To determine 6;, we insert B in Equation [22], take a 
bisealar product of the left and right halves of [22] with the 
tensor (m/kT)C°C, and sum and integrate the resultant 
equality over allv and c. Since C°C:U = 0, the terms with 
bs drop out from the equation. The terms that take into ac- 
count collisions of the second kind cancel each other out. 
Neglecting small corrections, which take into account colli- 
sions of the third type, we obtain 


: 16 no® 


To determine bs; we multiply the tensor Equation [22] bi- 
scalarly with the unit tensor 


kT kT Ci \2kT 2 

and sum and integrate the expression term by term over all v 
andc. ‘Terms with }, drop out from the equation. The terms 
that take into account the collisions of the first and second 


type cancel each other out. As a result b3 will be determined 
only by the inelastic collisions 


C, \mkT) no®\ ho 
[23) 


where 


3 
“hd 


a* is a function of the temperature. 
The expression for the pressure tensor, Equation [21], can 
be reduced to the form 


P =nkTU — Qu 0% — SdiveoU 


Here yu is the shear viscosity and ¢ the bulk viscosity 
nkTh [24] 
= —nkT(C*,/C;)bs [25] 


(0/or)c% is a symmetrical tensor, the sum of the diagonal 
components of which is zero. 
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The expression for u coincides with the same expression as 
for a monatomic gas. Inasmuch as the probabilities of the 
collisions of third type are small up to dissociation tempera- 
tures, the influence of the vibrational degrees of freedom on 
the shear viscosity can be neglected. The rotational degrees 
of freedom, apparently, to modify Equation [24] to some 
degree at high temperatures. 

The bulk viscosity ¢ is determined only the elastic colli- 
sions. To explain the influence of rotations of molecules on 
the value of ¢ at ordinary temperatures, it is possible to use 
the phenomenological theory of bulk viscosity, developed by 
Mandel’shtam and Leontovich (10). Since at ordinary tem- 
peratures the probabilities of inelastic collisions with simul- 
taneous variation of both the vibrational and rotational state 
of the molecules are considerably less than the probabilities of 
inelastic collisions, in which either the vibrational or the ro- 
tational state of the molecule changes, we employ that part 
of the theory, in which the concepts developed by these 
authors are extended to include the case when the gas con- 
tains several relaxation processes, which are not interrelated 
to each other. It can be shown that when 


wTn K 1 
where 
Tn = relaxation time of any relaxation process 
w = angular frequency of ultrasound 
we have 
nk?T 
f= [26] 


In Equation [26], C*, is the heat capacity of the internal de- 
grees of freedom of the nth kind per molecule; C, is equal to 
Ci + 

It follows from Equation [26] that the bulk viscosity is 
made up additively of terms that characterize different relaxa- 
tion processes. If the heat capacity of the vibrational de- 
grees of freedom is not too small, and the relaxation time for 
the oscillations exceeds considerably the relaxation time for 
rotations, then ¢ is determined in practice only by the in- 
elastic collisions with change in the vibrational state of the 
molecules. At not too high temperatures, this holds for most 
gases. 

The applicability of Equation [24] for different tempera- 
tures in addition to the foregoing remark, is restricted by the 
method with which the Boltzmann equation is solved and by 
the methods used to calculate the integrals through which 6; 
is expressed. 

The limits of applicability of Equation [25], from the 
point of view of the method used to solve the Boltzmann 
equation, are determined by the inequality 


nkT > div 27] 


For the case of propagation of ultrasonic waves in a gas, con- 
dition [27] becomes 


p 


The limits of applicability of Equation [25] from the point 
of view of the methods used to calculate the integrals are 
determined as follows. 

These integrals can be reduced to the form 


J= (m/4kT) g? ¢ + (gn)g? sin cos 


where 


K(v, 1 > v + 1, 1) 
+1) 


gt= 
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After substituting x for cos y we obtain 


1 
J= f, (m/ART)9* + (¢x)adaxdg 


When going to the variables y = gx and x = t, domain I of the 
(g, z) plane, bounded by the lines g = gi, g = go, x = 0 and 
x = 1, is transformed into a domain in the (y, t) plane bounded 
by the rays y = git and y = get and the lines ¢ = O and ¢t = 1; 
we have 


Jt = Ju + Jin 


where J1 is the integral over domain I in the (y, x) plane; Ji 
is the integral in the region of the (y, t) plane bounded by the 
lines y = gi, = 1 and the ray y = got; Jt is the integral 
over the domain in the (y, ¢) plane bounded by the rays y = 
get, y = git, and the straight line y = g:. Making the sub- 
stitution z = 1/t? and y = y in Jit, and going to the limits 
a1 Oand go—> ©, we obtain 


0 


kT 
— f, (m/AkT)u* (y)dy 


This integral was calculated by Schwartz (7) by the steepest 
descent method under the following assumption 


e(22?m)/(a*h) > 1 


4 4 


where 


a*hm a*hm [29} 


m 


1/6 
(1 + + Dy 


1 
2 2 


(mgn*?/4) re 


e and 79 are the constants of the Lennard-Jones potential. 
Conditions [28] determine the limits of applicability of Equa- 
tion [25] from the point of view of the methods used to cal- 
culate the integrals. 

Equation [25] gives the correct order of magnitude of ¢. 
Using the constants of the Lennard-Jones potential for COs, 
namely ¢/k = 190 and rp = 3.996 A, we obtain a* 4.98 108 
em for T = 298 K by successive approximations from Equa- 
tions [29 and 30]. Assuming that o equals r;, where 7; is the 
classical distance of the closest approach of the molecules, we 
obtain = 3.54 A; ¢ = 1.15 g/cem-see; v/p< 0.1 cem/atm; 
aA = 17.5v/p, where v/p is in ms/atm; \ and as are the 
wave length and the coefficient of absorption of ultrasound 
respectively. According to data by Fricke (11), we have 
aA = 0.1 when p/y = 0.01. 

In conclusion, the author expresses his gratitude to Profes- 
sor E. V. Stupochenko for useful remarks concerning the 
present work. 
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Reviewer’s Comment 


The approach to the problem used in this paper is essentially 
the approach of Wang Chang and Uhlenbeck. The present 
authors, however, have carried the formal kinetic theory treat- 
ment somewhat further. The expression for the coefficient 
of thermal conductivity obtained by these authors is similar 
to that obtained by other methods by Hirschfelder (1). The 
authors use the phenomenological theory of bulk viscosity de- 
veloped by Mandel’shtam and Leontovich. It is interesting 
that the expressions that they obtain are very similar in form 
to those obtained by the kinetic theory approach by Wang 
Chang and Uhlenbeck. In the introductory comments the 
authors refer to earlier studies of the effects of internal struc- 
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ture on the transport coefficients. The authors, however 

apparently were unaware of the formal treatment and evalua- 

tion of the transport coefficients for the idealized model of a 
gas made up of rigid ovaloids (2-4). 

—C. F. Curtiss 
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Method of Calculation of the 


Transport Coefficients of Air 
at High Temperatures 


HIS paper contains theoretical calculations of the trans- 

port coefficients of air—the viscosity and heat conduc- 
tion—in the temperature range from 2000 to 8000 K and in the 
pressure range from 0.001 to 1000 atm. As in the ease of the 
thermodynamic quantities, there are still no well-developed 
experimental methods available for the direct study of the 
transport properties of gases at such high temperatures. Yet, 
in solving the most important problems of modern engineering, 
particularly the problem of jet propulsion, it is essential to 
know the transport coefficients of the air at high tempera- 
tures and over a broad range of pressures. Therefore the 
theoretical approach, as in the case of the thermodynamic 
properties, becomes particularly important here. However, 
there is a great difference between problems in statistical 
thermodynamics and problems in kinetic theory. In the 
former (in the case of ideal gases) the mechanism of inter- 
molecular (or interatomic) interactions does not play a sub- 
stantial role, and all that need be known is the energy levels 
of the particles (molecules, atoms, etc.) and their statistical 
weights. In kinetic theory on the other hand, the laws of 
interaction between particles, the probabilities of transitions 
upon collision, play a decisive role; they enter in a substantial 
manner into the gas kinetic equation. The particular difficul- 
ties involved in the calculation of the transport coefficients 
are essentially due to this circumstance. 

As is well-known, the statistical methods of calculating the 
transport coefficients, which Enskog and Chapman de- 
veloped by extending and improving the method of Hilbert, 
were derived only for central forces diminishing with sufficient 

apidity with increasing distance between molecules. Coulomb 

forces, and also inelastic collisions, particularly those ac- 
companied by chemical reactions and ionization, entail 
certain complications. In a considerable portion of the in- 
vestigated region of T and p in air, far-reaching dissociation 
and ionization processes take place, not to mention the excita- 
tion of the internal degrees of freedom. Therefore the cal- 
culation of transport coefficients is differently solved in differ- 
ent intervals of 7 and p. 

The probability of molecule collisions accompanied by a 
reaction is relatively small. Equally small is the probability 
of transition of translational energy into vibrational energy, 
and vice versa. The transition of translational energy into 
rotational energy and vice versa is found to be more probable, 
but even these transitions do not affect the momentum trans- 
fer noticeably. Therefore, the viscosity can be calculated 
from the formulas for gas mixtures, neglecting the inelastic 
collisions even in the presence of substantial dissociation. 
The important problem in this case is that of choosing the 
form of the potential of interaction and the numerical values 
of the constants that enter into the expression for this poten- 
tial. 

Section I of this paper is devoted to the principles of caleu- 
lating the viscosity of air under the conditions indicated in the 
foregoing. There are several known expressions for the in- 
Translated from Fizicheskaya Gazodinamika (“Physical Gas- 
dynamies’’), USSR Acad. Sci., pp. 39-58. 
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termolecular potential of the molecule interaction, and these 
expressions will be in good agreement with experiment if the 
constants are suitably chosen. However, the use of most of 
these expressions would not be correct for our calculations, for 
it would mean an excessive extrapolation of the experimental 
data to the region of very high temperature. For example, 
the term in the Lennard-Jones potential, proportional to r~'” 
(where r is the distance between molecules), which takes into 
account the repelling forces, is purely empirical in character. 
Yet, the role of repelling forces is large at high tempertures. 
These considerations have led us to choose the potential pro- 
posed by Buckingham, the form of which is theoretically justi- 
fied from a quantum-mechanical examination of the problem 
of interaction between two molecules. 

As regards the gas kinetic diameters of the atoms, the 
available data—experimental and _ theoretical—are quite 
scanty. The same atoms may interact with each other in 
accordance with entirely different laws, depending on the 
relative orientation of the spins of the outer electron shells 
of the colliding atoms. In particular, interactions of purely 
repulsive type or with a small potential well, may be realized. 
Since such interactions do not lead to the production of 
stable molecules, then naturally there are no spectral data 
that could be used to calculate the curve of the interaction 
potential. Yet most collisions do obey repulsive laws. 

To obtain the necessary values of the nitrogen and oxygen 
atomic diameters, which were not measured experimentally, 
we used interpolation of the experimental data on the diame- 
ters of the atoms of other elements, based on the regularity 
of the periodic system. 

An analysis of heat conduction is made complicated by the 
need of taking into account the excitation of the internal de- 
grees of freedom as well as the chemical transformations 
occurring during the collisions between molecules. In Section 
2 we analyze transport phenomena in a mixture of reacting 
gases. Recently many authors (Onsager, Haase, Eckert, 
Meixner, de Groot and others) have developed their own 
methods of investigating various macroscopic processes, such 
as diffusion, heat conduction, electric conduction, ete. These 
methods are known in the literature under the name of ther- 
modynamics of irreversible processes. The most interesting 
and newest results are obtained by the thermodynamics of 
irreversible processes in which several processes interact, for 
example, heat conduction and diffusion, heat conduction and 
electric conduction, etc. In gases and liquids, where diffusion 
proceeds freely, heat transfer is always connected with a 
diffusion flow of matter. However, in the absence of a reac- 
tion, the fraction of the heat flow due to diffusion is insig- 
nificant and exists only until the establishment of a steady 
state, in which diffusion and thermal diffusion balance each 
other. In the presence of a reaction, the amount of heat 
transferred with the diffusing substance increases considerably 
(by the amount of the “energy of chemical affinity’’), and, in 
addition, this component of the heat flow does not vanish 
with time, since there exists a nonvanishing diffusion flow of 
matter in the state of local chemical equilibrium. 
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The problem consists of determining the specific type of 
law of heat transfer under conditions of local chemical equi- 
librium, and in finding the phenomenological coefficients that 
enter into this law. The first part of this problem is solved 
in the present paper by methods of the thermodynamics of 
irreversible processes. We succeed in expressing the law of 
heat conduction in terms of those transport coefficients that 
can be determined by methods of modern kinetic theory 
(5)1, making it possible in turn to calculate these coefficients 
by means of one molecular model or another. In addition, 
no assumptions are made in this paper concerning the absence 
of a pressure gradient in the medium, as was done, for ex- 
ample, by Haase (4). This is necessary because heat exchange 
occurs in gasdynamies, as a rule, under nonisobarie condi- 
tions. 


I Viscosity of Air in the Presence of Molecular 
Dissociation 


1 Preliminary Remarks 


In this chapter we report the results of an examination of 
the viscosity of air in the temperature range from 2000 to 8000 
K, for pressures from 0.001 to 1000 atm. ‘The viscosity of a 
gas in which no chemical reactions take place depends on the 
temperature and is almost entirely independent of the pres- 
sure. At high temperatures, on the order of several thousand 
deg K, the molecules of the air begin to dissociate into atoms, 
and consequently chemical reactions take place in the air and 
its viscosity becomes dependent on the pressure. 

It is possible to obtain experimental values of the viscosity 
of the air only up to 2000 K. We are interested, on the other 
hand, in the viscosity of air in the temperature range from 
2000 to 8000 K. We must therefore resort to theoretical 
methods of determining the viscosity. Among the papers 
devoted to this problem are those of Chapman and Enskog, 
who have derived approximate expressions for the viscosity 
and heat conduction of gas in terms of the so-called ‘collision 
integrals,” which depend on the law of interaction between 
the gas particles. 

The laws of interaction between atoms and between 
molecules are known only approximately, and this in turn 
makes for further inaccuracy in the expressions for the 
transport coefficients. 

There are many equations for the dependence of the vis- 
cosity on the temperature, in accordance with the various 
representations concerning the character of the intermolecular 
interaction. Thus, the power law dependence of 7 on the 
temperature 


n = constant [1] 


corresponds to a power dependence of the interaction poten- 
tial on the distance. However, when 7 < 1000 K, a better 
agreement with experiment is given by the Sutherland 
formula 

aT’? a = constant 


1+ C/T C = constant [2] 


which takes into account the weak forces of attraction be- 
tween molecules. But at high temperatures (~2000 K) it is 
possible to neglect the term C/T’, since C is on the order of 
~100 K. Then Equation [2] becomes 7 = aT’/?, which is 
in poorer agreement with experiment than the power law 
formula. Other empirical equations are also known (5). 

The foregoing formulas were obtained from very simple 
concepts regarding the intermolecular forces, not correspond- 
ing to the present-day status of this problem. Therefore, to 
obtain a more accurate expression for this viscosity, we turn 
to later research data on intermolecular forces. 


1 Numbers in parentheses indicate References at end of paper. 
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It must be noted that at high temperatures the viscosity of 
air can be determined only with small accuracy, since at such 
temperatures free atoms appear in the air, and the laws of in- 
teraction between these atoms are still unknown. It can be 
considered established, however, that the gas kinetic diameters 
of the molecules of nitrogen and oxygen are close to those of 
the atoms of nitrogen and oxygen, and therefore in our ap- 
proximate calculations of the transport coefficients we can 
assume that at such temperatures air represents a binary 
mixture of gas components: Molecular (No, Ox, NO) and 
atomic (N, O, A). 


2 Interaction Between Molecules and Viscosity of the 
Molecular Component of the Mixture 


A quantum mechanical analysis makes it possible to 
establish that at large distances (compared with the radius of 
the molecule) nonpolar molecules attract each other with a 
force, the potential of which can be represented in form of a 
series 


ro 

Obviously, at large distances the first term is the only one of 
importance. Without any theoretical justification, the re- 
pulsion forces were also expressed by means of a power law, 
with a suitable choice of the power exponent. A well-known 
potential for the intermolecular interaction is the Lennard- 
Jones potential 


U = — [4] 


where 


€ = depth of potential minimum 
= 
rm = distance between centers of molecules 


For the molecules O2 and No», the value of €/k (where k is 
the Boltzmann constant) is approximately on the order of 100 
K, rm is approximately 4 A, and € and rp, can be determined 
from the experimental data, both from the transport co- 
efficients and—independently—from the values of the virial 
coefficients. It is important to note here that the values of 
rm and € determined by such independent methods were found 
to be practically equal (Table 1). 


Table 17 
Virial Virial 
Gas Viscosity coefficients Viscosity coefficients 
air 97.0 99.2 3.617 3.692 
Oo 91.46 95.9 3.681 3.72 
No 113.2 117.5 3.433 3.58 
@ See (6). 


“Collision integrals” for the Lennard-Jones potential were 
tabulated by Hirschfelder, Bird and Spotz (6). The values of 
the viscosity, calculated from the data of these calculations, 
are in good agreement with the experiment. It turns 
out, however, that at high temperatures the dependence 
of the viscosity so calculated on the temperature obeys a 
power law 

n~ [5] 


where n = 0.666. 

This is explained by the fact that at high temperatures the 
principal role is assumed by the repulsion forces, correspond- 
ing in the Lennard-Jones formula to the term ~r~!*, which 
does not reflect correctly the nature of these forces. Conse- 
quently, the use of the Lennard-Jones potential in the calcu- 
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lation of the viscosity is not reliable at high temperatures. It 
can be considered as proved (7) that the repulsion forces de- 
pend on the distance exponentially. One of the most 
useful expressions for the law of intermolecular interaction, 
with account of this character of repulsion forces, is the modi- 
fied Buckingham potential 


€ 6 
1 — 6/a | (6) 


Here q@ is a dimensionless parameter. Depending on the type 
of molecule, a may range from 12 to 17 or higher. The 
values of a, fm and € can be obtained by comparing the results 
of theoretical calculation with experimental data on viscosity, 
the second virial coefficient and the interaction of molecules 
in crystals at temperatures near 300 K. 


Table 2° 
Gas a Pin; A €, erg 
Ne 14.5 3.147 38 
A 14 3.866 123.2 
Kr 12.3 4.056 158.3 
Xe 13.0 4.450 231.2 
CH, 14.0 4.206 152.1 
Ne 4.040 113.2 
117.0 4.011 101.2 
rare) 4.10 132 
(17 3.937 119.1 
@ See (8). 


The discrepancy between the values of the parameters ob- 
tained from different experiments for the same diatomic 
molecules is due to the fact that, strictly speaking, the poten- 
tial of Equation [6] is applicable only to spherically symmet- 
rical molecules. Diatomic molecules can be considered spheri- 
cally symmetrical only approximately. The nonsphericity 
influences least of all the equilibrium properties, i.e., the 
second virial coefficient. 

The error introduced by nonsphericity into the viscosity of 
N: is found to be on the order of 2 per cent at temperatures 
above 600 K; i.e., it is fully within the permissible accuracy. 

As yet we have no corresponding data for O2 molecules but, 
using the similarity between the Lennard-Jones potential and 
the potential of Equation [6] at not too small values of r(r > 
0.3rm), we can determine these parameters. The viscosity 
(n), calculated by the Enskog-Chapman method and Equa- 
tion [6], has the following form 


5 (mkT\'/2 fn 
where 
m = mass of molecule 
k = Boltzmann’s constant 
T = absolute temperature 
T* = (kT/e) 
€ = same meaning as in Equation [4] 


The functions 22: and f, are not expressed analytically. 
They were tabulated by Mason (7), and by Mason and Rice 
(5). It can be seen from Equation [7] that the ratio 


= ¥(T*) 


where 79* is a certain fixed value of 7*, is the same function 
of the reduced temperature 7’* for all gases. If we calculate 
a similar function ¥(7'*) for the model of elastic solid spheres, 
then at high temperatures (7 > 4000 K) its values will be 
1.5 or 2 times as small as the corresponding values of ¥(7"*), 
calculated with the use of the Buckingham potential. 
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3 Gas Kinetic Diameters of Atoms and Viscosity of the 
Atomic Component of Air 


There are only few known expressions for the interaction 
forces between atoms, from which to calculate the transport 
coefficients. Only those inert gases have been studied for 
which the interaction forces between atoms and their gas 
kinetic diameters? are known with sufficient accuracy. See 
Table 3. 


Table 3 

Gas or, A 
He 273 2.19 
Ne 293 2.59 
A 273 3.66 
Kr 273 4.18 
Xe 273 4.93 
Kn 288 5.2 


The outer shell of the atom of an inert gas is closed, and 
this facilitates considerably the theoretical and experimenta! 
investigations. In particular, the attraction forces betwee: 
the atoms are inversely proportional, in first approximation, 
to the sixth power of the distance. 

The total interaction between the atoms of inert gases is in 
good agreement with the Lennard-Jones potential, making it 
possible to determine the parameters ¢ and r,, for these gases 
with satisfactory accuracy. 

The situation is quite different for atoms with an unfilled 
electron shell. The energy of the interaction between such 
atoms either has a very strong minimum or corresponds to 
pure repulsive forces, depending on the relative orientation of 
the spins of the external electrons. For an interaction with 
a minimum, the dependence of the interaction energy on the 
distance is sufficiently well known from spectroscopic data. 
Several empirical formulas have been obtained, such as the 
Morse function, the Hulbert-Hirschfelder function, ete. 

However, spectroscopic data make it possible to determine 
the behavior of the curve only near the minimum, and only by 
extrapolation can one obtain information on the behavior of 
the curve at large distances (compared with r,,) or in the re- 
pulsion part. As to curves of purely repulsive type, there are 
as yet neither experimental nor reliable theoretical data, since 
in spectroscopy one deals with interactions that correspond 
to a curve with a minimum. The interaction energy was 
calculated theoretically for hydrogen; for more complicated 
atoms, the calculations entail considerable difficulties con- 
nected with solving many electron problems. 

Consequently, in calculating the coefficients for the atoms 
it becomes necessary to forego any dependence of the interac- 
tion energy on the distance and to use the model of elastic 
rigid spheres. However, even in this formulation of the 
problem we find that no reliable experimental data are as yet 
available for the solution of this problem. 

The effective cross sections for the viscosity can be deter- 
mined if one knows the effective cross sections for scattering, 
since a definite connection exists between the cross section for 
the viscosity and the cross section for scattering (see, for ex- 
ample, Mott and Massey, “Theory of Atomic Collisions’’). 
Specifically, if we denote by /i#dw the probability of scattering 
of a particle within a solid angle dw at an angle # to the initial 
direction, then the total effective scattering cross section Q is 
determined as follows 


Q = "1(8) sin (8) 
~ 2'The gas kinetic diameter or at a temperature 7’ is given by 


or = 
where 7(7’) is the viscosity at the temperature 7’. 
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and the effective cross section for the viscosity is 
Q = sin’ Sud [9] 


The angle # depends on the law of interaction between the 
atoms. Since there is still no exact theoretical expression for 
this law, it becomes necessary to resort to experimentally de- 
termined effective scattering cross sections. 

However, the presently available experimental data are 
scanty, and furthermore it is found that the values of the 
diameters of the same atoms differ greatly from each other 
when scattered by different gases. 

We need more reliable information, in the form of data on 
viscosity and diffusion. Data on viscosity, for example, yield 
the diameters of atomic hydrogen [9] [cited in (6) ] 


or = 2.52A at = 273K 
or = 248A at = 373K 


The diameter of deuterium was found to equal the diameter 
of ordinary hydrogen. If the temperature dependence of or 
is presented in the form 


or? = 6X1 + C/T) [10] 


we obtain for different atoms the values of o, and C shown in 
Table 4. 


Table 4 
He 1.82 173 
Ne 2.52 128 
A 2.99 142 
Kr 3:22 188 
Xe 3.55 252 
H 2.39 31 
D 2.39 31 


From experiments on diffusion we know the gas kinetic 
diameter of sodium atoms 


or =3.7A at 7 = 291K 
In caleulating this value it was also assumed that 
Oy = (1/2)(0, + oe) 


In addition, there is available information on the diameters ¢ 
for atoms of zine, cadmium and hydrogen. Thus, the in- 
formation on gas kinetic diameters is extremely scanty. 
Consequently, one can speak only of an approximate estimate 
of this quantity for the oxygen and nitrogen atoms, which are 
of interest to us. 

The meager data cited show that the gas kinetic diameters 
of the atoms are definitely related to the number of the atoms 
in the periodic system of the element. The most important 
fact is that as the atomic number increases, we have not a 
monotonic increase in the atom diameter, but a periodic varia- 
tion. Furthermore, the diameter does not increase with in- 
creasing atomic number within each period, but, on the con- 
trary, decreases. Thus, the diameter of the hydrogen atom at 
T = 273 K is 2.52 A, which is noticeably greater than the 
diameter of helium (2.19 A) and is approximately equal to the 
diameter of neon (2.59 A at 293 K), in spite of the fact that the 
latter has 10 electrons as opposed to the single electron of hy- 
drogen. Furthermore, 7 of hydrogen is 2.39 A, whereas that 
of helium is 1.82 A and neon 2.25 A. 

The diameter of the sodium atom, the first in the third 
period of the periodic system, is 3.7 A at 291 K, considerably 
greater than the diameter of neon which has the next lower 
atomic number, and is somewhat greater than the diameter of 
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argon (3.66 A), which is at the end of the third period. Finally, 
the diameter of lithium, located vertically between hydrogen 
and sodium, should also be intermediate between the two, 
i.e., approximately 3 A, considerably in excess of the gas 
kinetic diameter of neon (2.52 A), with which the second 
period ends. 

The foregoing data indicate that, at least in the first periods, 
the gas kinetic diameter of the atoms decreases somewhat 
from the start to the end of the period, and then, in going from 
the inert gas which ends the period, to the alkali metal which 
begins the new period, the diameter increases sharply. 

There is another physical quantity which also depends on 
the atomic number of a given atom in the periodic table. 
This is the ionization potential, with which the atomic radius 
is directly connected (quantum mechanically). The values 
of the ionization potential have been well investigated for al- 
most all gases. The ionization potential increases with the 
atomic number within the period and decrease sharply in 
going to the next period. It is particularly significant to us 
that the change in potential within the period is practically 
monotonic, which apparently indicates a monotonic variation 
in the gas kinetic diameter, too. 

In connection with this, one must note the following fact. 
The products of the gas kinetic radii o7 of the atoms of inert 
gases times the ionization potentials V; of the same atoms are 
nearly equal to each other (Table 5). We can therefore as- 
sume that for each group (more accurately, subgroup) of ele- 
ments in the periodic system, the products of the gas kinetic 
radii of the atoms times the values of the ionization potentials 
of the same atoms are nearly equal to each other, i.c., that this 
product is practically constant for a given group (subgroup) 
of the periodic table. An indirect confirmation of this group 


Table 5 
Gas Vi, volts or(A) 
He 24.46 2.19 53.6 
Ne 21.47 2.59 55.6 
A 15.68 3.66 
Kr 13.93 4.18 58.3 
Xe 12.08 4.93 59.5 
Rn 10.7 5.2 55.7 
* Hodgeman, Ch. D., “Handbook of Chemistry and 
Physies”’ 1951. 


is the fact that the magnitude of the ionization potential— 
the energy necessary to remove an electron from the atom— 
is inversely proportional to the distance between the given 
electron and the nucleus of the atom, and this distance, or the 
quantum-mechanical radius, should be directly proportional 
to the gas kinetic radius of the atom. One can assume further- 
more that this product changes over the period as a function 
of the number of the electrons on the outer orbit or as a func- 
tion of the total spin of the outer electron shell of the atom. 

However, since the changes in the collision diameter over the 
period are not very large, we shall not introduce a large error 
if we assume them to be linear within each of the periods. 
To find the diameter o,; we can take the arithmetic mean of 
the diameters of hydrogen and sodium, between which 
lithium is located 

2 

We next locate the diameter of neon at high temperatures. 
The value of o~ from Table 4 cannot be used, since Equation 
{10] cannot be applied at such high temperatures. Nor are 
there any experimental data for temperatures from 3000 to 
10,000 K. 

However, since the electron shell of the neon atom is filled, 
we can apply to it the Buckingham molecular potential. 
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Therefore the gas kinetic diameter of neon can be deter- 
mined at the required temperatures from a knowledge of the 
“eollision integrals” and of the transport coefficients for neon 
at the same temperatures. The gas kinetic diameter obtained 
in this manner is 


One = 1sA 


If we assume that the gas kinetic diameter of the atom 
varies linearly along the period we obtain for the oxygen atom 
go = 2.16 A, and for the nitrogen atom ox = 2.3 A. It is 
obvious that within the limits of the presently attainable 
accuracy, one can assume for the general gas kinetic diameter 
for the atomic component 


o, = 2.2 K (or 2.3 A) [11] 
Knowing the gas kinetic diameter, we can calculate the 
viscosity by the Chapman formula 


n = [12] 


where m, is the mass of the atom. 


4 Viscosity of Air 

To calculate the viscosity of air, which is a mixture of the 
atomic and molecular components, it is necessary to know 
the law of interaction between the atoms and molecules. 
Information on this interaction is also practically nonexistent, 
and consequently we must again use the model of elastic 
spheres. The collisions between the molecule and the atom 
are considered here as collisions between elastic spheres with 
effective diameter 


O12 = (1/2)(o1 + o2) [13] 
where 


o; = diameter of atoms as determined in the foregoing 


o2 = diameter of the molecules, calculated from the re- 
lation 
= o( [14] 
n(T) = viscosity of the molecuarl component at temperature 


Obviously, the diameter of the molecule depends on the 
temperature. 

Expressions for the viscosity of a binary mixture of gases 
were derived by Chapman and Enskog, but since the viscosity 
of the atomic component has been calculated thus far with 
low accuracy, it is logical to employ the earlier formula 
derived by Enskog on the basis of an examination of the 
transport processes in a gas consisting of molecules, with an 
interaction that diminishes rapidly with distance (10) 


= 
Ne 
1 + 1 + 
2d(qine)'/? 
1 — + [1 + [15] 
where 
™ = viscosity of the atomic component 
Ne = viscosity of the molecular component 


relative molar fractions of the first and second com- 
ponents 


1 2m 5m, + 
3 1 + Me m, + 


1 /2 3m, + 
3 2 m, + Me m + Mm 
8\'/2 (mm:)*/* 
a = 
9) (m + 
o1 = effective diameter of atom 
C2 = effective diameter of molecule 
C12 = effective diameter for collisions between atoms ani 
molecules 
m,M2 = masses of atom and molecules 


We note that if we calculate the viscosity of air with the 
aid of the exact Chapman-Enskog formula, which has been 
derived for elastic spheres, the result deviates from that 
calculated by the approximate formula [15] by approxi- 
mately 2 or 3 per cent which obviously includes the inaccurac) 
in the calculation of the viscosity of the atomic component. 

Data on the composition of air at high temperatures, whic! 
are essential for the calculation of x: and 22, can be found ini 
the report of the Laboratory for the Physics of Combustion, 
the Power Institute for 1955 (5). This report contains data 
on the composition and thermodynamic properties of air at 
temperatures from 2000 to 8000 K and pressures from 0.001 
to 1000 atm. Itisseen from these data that total dissociation 
of the molecules of air occurs only at pressures less than | 
atm. For higher pressures at 8000 K a certain amount of 
undissociated molecules still remain; at pressures of several 
hundred atmospheres, approximately 50 per cent of the 
molecules remain undissociated. 

It must be noted further that at temperatures near 8000 or 
10,000 K, and for very low pressures even somewhat earlier, 
the degree of ionization becomes so substantial, that further 
calculation of the transport coefficients without allowance for 
ionization becomes meaningless. 

The results of the calculations of the viscosity of air in the 
corresponding temperature intervals and pressure intervals 
are contained in the reports of the Laboratory for the Physics 
of Combustion, the Power Institute, for 1955 (5 and 11), de- 
voted to calculation of the transport coefficients (Fig. 1). 


Il Heat Conduction in a Binary Mixture of 
Chemically Reacting Gases 


1 Heat Flow in a Binary Mixture of Gases 


Unlike the viscous stress tensor, which always depends only 
on the velocity gradient tensor, heat flow in liquids and gases 
in which diffusion is possible is closely related with the diffu- 
sion flow of the material, and depends therefore not only on 
the temperature gradient, but also on the gradient of the 
chemical potential. In the presence of a chemical reaction, 
the role of those heat flow components, which arise as a result 
of diffusion of molecules, increases greatly. 

To investigate the problem of heat conduction under 
conditions of chemical reactions it is necessary first to agree on 
the definition of the concept “heat flow.” The heat flow is 
determined by the form of the energy equation employed. 
In this paper we shall write the energy equation in the form 


oF oE Ou; 
ol Ory OX; 
where 
E = internal energy per unit mass of mixture 
Tix = stress tensor 
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uz = hydrodynamic velocity 
p = density of the mixture 
1 = heat flow vector 


Repetition of an index denotes summation over this index. 
In a binary mixture of A; and Ag 


E = mE, + mE: [2] 
where 
pi _ Pi 
ni 
pit ps 
™, M2 = mass fractions of components A; and Az in mixture 
pi, p» = densities of components A; and Z, 


E,, EF; = internal energies per unit masses of A; and A» 
The diffusion flux I™ of the component A; is defined as 


= p(u uw) = — (u® — = I% [3] 
p 
where 
u = hydrodynamic velocity, determined from the re- 
lation 
u= nu + nou [4] 
u™, u® = velocities of the components A; and A» 
I? = diffusion flux of component A» 


We obtain from the thermodynamics of irreversible proc- 
esses (1 and 2) 


= — = TY ( ) — In (6) 
where 
Ia, Lu, = phenomenological coefficients 
ge), oi = partial specific thermodynamic potentials of 
components A; and A» 
T = absolute temperature 
= grad = — + — + — 
On; Ox2 


| 

a 6000 9000 12000 
Fig. 1 Coefficient of viscosity in millipoise (1 poise = gm/cm- 
sec) asa function of the temperature for various pressures p, atm 
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Equations [5 and 6] have been derived assuming no mass 
forces. 

Onsager’s principle of the symmetry of the kinetic coeffi- 
cients was taken into account in the derivation of Equations 
[5 and 6]. 

The coefficients Lg, Ly and La determine fully the heat 
flow and the diffusion flow and consequently all other co- 
efficients introduced for the description of diffusion, and heat 
conduction should be expressible only in terms of these three 
coefficients. 

Let us introduce a new phenomenological coefficient of 
“heat transfer” [see (1 and 2) for more details] 


La = 7] 
We then have instead of [5 and 6] 


— 


1 
1 = — { Leg Lu Q*} VT [9] 


or, denoting 
— Lg Q*)/T (10] 
we obtain the following expression for the heat flow 
1@ = Q*I — (11) 
The stationary state, in the absence of a chemical reaction, is 
determined by the condition 
ID) =0 [12] 
Thus we have in the stationary state 
I@ = —\AVT [13] 
and consequently X is the coefficient of heat conduction in 
the stationary state. 


Next, assuming that ¢ = ¢(T7, p, m), and employing the 
Gibbs-Duhem equation 


dg?) 
+m- 
om on 


m =0 [14] 


we obtain from Equation [8] the following expression for J“ 


1 (1) 
p,T 


m \ Om 
Vm — (Q* — +h) [15] 
where 
Pp = hydrostatic pressure 
hi, ho = specific enthalpies of A; and A», 
v1, = specific volumes of A; and A, 


Comparing Equations [11 and 15] we find that even in the 
absence of a temperature gradient the heat flux does not 
vanish if the concentration gradient or the pressure gradient 
does not vanish. Inserting [15] into [11], we find that the 
coefficient of VT is 


K = + LnQ*(Q* + he)(1/T) [16] 


x is the coefficient of heat conduction, which determines the 
heat flow in the state of uniform mixing, i.e., when Vp = 0 
and Vm = 0; it is obvious that « does not equal A. If there 
is no reaction, the initial levels from which the energies of the 
molecules of A; and Az are reckoned can be chosen independ- 
ent of each other. In particular, we can reckon the energy 
of each molecule from the lowest energy level of the cor- 
responding molecule. Then the difference h; — he is small. 
Then, as seen from theoretical calculations based on kinetic 
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theory and from experiment, the second term in Equation 
[16] is small compared with the first, and only one term, 
—d grad T, need be retained in Equation [11]. The calcula- 
tions for this case are known from the kinetic theory of gases 
(3), and lead to the expression 


I@ = pkr(C: — C2) — AVT [17] 
where 


kr thermal diffusion ratio 
C:, C: = mean thermal velocities of A; and As 


Inasmuch as C; — C2: equals u“ — u™, we can rewrite 
Equation [11] with the aid of Equation [3] in the form 


I@ = Q*(pip,/p)(C: — C2) — AVT [18] 


Comparing [17 and 18] we see that the coefficient of heat con- 
duction, calculated in the Enskog-Chapman theory for a 
binary mixture of gases, coincides with the coefficient of heat 
conduction for the stationary state, Equation [10]. 


2 Heat Flux ina Reacting Mixture of Gases at Equilibrium 


We consider the case when a chemical reaction takes place. 
From the molecular-kinetic point of view, the heat flux in the 
reacting mixture contains additional terms, over and above 
those contained in the absence of reaction. Indeed, the heat 
flux is expressed as 


+ 
where 


E, = (1/2)mCo? + UM 
(1/2)m.C2 + U® +e 


and are the energies of the molecules A; and Az respectively; 
U® and U™ are the energies of the internal degrees of free- 
dom of molecules A; and A; € is the difference in the lowest 
energy levels of molecules A; and Ae. 

We shall assume henceforth that equilibrium exists between 
the translational and internal degrees of freedom of the 
molecule. 

In the absence of a reaction, the initial levels from which 
the energy is reckoned can be chosen independently, and 
consequently we can assume € = 0. This cannot be done in 
the presence of a chemical reaction, which leads to the ap- 
pearance of an additional heat flux term eC, which can also 
be written in more symmetrical form as aC™ + eC. 

The physical meaning of these terms is obvious: In the 
presence of reaction, the heat flux includes also the transfer 
of chemical-affinity energy, caused by the thermal motion of 
the molecules. The mechanism of the transfer of this ene-gy, 
like the mechanism of transfer of the internal energy, is of the 
diffusion type. 

Assume that the reaction begins at ¢ = 0 and there was no 
reaction before that time. If 7 is the relaxation time of the 
reaction, then at t < 7 the reaction does not influence the 
heat flow or any other macroscopic process. Let us examine 
what takes place at t >> 7. The presence of a reaction does 
not change the form of the general thermodynamic equations 
for the heat flux, Equations [5 and 11]. This is due to the 
fact that the reaction is determined by scalar quantities (the 
speed of reaction and the chemical activity), whereas the 
heat flux is a vector. A vector can be linearly related with a 
scalar only if the coefficient of proportionality is a vector. 
Thus, the phenomenological coefficient that characterizes 
the connection between the heat flux and the chemical reac- 
tion should be a vector. But, on the other hand, this co- 
efficient should be isotropic, since it pertains to an isotropic 
medium. The only isotropic vector is a null vector, and this 
indeed explains the absence of a direct connection between 
the heat flux and chemical reaction. 

However, the presence of reaction leads to a change in the 
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stationary state. If the stationary state is defined by the 
condition 7 = 0, in the absence of reaction, then, in the 
presence of reaction, if t > 7, the stationary state is a state of 
local chemical equilibrium, determined by the condition 


— = — = 0 [19] 
where 


g™, @® = partial molar thermodynamic potentials of com- 
ponents A; and A, 


Vi, Ve = stoichiometric coefficients of reaction in a binary 
mixture of gases 
— = 0 [20) 
Finally 
veo = Vip v2 = _ [21] 


where y; and ye are the molecular weights of A; and A». 

Since Equation [19] holds at each point, then 

Vine™ — =0 (22) 

thus, in the case of a chemical reaction, the diffusion flux does 
not vanish in the stationary state. 

To calculate the heat flux and the conditions of local chemi- 
cal equilibrium, we make use of Equations [11 and 15}. 
Inserting J from Equation [15] into [11] we find 


IM = — )Vp — X 


Ne om p,T 7 


We can use Equation [19] to rewrite the conditions of local 
equilibrium [22] as 
hg 4 (6%) Vn =0 [24] 
7 Ne dm 
where we introduce the symbols 
vp = Vo Ve [25] 
hr = = vhs [26] 


After obtaining Vm from Equation [24] and substituting 
into [23], we find 


po — 
Bip 
TI + Vp 
Q* +h hr 1 
4 £1Q* A> V1 


where v = mt; + 72v2 is the specific volume. Since, because 
of the conservation of mass, v4; equals veye, or 


[27] 
we have 
hp = (hi — he) [28] 
vim + = [29] 
that is 
+ vr [30] 


Thus, if the stationary state, under the conditions of a 
chemical reaction in a binary mixture of gases, is the state of 
local chemical equilibrium, then the term proportional to 
grad p drops out of the expression for the heat flux (exactly as 
the stationary state without a reaction). Equation [30] is 
the simplest expression for the heat flux under conditions of 
chemical equilibrium, but the coefficients Zy and the heat 
transfer Q* are still not in general use. In connection with 
this, it is advantageous to obtain an expression for the heat 
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flux in terms of the better known coefficients of diffusion and 
thermal diffusion. We introduce the coefficient of ordinary 
diffusion Dy: by means of the equation 


= —pDeVm [31] 
at VT = Oand Vp = 0, and a coefficient of thermal diffusion 

I = [32] 


at Vp = 0, Vm = 0. 
Then the coefficients Ly, and Q* are expressed in terms of 
Dy and Dy: in the following manner 


= [33] 
7 
12 1 


From Equations [30, 33 and 34] we find 
(Di2")? 
= + p — | -— 
Om 
— he)? 
(dg /Om)p, rT 
This expression for the heat flux is equally useful for a react- 
ing binary mixtures of liquids or of gases. Let us consider 


next the case when the reacting mixture can be considered as 
a mixture of ideal gases and then 


+ 2pmmDi27(hi — ht VT = [35] 


RT P(m/ 
1 = go(T) + — In 36 
(1) RT 1 
(=) RT [37] 
om p.T ™ Me + 


From this we obtain, for a binary mixture of reacting gases, 
the following expression for the heat conduction 


(D127)? pkT 
=At+ + 
hy — he)? 
Domne + Nop) + (hi — he) [38] 
or, introducing 
kp = [39] 


the thermal diffusion relation, and 


= [39a] 
Mem + Mite 


the mean molecular weight, we obtain finally 


Mike 


Pp h, — ho\? — he 


This expression determines the coefficient of heat conduc- 
tion of a reacting binary mixture of gases in the state of local 
equilibrium. The thermal flux is then expressed as 


I@ = — [41] 


The coefficient Apo is expressed in terms of the following trans- 
port coefficients 


AX = coefficient of heat conduction of stationary state of the 
nonreacting gases 
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Dz = coefficient of ordinary diffusion, defined by Equation 
[31] 
ky = thermal diffusion ratio, see Equation [39] 


The definition of kr includes also the coefficient of thermal 
diffusion D,,7, which is defined by Equation [32]. 

If the number of inelastic collisions accompanied by the 
reaction is small compared with the total number of collisions, 
then all these coefficients can be calculated by the ordinary 
techniques of modern kinetic theory of gases (the Chapman- 
Enskog method). The coefficient » includes corrections for 
the heat conduction of the internal degrees of freedom (i.e., 
corrections of the same type as introduced by Euken). 

Inasmuch as [(Ai — he)/RT]u is a large quantity, and kr 
is a small one, the second term of Equation [40] is small com- 
pared with the remaining terms. Consequently, the third 
term is greater than the fourth, but not enough so as to permit 
the neglecting of the latter in calculating the heat conduction. 
The third and fourth terms are comparable with the value of 
A, and in some cases exceed it considerably. 

The problem of heat conduction in a binary mixture of 
reacting gases was the subject of a paper by Nernst at the 
beginning of this century [this work is mentioned in (4)]. 
Nerst considered a mixture of ideal gases and did not take 
thermal diffusion into account; he therefore obtained only 
one of the three terms that supplement the ordinary heat 
conduction (A), namely, a term proportional to the square of 
(hi — he)/RT. This problem was analyzed by modern 
thermodynamic methods by Haase (4). Haase arrived at an 
expression similar to that obtained in this paper, but without 
the first correction term (proportional tokr*). This difference 
is due apparently to the difference in the definition of the con- 
cept “heat flux.” In addition, Haase considered a case in 
which the reaction occurs in a mixture without a pressure 
gradient, i.e., a case having little interest in hydrodynamics. 

We note that the transport coefficients for a binary mix- 
ture, calculated in (3), are related to the coefficients in Equa- 
tion [40] by the relations 


Al =X [42] 
D2} Dy [43] 
D7! = DyT [44] 
kr! = [45] 
where 
M = heat conduction coefficient for a binary mixture of 
gases, after Chapman 
D+ = diffusion coefficient, after Chapman 
DT! = coefficient of thermal diffusion, after Chapman 
ky! = thermal diffusion ratio, after Chapman 
&, & = molar fractions of the components A; and A 


The connection between the molar (£) and mass (7) frac- 
tions is given by the equations 


(m/bi) + + Nev 
ni [47] 


Midi + bake & + (&/v) 


where v = p,/v2 is the ratio of the stoichiometric coefficients of 
the reaction. 

To solve any hydrodynamic problem it is necessary to 
write also, in addition to the energy equation, the continuity 
equation and the equation of motion. Neither equation 
changes in the presence of a reaction. Under conditions of 
local equilibrium, in order to make this a closed system of 
equations it is enough to supplement it with the system of 
equations of chemical equilibrium, which, in the case of a reac- 
tion in the binary mixture of gases, reduces to one equation, 
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Fig. 2 Coefficient of heat conduction in erg/cm-sec-deg as a 
function of the temperature for various pressures p, atm 


expressing the law of effective masses. Consequently, the 
resultant system of equations includes the following kinetic 
parameters: Lg, Lg, v’ and v”; here v’ and v” are shear and 
bulk viscosities, respectively, which enter into the expression 
for the stress tensor. In the general case the heat flow de- 
pends on two transport coefficients Z,, and Lg, but on the basis 
of the foregoing discussion it can be expressed, in the case of 


local equilibrium, in terms of one coefficient A» from Equation 
[40]. 

In the case of a boundary layer, the stress tensor terms 
proportional to the bulk viscosity drop out of the equations, 
and the system of the phenomenological equations then de- 
pends only on two transport coefficients: The shear viscosity 
v’ and the coefficient of heat conduction of equilibrium react- 
ing mixture Apo. 

Calculations of the coefficient of heat conduction in a react- 
ing binary mixture of gases (Ao) in equilibrium, performed on 
the basis of formula [40] are given in the 1955 report of the 
Laboratory of Physics of Combustion, Power Research In- 
stitute (11). See Fig. 2. 
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Reviewer’s Comment 


A calculation of the coefficients of viscosity and thermal 
conductivity of air in the temperature range 2000 to 8000 K 
and the pressure range 0.001 to 1000 atm is discussed in this 
paper. Clearly, the problem is one with important technical 
applications. It is disappointing that the present paper 
presents only rough results obtained some time ago (1955), 
and these results are given only in graphical form. 

The treatment is based on a model in which air is idealized 
as a binary mixture of “atoms” and “molecules.” The prob- 
lem of predicting the three intermolecular potentials is dis- 
cussed at some length. In the actual calculations the 
molecule-molecule interaction is taken to be of a modified 
Buckingham form and, because of a lack of information, the 
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atom-atom in atom-molecule interactions are taken to be of 
the simple rigid sphere form. The choice of the rigid sphere 
diameters is made by an interesting interpolation in the 
periodic table based on a few known values. 

The transport of energy in a reacting gas mixture is dis- 
cussed from the point of view of the thermodynamics of ir- 
reversible processes. It is concluded that, if the reaction 
rates are sufficiently rapid, the transport of energy is de- 
scribed by a single effective coefficient of conductivity. The 
calculation of this effective coefficient of conductivity is dis- 
cussed and the results are given graphically. 

—C. F. Curtiss 
University of Wisconsin 
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Methods of Combatting Interfering 
Currents That Arise at the Input of an 
Electrostatic Fluxmeter Operating 


in a Conducting Medium 


NE of the authors (3)! proposed the use of an electro- 
static rotary-type fluxmeter to measure the charge ac- 
quired by a satellite through a variety of processes, starting 
with diffusion charging due to difference in thermal velocities 
of the ions and electrons in the plasma of the ionosphere, and 
ending with charging due to ultraviolet rays. 

The operation of an electrostatic fluxmeter mounted on an 
insulated body in a plasma can be briefly reduced to the fol- 
lowing. The stator of the fluxmeter, a detailed description 
of which can be found in the monograph (1) by one of the 
authors, is a section of the surface of the body, but electrical 
contact with the remaining surface (the outer shell of the 
body) is through resistor R, which serves as the input re- 
sistance for the instrument. It is obvious that under steady 
state conditions when the rotor is stationary, the surface 
of the stator has the same potential and the same surface- 
charge density as would be possessed by that part of the 
surface which the stator replaces. Because of the periodic 
screening of the stator by the rotor, which is very close to 
the stator, the space charge, which always surrounds a 
charged surface in a plasma, causes periodic flow of surface 
charge to and from the stator. This displacement current 
flowing through the resistance R produces a voltage drop, 
which is measured with the vacuum-tube circuit. 

However, along with this displacement current, produced as 
noted by changes in charge as the stator is being screened, 
additional currents can flow to this plate in a plasma and 
thus distort the measurement. These currents result from 
the fact that the medium in which the meter is placed is 
conducting. 

Such conduction currents can be caused by the following 
factors. 

1 The satellite may move with a speed comparable to 
the thermal velocities of the ions and electrons; this causes 
current to flow to the satellite surface. 

2 The potential of the satellite at the place where the elec- 
trostatic fluxmeter is located may deviate from the equilibrium 
value, i.e., from the potential that an isolated body placed in 
the plasma would acquire at rest. This is caused by the 
motion of the satellite in the magnetic field, by the flow of 
displacement current in the input resistance, or by a forced 
change of potential in some section of the satellite surface. 

3 The satellite surface may be exposed to ultraviolet, 
x-rays, and other types of radiation capable of producing 
photoemission. 

4 Current may be produced by motion through the iono- 
sphere. 

The major distinguishing feature of the conduction current 
is that its instantaneous value is proportional to the area of 


Translated from Iskusstvennye Sputniki Zemli (‘Artificial Earth 
Satellites’), No. 3, USSR Acad. Sci. Press, Moscow, 1959, pp. 
77-83. 

1 Numbers in parentheses indicate References at end of paper. 
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the exposed portion of the stator. It is therefore 90 deg 
out of phase with the displacement current, if the flux- 
meter itself can be considered an active load (1). An esti- 
mate of the order of magnitude of the conduction current 
shows that for a charged particle concentration N = 108/cm‘, 
the limiting value of the current density does not exceed 
107 amp/cm*. At the same time, the photocurrent density 
does not exceed amp/cm?. 

The appearance of conduction current (3) at the input of 
the measuring circuit of the electrostatic fluxmeter raises the 
lower limit of the electrostatic field measured by this fluxme- 
ter, and consequently also the lower limit of the measured 
satellite charge. A way must therefore be found of combat- 
ting the conduction current. 

The fact that the conduction current is 90 deg out of phase 
with the displacement current when the fluxmeter operates as 
an active load (1) affords an opportunity of greatly reducing 
the influence of the conduction current by using a synchronous 
detector at the output of the fluxmeter measuring circuit.? 
(See Fig. 1.) 


In fact 

= Uin-k-cos 

U'out = U'in'k-cos (¢ + 90 deg) = |U’in-k-sin ¢| 

where 

Uin = effective value of the displacement signal input voltage 

Uout = signal voltage at the output of the synchronous 
detector 

U'in = effective value of the conduction input voltage 

U'out = conduction voltage at the output of the synchronous 
detector 

k = gain of the measuring circuit 


= phase shift between the reference (commutator) 
voltage of the synchronous detector and the signal 
voltage of the detector, due to the operation of 
the measuring circuit 


S; 
ynchronous 


Amplifier Detector 


Fig. 1 Measuring circuit of the fluxmeter 


2 The problem of the use of synchronous detectors in electro- 
static-fluxmeter circuits is considered in detail in (1). 
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In the ideal case, when g = 0 deg, the use of the syn- 
chronous detector would permit complete elimination of the 
conduction voltage. In an actual circuit, however, it is im- 
possible to obtain a zero phase shift in the measuring circuit 
over the entire operating frequency range; furthermore, 
variation in the temperature characteristics of the components 
used in the circuit also prevents a zero phase shift to be main- 
tained over the entire temperature range. Finally, difficulties 


t 
lu, | 


Fig. 2 Wave forms of input voltages 


a 


Fig. 3 Diagram of experiment on the determination of the ratio 
of the currents flowing on the grid and on the collector 


200 


Fig.4 Overall diagram of the experiment: 1—ignitor, 2—main 
anode, 3—probe, 4—auxiliary anode, 5—grid, 6—collector, 7— 
cathode 
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lie also in the fact that the actual wave forms of the input 
voltages (Fig. 2) contain higher harmonics, as a result of 
which it is impossible to suppress the conduction voltage 
completely. 

Let us dwell on another method of reducing the conduction 
current. We recall that in ordinary electrostatic fluxmeters 
a-c amplifiers are usually employed. Therefore the conduc- 
tion current is amplified only when it is modulated by the 
rotation of the rotor. This brings to mind the following 
method of reducing the conduction current—construct the 
rotor and stator plates in such a way as to reduce, on the one 
hand, the modulation of the flux of charged particles incident 
on the stator and on the other hand, the absolute magnitude 
of the current incident on the surface of the stator. This 
can be done by making the rotor and stator in the form of 
metallic grids of certain electrical and optical transparency. 

To verify the ability of a grid to transmit charged particles 
without stopping them, we determined experimentally the 
ratio I,/I, of the electric currents flowing in the grid and 
collector of a mercury gas discharge (Fig. 3). The general 
experimental setup is illustrated in Fig. 4. The parameters 
of the plasma were determined from the probe characteristics 
and were as follows for the condition of the experiment: Elec- 
tron temperature, 20,000 K; electron concentration, 10° to 
at a mercury vapor pressure of 1.3 107% to 5 
mm mercury. 

Fig. 5 shows the dependence of J,/7.on the thickness d of the 
space charge layer near the grid (calculated theoretically as 
a function of the potential applied to the grid and to the col- 
lector). In the experiment we used a metal grid with square 
(0.5 & 0.5 mm) meshes made of 30 » wire. As seen from Fig. 
5, at space charge layer thicknesses considerably in excess of 
the linear dimensions of the grid cells, the grid current is 
several times smaller than the collector current, meaning that 
the charged particles pass through the grid meshes without 
being stopped. This results from the fact that when the 
space charge layer is considerably thicker than the dimensions 
of the individual mesh, the distribution of the electric field at 
the grid has approximately the same character as in a parallel- 
plate capacitor. Consequently the particles entering the 
plane of the grid have a velocity component directed perpen- 
dicular to the plane of the grid, and therefore penetrate freely 


8 
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Fig. 5 Dependence of J,/IJ, on the thickness of the space 
charge layer near the grid 
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into the space behind the grid where they are trapped by the 
collector. 

The use of a grid instead of a solid rotor reduces the 
modulation of the charged-particle current, and the use of a 
grid-like stator reduces the absolute magnitude of the current. 
However, the parameters of the grid, the dimension of their 
meshes and the diameter of the wires should be chosen such 
that the electric field produced by the space charge at the 
surface does not penetrate beneath the grid; in other words, 
the grid must be equivalent to a metal plane as regards the 
distribution of the electric field at the surface. There are 
known formulas for the calculation of the electric constant of 
grids in a vacuum (2), which can be used in the first approxima- 
tion to calculate the grids of the rotor and stator for the case 
when the space charge layer is considerably thicker than the 
linear dimensions of the grid meshes. 

The construction of the receiving part of the electrostatic 
fluxmeter is shown schematically in Fig. 6, where 1 is the 
rotor, 2 is the stator and 3 is the collector. Parts 1 and 2 are 
made, as usual, in the form of multisector plates. The col- 
lector serves to trap the charged particles passing through 
the grids of the stator and the rotor. 

Rotor 1, made of wire mesh, should have a lower electric 
penetration factor than stator 2, which in turn should have a 
lower penetration factor than collector 3. Let us denote the 
distance between the rotor and stator by a2 and that between 
the stator and the collector by 23 (naturally, the spacing be- 
tween the mesh wires should be chosen less than 22 or 22). 

To calculate the electric penetration factor of the rotor 
(assuming the stator to be a metallic surface) we use the 
formula 


Ca, In coth 
= [2] 
Ca, — In cosh 70 
where 
o = filling factor of the grid 
Lg = (2/p){1 — (dg/p)] for a square-mesh grid 


To calculate the electric penetration of the stator we used 
the following formula 


R 


In coth ro 
ag = 


Ca, — In cosh To 


[3] 


The values of a; and a, determine the fraction of the source 
field that penetrates through the grids. 

Let us note that the formulas underestimate somewhat the 
values of a, and az, since the formulas can be used only when 
the source is located at a distance greater than the dimension 
of the grid mesh (2). The field produced by the space charge 
located near the grid, at distances less than the spacing of the 
grid, will penetrate much more strongly than indicated by 
the foregoing formulas. However, in view of the fact that at 
ion concentrations and temperatures observed in the plasma 
of the ionosphere the thicknesses of the space charge layers 
at negatively charged values will be on the order of a centime- 
ter, it is possible to select the wire spacing in the grid such 
as to make the error, due to disregarding the space charge 
near to the grid, less than 10 per cent. 

By way of an example, let us perform calculations based on 
formulas [2] and [3] for 212 = 0.3 cm, 223 = 1 cm, and grids 
with square meshes and the following parameters: Spacing 
0.5and 1mm, wire diam 304. The calculation data are given 
in the following 


Spacing, 
mm a2 
0.5 0.12 0.01 0.002 
1 0.06 0.05 0.01 


j 


L 
| a f 
= = 


Fig. 6 Diagram showing receiving portion of the electrostatic 
fluxmeter 


These calculations show that the grids, while almost com- 
pletely shielding the electric field, permit the passage of the 
charged particles, trapping only a small fraction of the latter. 
If the electric data are available on the circuit actually em- 
ployed, and if the coefficients a, az and b are determined (see 
material which follows), it is possible to choose the necessary 
combination of grid and collector parameters. 

The method proposed for combatting the conduction cur- 
rent due to the charged particles is equally applicable for 
combatting photocurrent, since the use of a rotor in the form 
of a grid greatly attenuates the modulation of the radiation 
fluxes that create the photocurrents. The coefficient of at- 
tenuation in this case is equal to ¢, the filling factor. In addi- 
tion, the use of a stator in the form of a grid reduces greatly 
the area subjected to the action of the radiation, and conse- 
quently decreases the magnitude of the photocurrent by a 
factor 1/cer. However, since the collector is close to the 
stator, photoelectrons may strike the latter and consequently 
gett May increase. In general, however, the conduction cur- 
rent J’paot produced at the input by a photocurrent Jphot is 

Let us consider still another way of combatting conduction 
current—the use of negative feedback proportional to the 
conduction voltage. 

Indeed, by using two detectors in parallel at the output of 
the measuring circuit, one tuned to the displacement signal 
voltage and the other to the conduction voltage, and by feed- 
ing back the detected conduction voltage to the input (4), 
it is possible to reduce greatly the amplitude of the a-c con- 
duction voltage at the input of the synchronous detector. 
The effectiveness of such negative feedback depends greatly 
on the extent to which the wave form of the input a-c voltage, 
produced by the feedback, coincides with the wave form of the 
input conduction voltage, i.e., it will depend on the purity of 
the wave form of the input conduction voltage. 

In concluding the survey of the methods that can be used to 
combat conduction current at the input of the fluxmeter, let us 
estimate the minimum value of the electrostatic field in- 
tensity, determined by the charge density on the surface of 
the satellite, which could be measured with an electrostatic 
fluxmeter without the use of additional data on its orientation. 

Let us introduce the following notation 


[4] 
u 
where 
u’ = value of minimum measurable signal voltage at 
output 
u” = value of the conduction voltage at the output 


u'in = effective value of the minimum measurable signal 
voltage at the input 
u"in = effective value of the input conduction voltage 
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The remaining symbols are indicated in connection with 
Equation [1]. 

From the theory of the electrostatic fluxmeter it follows 
that 


Emin-f-R in bR 
where 
Emin = minimum measurable value of the electric field in- 
tensity 

i = commutation frequency 
input resistance 
jo |= maximum density of conductioncurrent per unitsurface 
b = attenuation coefficient, equal to the product bib: 
b; = coefficient of attenuation of the modulation of the 


charged-particle current by the rotor grid 
bs = coefficient of attenuation of the charged-particle cur- 
rent by the stator grid 
In Equation [6], the area of the stator S is taken to be 
unity, and the entire calculation is carried out for the first 
voltage harmonic. 
Inserting the values of u’;, and u";, into Equation [4] we 
obtain 


Eninf 


Putting n = 1, we obtain for the value of the minimum 
measurable current, Emia 


Enin = [8] 


Using f = 1500 cps, 72 = 10-7 amp/cm?, and tan ¢ = 0.03, 


i.e., @ = 2 deg, we obtain several values for Emin as functions 
of the product bib. 


de 0.1 0.04 0.01 
Emin (v/em 0.54 0.22 0.05 


We emphasize that Enis is the minimum measurable elec- 
trostatic field intensity at the surface of the satellite, caused 
by its own charge, and determined under the conditions of 
maximum possible conduction current to the surface of the 
satellite. Inasmuch as the magnitude of the conduction 
current depends greatly on the orientation of the given portion 
of the surface relative to the satellite course, the values of 
Emin also depend on the same parameters. 

For comparison with Emin, we list the values of the in- 
tensity of the electric field on the surface of a stationary 
sphere, produced by its own charge under the following 
plasma parameters: N = 10 and 10° cm’, T; = 1000 K, 
T. = 5000 and 10,000 K (we assume that the radius of the 
sphere is considerably greater than the thickness of space 
charge layer around it). 


T, = 5000 K T, = 10,000 K 
N = 108 1.0 i2 
N = 105 0.32 0.38 
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Reviewer’s Comment 


The senior author of this excellent article has been as- 
sociated with the USSR program in atmospheric electricity 
for more than a decade. He is the leading Soviet proponent 
of the use of the rotating electrostatic fluxmeter, or electric 
field meter, to measure the electrical charge acquired by a 
satellite in the ionosphere. 

The description of the operation of the fluxmeter is concise 
and conventional. The 90-deg phase displacement between 
the conduction and displacement currents was also recognized 
by American scientists in the upper atmosphere rocket 
research program at least as early as 1947 [Clark, 1949 (2)], 
although no detailed description of the theory of operation 
of this mechanical-commutator instrument was published. 
Later versions of this American device [Clark, 1957 (3)] used 
an electromagnetic commutator and synchronous detector 
very similar to the Soviet system. British workers have 
utilized the same quadrature phase relationship to eliminate 
precipitation current effects from measurements of the electric 
field strength at the Earth’s surface [Mapleson and Whit- 
lock, 1955(4) ]. 

For some unaccountable reason, Imyanitov and Shvartz 
apparently do not use a variable phase adjustment in their 
commutator voltage source. This American technique 
gains a factor of about 10 below the USSR value of 0.03 for 
tan ¢. 

But the core of the present Soviet contribution is the 
ingenious metallic mesh construction of both the rotor and 
stator assemblies, which simultaneously decreases the photo- 
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current, the ion current, and the modulation of both currents, 
while maintaining high efficiency for the electric field meas- 
urement. As shown on the table following Equation [8], 
the factor of improvement varies from 10 to 100, depending 
on the percentage of open area of the mesh. 

Combining these two factors, one favorable and the other 
unfavorable to the USSR, we find that the best Soviet flux- 
meter is capable of measuring an electric field from six to 
ten times smaller than that which can be measured by the 
best solid-plate American model, for the same magnitude of 
interfering conduction current. 

It should also be noted, from a comparison of the last two 
Soviet tables, that their best instrument is capable of meas- 
uring an electric field strength some six times smaller than 
that expected to occur at the satellite surface. This means 
that the American device is still adequate to make the 
measurement under these worst possible assumptions of 
electron density, satellite orientation and conduction current 
density [Bourdeau and others, 1960 (1)]. However, there is 
little doubt that the Soviet fluxmeter is at present potentially 
the best available for its purpose. 

—Joun F, 
NASA 
WASHINGTON, D. C. 
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" Some Results of the Determination 
| Of the Structural Parameters of the 
“ | Atmosphere Using the Third Soviet 
. Artificial Earth Satellite 


ace The problem of measuring the pressure and density of high layers of the atmosphere using an 
artificial Earth satellite was discussed in (1). An analysis is given in this work of the state of the 
ideas held prior to 1957 regarding the structural parameters of the upper atmosphere. By that time 
: the mean distributions of pressure, density and temperature up to a height of 100 km were known 
(2-4), and it had been established that at heights up to 90 km the atmosphere is intermixed, and that 
oxygen dissociates above 90 km (5). Up to 1956, direct measurements of the density, pressure and 
composition of the atmosphere above 100 km were very few in number, and therefore the ideas con- 
cerning the atmosphere at these altitudes were not well defined. In recent years experiments have 
been carried out on the determination of the density of the atmosphere at great altitudes (6-8). A 
particularly large contribution to the study of the upper layers of the atmosphere resulted from the 
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investigations made with Soviet artificial Earth satellites, which made it possible to determine the 
sign density of the atmosphere both from measurements with manometers installed on the third Soviet 
oon Earth satellite (9) and from the drag on the satellites (10-11). The present article is devoted to an 
sical analysis of a part of the results of the determination of the density of the atmosphere obtained from 
“— manometer measurements on the third Soviet artificial satellite. 
I Apparatus mum vapor In manometers 
iia oe ea in the satellite were first outgassed and evacuated. The ma- 
HE APPARATUS (Fig. 1) used on the artificial Earth nometers were uncapped with the aid of a breaking device dur- 
satellite to determine the density of the atmosphere in- ing the separation of the protective cone from the satellite 
cluded two ionization manometers (M; and M2) and one mag- after it wan in orbit 
netic electric discharge manometer The command for 
-" switching on the apparatus was transmitted from the program 
ii. timing mechanism (PTM) to the switching-connecting device II Method 
[8] (SCD). The current between the electrodes of the manom- an — 
line eter, proportional to the number of particles per unit volume, The density and pressure within a manometer chamber are 
g was fed to the input of a direct current amplifier (DCA). not equal to their values in the atmosphere, owing to the high 
7 The output of the amplifier, which depended on the size of the satellite velocity, temperature differences between the gas in 
naa input signal, was transmitted to Earth by a radiotelemetry the satellite and the free atmosphere, as well as changes of the 
=. system (RTS). orientation of the inlet aperture of the satellite relative to the 
: The position of the manometers on the satellite is shown in velocity vector of the satellite. 
A Fig. 2. The manometers were installed in conical recesses in At the altitudes of interest to us, the mean free path is very 
“e the head part of the satellite, in so-called sockets. The ma- much greater than the characteristic dimension of the appa- 
nometers were secured in the sockets with the aid of panels ratus (the satellite). Under these conditions there are no 
= equipped with hermetic contacts, through which the elec- collisions between particles in the inflowing stream and the 
~ trodes of the manometers were connected with the ampli- particles that remain in the manometer. If the velocity of a 
sone fiers and power supplies inside the satellite. body vw is significantly larger than the mean thermal velocity 
“the To eliminate the effect of atmospheric ions and electrons on of the particles in the medium, then NSu sin @ particles will 
igo the manometer readings, special traps were installed at their pass in a unit time through a manometer opening of area S (N 
aod inputs (1,9). Besides this, a screen in electrical contact with is the number of particles per unit volume of the atmosphere, 
as the body of the satellite was placed at the edge of the socket and @ is the angle between the velocity vector and the plane 
i. protecting the collecting outlet of the manometer from incident of the aperture). Undergoing a series of collisions with the 
ially charged particles. walls of the manometer, the particles attain a most probable 
— To decrease the effect of desorbed molecules on the ma- velocity 1, which is determined by the wall temperature 71. 
ASA nometer readings, the outer shell of the satellite was made of In one unit of time (Nin/2V; n)S particles flow out of the 
We materials having a maximum rate of desorption and a mini- manometer aperture. At equilibrium both particle fluxes are 
8 equal. From the equality of the fluxes, the concentration of 
rom Earth particles Nata given point in the atmosphere is given bythe 
1360. | 97- Translation provided by the U. S. Joint Publication Re- formula 
search Service. 
1 Numbers in parentheses indicate References at end of paper. Na,-P; 
ol. 7, 2 The construction and operating principle of the manometers N= (0>7-8 deg) 
are described in detail in (9, 12 and 13). V 2rMR*T\u sin 0 
| ApriL 1960 407 
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Fig. 1 Block diagram of the manometric apparatus mounted on the satellite 


provided the manometer is graduated not in particle concen- 
tration N, but in pressure P; = N,kT,. Here M is the mo- 
lecular weight of the gas, R is the universal gas constant, and 
N,, is Avogadro’s number. 

The velocity u, which is needed for the calculation of N, is 
known with great accuracy from observations of the orbit. 
T; is measured directly on the satellite and is telemetered to- 
gether with the manometer reading P,. 

The molecular weight was taken equal to that given in (14); 
see Table 1. The value of M above 600 km was obtained by 


extrapolation. : 
re seen from the formula presented, a possible change in Fig. 2 Placement of the manometers on the satellite 
the molecular weight does not lead to a significant error in the b’ = —sin@dsing + cos dcos¢ cos % 
determination of N. For example, a change of M from 18 to c’ = cosgsin & 
24 gm/mole causes N to decrease less than 20 per cent. a” = sin dain & 
The orientation of the inlet aperture of the manometer rel- b” = cos sin & 
ative to the velocity vector of the satellite (sin #) was deter- c” = cosh 


mined from a method established in (15) 


sin = m {v [Aya + + Cyc] + [Ara + Bob + Cac] + + Bab + Circ) } + n{v[ — Ara’ + Bib’ + Cye’] + w Ara’ + 
+ Crc'] +£[ — Asa’ + + Cyc']} + & {v [A.a” — Bib” + [Asa” — Bob” + + & [Asa” — + 


where 
Here ¢, y and are the Eulerian angles; v, uw and &é are the 
a@ = cosdcosy — sin dsin ¢ cos % direction cosines of the velocity vector in an absolute system 
b = sin dcos¢y + cos dsin gy cos & of coordinates (the Z axis is directed toward Polaris, the X axis 
c =singsin & toward the vernal equinox, and the Y axis completes a right- 
a’ = cos dosing + sin cos ¢ cos | handed coordinate system) as in Fig. 3; A,;, B;and C; (¢ = 1, 2, 
Table 1 Variation of molecular weight with height 
Molecular Molecular Molecular Molecular 
Height, weight, Height, weight, Height, weight, Height, weight, 
km gm/mole km gm/mole km gm/mole km gm/mole 

225 21.28 295 18.63 365 17.25 435 16.53 

230 21.04 300 18.50 370 17.18 440 16.50 

235 20.80 305 18.37 375 17.12 445 16.46 

240 20.57 310 18.25 380 17.06 450 16.43 

245 20.36 315 18.14 385 17.00 455 16.40 

250 20.15 320 18.03 390 16.94 460 16.37 

255 19.94 325 17.92 395 16.89 465 16.34 

260 19.75 330 17.82 400 16.84 470 16.31 

265 19.56 335 17.73 405 16.79 475 16.30 

270 19.40 340 17.64 410 16.74 480 16.27 

275 19.22 345 17.55 415 16.70 485 16.24 

280 19.07 350 17.47 420 16.65 490 16.22 

285 18.91 355 17.39 425 16.61 495 16.19 

290 18.77 360 17.32 430 16.57 500 16.16 
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Table 2 Structural parameters of the atmosphere at heights of 225-500 km 

Height, 
km N, P, gm/cm* H, km P, dyne/cem* P,mm Hg 
225 6.01-109 2.12-10-18 40.0 936 7.76-10~4 6.25-107 
230 5.31 1.79 40.6 938 6.88 5.54 
235 4.7 LF 41.3 941 6.1 4.92 
240 4.17 1.42 42.0 946 5.44 4.4 
245 3.71 1.25 42.8 952 4.88 3.94 
250 3.3 i | 43.5 958 4.36 3.54 
255 2.94 9.73-10-%4 44.3 964 3.91 3.17 
260 2.64 8.66 45.2 971 3.54 2.88 
265 2.36 7.80 46.0 979 3.19 2.6 
270 2.12 6.83 47.0 987 2.89 2.35 
275 1.91 6.1 47.9 996 2.63 2.14 
280 1.72 5.44 48.8 1005 2.39 1.95 
285 1.55 4.87 49.7 1015 2.17 1.78 
290 1.4 4.36 50.7 1026 1.98 1.62 
295 1.27 3.93 51.7 1037 1.82 1.49 
300 1.15 3.53 52.7 1048 1.66 1.37 
305 1.07 3.26 53.7 1059 1.56 1.29 
310 9.57-108 2.9 54.5 1072 1.42 117 
315 8.73 2.63 55.9 1084 1.31 1.008 
320 7.98 2.39 57.0 1097 1.21 1.0 
325 7.31 2.17 58.1 1110 1.12 9.28-10-3 
330 6.7 1.98 59.2 1124 1.04 8.62 
335 6.17 1.82 60.3 1138 9.69-10-5 8.06 
340 5.68 1.66 61.5 1153 9.04 7.52 
345 5.22 1.52 62.8 1169 8.96 7.46 
350 4.82 1.4 64.8 1185 7.88 6.58 
355 4.46 1.29 65.2 1200 7.39 6.18 
360 4.13 1.19 66.7 1219 6.95 5.82 
365 3.86 i a 68.1 1237 6.59 5.53 
370 3.56 1.02 69.5 1257 6.18 5.19 
375 3.31 9.41-10-% 70.9 1276 5.83 4.9 
380 3.08 8.72 72.4 1295 5.51 4.64 
385 2.92 8.24 73.9 1305 5.26 4.44 
390 2.69 7.56 75.2 1335 4.96 4.19 
395 2.52 7.07 76.7 1353 4.71 3.98 
400 2.36 6.6 78.9 1373 4.47 3.79 
405 2.21 6.16 79.7 1393 4.25 3.6 
410 2.08 5.78 81.2 1417 4.07 3.46 
415 1.95 5.41 82.9 1440 3.88 3.3 
420 1.84 5.09 84.6 1465 3.72 3.17 
425 1.73 4.79 86.3 1489 3.56 3.04 
430 1.64 4.51 88.1 1514 3.43 2.93 
435 1.55 4.25 90.0 1539 3.29 2.82 
440 1.47 4.03 91.7 1563 3.17 2.72 
445 1.39 3.8 93.6 1589 3.05 2.62 
450 1.32 3.6 95.5 1614 2.94 2.53 
455 1.25 3.4 98.6 1643 2.84 2.44 
460 129 3.23 99.9 1675 2.75 2.37 
465 1.13 3.06 02.0 1709 2.66 2.3 
470 1.08 2.92 104.5 1745 2.6 2.25 
475 1.03 2.79 107.0 1781 2.53 2.19 
480 9.82-107 2.65 109.3 1810 2.45 2.13 
485 9.4 2.53 111.5 1845 2.39 2.08 
490 8.97 2.42 113.9 1880 2.33 2.02 
495 8.61 2.31 116.5 1917 2.28 1.98 
500 8.24 2.21 119.0 1953 2.22 1.94 


3) are constants for a given loop 


A, = COS Pp COS Yo A? = — sin py As = COS Pp BIN Yo 
By = sin Ye B, =0 Bs; = — cos Yo 
C; = sin po Cos Yo C: = cospo C; = sin p sin Yo 


po and Yo are angular coordinates of the angular momentum 
vector; m, n and k are the direction cosines of the manometer 
axes relative to the principal axes of the satellite. 

After determining N (the number of particles per unit 
volume), the density p was calculated and also the height H of 
the homogeneous atmosphere*® 


Here N” and N’ are the number of particles per unit volume of 
the atmosphere at two points differing in altitude by Ah = 10 
km 


Next were calculated the temperature 
T = MGH/R* 


and the atmospheric pressure 
P=kNT 


when k is the Boltzmann constant. 


’ The value of H determined with this formula leads to an error of about 30 per cent compared with calculations based on the formula 


Ah 
~ (log N*/N’ — log T°/T’) 


Aprit 1960 


Ah 
2.3 (log N’ — log N’) 


e the 
stem 
axis 
right- | 
1, 2, | 
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It should be noted that in the determination of the density 
of the atmosphere we did not use the hypothesis that the local 
distribution of the particle velocities is Maxwellian at great 
heights. The height of the homogeneous atmosphere, the 
temperature and the pressure were calculated under the as- 
sumption that the usual barometric formula is appplicable to 
the rarefied gas in the presence of temperature gradients. 


III Results 


The third Soviet Earth satellite was launched on May 15, 
1958. In accordance with the plan of the experiment, from 
the moment that the manometers were uncovered their dis- 
charge currents (corresponding to the density measurements) 
were continuously recorded, and periodic calibration of the 


Table 3 Values of the density (in gm-cm~*) at various heights from manometer measurements on rockets and from the dragon | 
satellites 
Drag data = 
— Rocket data 1957 a, 
Containers 1957 ae 
and rockets. 1957 a; 1957 By 
Mean latitude 1957 az 1958 a 
Height, of European Viking-7 Aerobee-Hi 1957 Bi 1958 Be 
km USSR 33° N lat. 59° N lat. (10, 11) 1958 + (19-25) 
100 4-10-” 2.5-10-” 7-10- 
110 9.8-10-1! 5-10-11 1.5-10-” 
120 2.2-10-" 6.9-1071! 
130 7.4-10-% 6-10-12 3.01-107"! 
140 3.2-10-12 1.2-10-12 3-10-12 1.49-10-"! 
150 1.6-10-% 2-10-12 8.07-10-12 
160 9.5-107-18 4.3-1078 4.70-10-12 
170 6.4-10-8 3.0-10718 1-10-12 2.89-10-2 
180 4.4-10738 2.3-1078 8-10-38 1.87-10-12 
190 3.3-10738 1.8-10718 7.5-10738 1.25-10-12 
200 2.7°10-# 1.4-10738 7 (3-8.63)-10—18 
201 +4 6.7-1078 
202 + 4 7.37°10-8 | 
206 + 4 5.4-10718 
210 2.0-1073 6.04-10738 
212 (4.4-4.8)-10-18 
215 4.7-10-43 
220 1.6-10-38 9.0-10-14 (3.5-5.7)+10-2 
230 1.25-10738 3.32-10—38 | 
232 1.5-10733 | 
233 2.2-10-8 
240 1.1-10-38 2.51-10—"8 
241 2.5°10733 
250 9-10-14 (1.5-1.6)-107—3 (1.1-1.9)-10718 
260 6.9-10-"4 1.51-10-8 
270 (9.4-10)-10-14 1.19-10-38 
275 8.5-10-14 
280 9.51-10-4 
290 (5.8-7.0)-10-4 7.68-10-"4 
300 (5-6 .27)-10—14 
310 (3.8-4.7)-10-"4 5.16-10-"4 
320 4,29-10-14 
330 (2.6-3.2)-10-™4 3.58-10—14 
340 3.02-10-%4 
350 (1.8-2.2)-10-"4 (2.1-3)-10-14 
360 2.18-10-14 
368 (1.4-1.5)-10-" (1.4-1.5)-10-"4 
370 1.87-10-4 
400 1.5-10-"4 
9.3-10-% 
405 (978)-10-% 
450 (1.0-4.5)-10-% 
500 (2.3-6)-10-% 
550 (2.2-4)-10-% 
600 2-10-* 6.8-10-% 
650 1-10-45 3.8-10-% 
700 7-10-16 
720 (1.2 + 0.3)-10-% 
Conventional designations of the satellites: 1957 a;—Carrier rocket of tue first Soviet satellite (launched Oct. 4, 1957); 1957 ag— 
First Soviet artificial satellite (launched Oct. 4, 1957); 1957 a;—Nose cone of first Soviet satellite (launched Oct. 4, 1957); 1957 
B,:— Record Soviet artificial satellite (launched Nov. 3, 1957); 1957 B.—Nose cone of record Soviet satellite (launched Nov. 3, 
1957); 1958 a—Explorer satellite (launched Feb. 1, 1957); 1958 B:—Carrier rocket of Vanguard (launched March 17, 1958); 
1958 B-—Vanguard I satellite (launched March 17, 1958); 1958 V—Explorer III satellite (launched March 26, 1958). 
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Fig. 3 Illustration of the orientation of the manometer in space 


amplifiers was made, as were checks of the emission current 
and the manometer wall temperature. 

The results obtained provide a description of the atmos- 
phere at altitudes of 225-500 km for May 16, 1958. Meas- 
urements at various altitudes refer to different times of day 
(13:00-19:00 local time) and different geographical latitudes 
(57 deg N-65 deg N). 

The values of the density obtained in the present experi- 
ment practically agree, in the region of perigee of the Soviet 
satellites (225-230 km), with the density values found from 
observation of the change of orbit (see Fig. 4a and Table 3). 

At the same time, the values of density and pressure ob- 
tained are about 1.5-10 times greater than the values proposed 
by numerous models (14,16,17), with the exception of the 
model introduced in (18). See Figs. 4b and 5. In the major- 
ity of cases the models were based on generalizations of a very 
few observations on the upper layers of the atmosphere, ob- 
tained by various means, at different times of day and sea- 
son, at different latitudes, etc. The experimental investiga- 
tions of the upper layers of the atmosphere were isolated. 
The assumptions used in the construction of various models 
were distinct and often contradictory, for instance, the use of 
diffusion separation in some and intermixing in others. All 
this causes the models to be noticeably different from each 
other, so that they neither characterize the actual state of 
the atmosphere (since recent investigations indicate the pres- 
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Fig.4 Variation of density with altitude. (a) Experimental data; 
(b) based on several models. 


ence of seasonal, daily, latitude and other variations, and in- 
dicate that it is incorrect to unify obtained data) nor yield 
statistical average values of the parameters, since the amount 
of data is insufficient for this. 

Some of the differences in the data obtained from the ma- 
nometer readings and from the satellite drag can be explained, 
obviously, by the presence of regular daily, seasonal and lati- 
tude variations as well as variations associated with solar dis- 
turbances. 

It is known that the drag of the second and third Soviet 
satellites was not constant (27-30). In January 1958 the 


Number of revolutions (from Nov. 4, 1957) 
200 400 G00 800 1000 1200 
T T T T T T T 


-50 
-48 


> 
y 

r 


> > 
+ 


T 


- 30 

Oec.1,1957 Janl, 1958 
Fig. 6 Variation of period of the second Soviet satellite 

411 


Rate of change of period, sec/day 


1958 


15, 
rom = 
nts) ) 
the 
| 
| \ 
20. 
x \ Third Soviet Satetlite 
| 400 Gerson, 1957 (summer) 
| \ 
| Mitra 
| 1952 ™. 
| 
| 
| 
\ | 
! 
j 


rate of decrease of the period of the second satellite fluctuated 
between 4.4 and 3.9 sec per day (Fig. 6). In February an in- 
crease of from 3.9 to 5 sec per day alternated with a decrease 
of as much as 4.5 sec per day. The probable error in the de- 
termination of the period was found to be about 0.03 sec per 
day (27). 

In the first days of the existence of the third satellite, the 
decrease in the period of rotation amounted to 0.97 sec per 
day, and after 450 revolutions (June 17) to only 0.59 sec per 
day (30). A retardation in the decrease of the period is also 
noticeable in a curve of the change of period with time (Fig. 
7). Each point represents an average taken over 54 revolu- 
tions. It is noted in (30) that in the period from May 15 to 
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Fig. 7 Change of the period of rotation of the third Soviet 
satellite 
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Fig.8 Correlation between the rate of change of the period of ro- 
tation of the second satellite and the intensity of stellar flares 
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Fig. 9 Variation of density of atmosphere with latitude 
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July 17 the atmosphere contracted so to speak by an amount 
equal to 25 km. 

In addition, during observations of solar activity and of the 
change of the period of rotation of the second Soviet satellite, 
it was observed that the strong magnetic storm in the middle of 
February 1958 and the 10 deg extension of the zones of the 
aurorae toward the Equator were accompanied by an increase 
in the change of the period of the satellite from 3.9 to 5 sec 
per day (31). Changes in the intensity of flares on the sun 
(i.e., the product of the area of a flare by its duration) also 
exhibit a correlation with the rate of change of the period of 
the satellite (32). See Fig. 8. 

The results of rocket investigations point to the existence 
of daily, seasonal and latitude variations (33,8). Thus, ac- 
cording to the data of these articles, at a height of 200 km 
(Fort Churchill, 59 deg N) the summer daytime value of the 
density (6.7 + 2.0) X 10-"* gm/cm‘ is approximately twice as 
high as the corresponding winter value 


(36 X 10-13 gm/cm? 
and five times larger than the winter night-time value of 
the density (1.3 + 0.6) X 10-% gm/cm~*. At the same 
height the density at latitudes 33 and 59 deg N on a 
summer day proved to be equal to (1.4 + 0.5) X 107! and 
(6.7 + 2.0) X 10-8 gm/cm', respectively. The presence of 
a latitude effect is also shown by the results of investigations 
made with meteorological rockets (34), by the method of fall- 
ing spheres (35), and in other ways. 

The available experimental data are insufficient to provide 
the hoped-for determination of the functional dependence of 
density on latitude in the upper layers of the atmosphere (Fig. 
9). However, all the results of rocket density measurements 
made with manometers on summer days (curve 1) display an 
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Fig. 10 Time variation of the period of rotation of the second and 
third satellites and the rocket of the third satellite with time 
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Fig. 11 Temperature of the atmosphere 


increase in density with increasing latitude in the range 33- 
59 deg N. For winter days, however, there is observed a 
somewhat less pronounced increase of density with latitude 
[see (33) and Table 3, column 2]. 

The latitude effect on density in the fall and winter season, 
as gathered from observations of the change in orbit of the 
first and second Soviet satellites, is much weaker [see curve 2 
in Fig. 9 and (36,37)], possibly because the same value for the 
height of the homogeneous atmosphere, H = 30 km, was used 
for all latitudes in extrapolation of the data to 206 km. In 
addition, it is possible that these data refer to different times 
of the day. 

The observed kink in the curves of the periods of rotation of 
the second and third satellites and the carrier rocket of the 
third satellite vs. time as the latitude 30 deg is approached 
probably indicates an increase in density on approaching the 
Equator, owing to the equatorial bulging of Earth. 

The facts enumerated attest to the fact that the state of 
the upper atmosphere may vary markedly and that it depends 
both on external effects on the atmosphere and also on the 
season, time of day and geographical latitude. 

Values obtained from measurements on the third Soviet 
artificial Earth satellite, which give the density, height of the 
homogeneous atmosphere, temperature and pressure, char- 
acterize the state of the atmosphere at altitudes of 225-500 
km for the daytime hours on May 16, 1958 at latitudes of 57- 
65 deg N. In particular, on May 16, 1958 the temperature of 
the atmosphere appeared to be relatively high (Fig. 11). 


As a result of further reduction of all the experimental data, 
additional material will be obtained on the density, pressure 
and other parameters of the atmosphere, which will enable a 
clearer explanation of the dependence of the structural pa- 
rameters on the latitude and time of day. 
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Reviewer’s Comment 


This article is a competent account of experimental pro- 
cedures and results obtained by USSR scientists in the 
measurement of upper atmosphere densities, using manom- 
eters in Sputnik III as well as in the analysis of the drag 
on the Sputniks. The discussion of Soviet work is completed 
by a general review of density determinations by rockets and 
satellites, and their theoretical interpretation, both in the 
U.S. and in the USSR. 

The values of density and pressure quoted in the USSR 
paper agree with the results of other investigations in their 
indication that upper atmosphere densities are greater by an 
order of magnitude than the values predicted in earlier at- 
mospheric models. Among the interesting aspects of the 
USSR results is the indication for a significant latitude de- 
pendence of density in rocket manometer measurements. 


Apri 1960 


Mikhnevich et al. also report the detection of a latitude 
effect by the measurement of the drag on Sputniks I and II, 
although most U. S. and British analyses of satellite drag 
have not shown a significant latitude variation. It is clear 
that additional data, preferably obtained in situ by gages 
rather than by study of satellite orbits, will be required 
for the resolution of this very interesting question re- 
garding the degree of latitude variations in the upper at- 
mosphere. This reviewer is in agreement with the position 
taken by the USSR investigators, who conclude that the 
“state of the atmosphere may vary markedly and that it 
depends on external effects on the atmosphere and also on the 


season, time of day, and geographical latitude.” 
—RosBertT JASTROW 


NASA 
Washington, D. C. 
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Motion of a Slender Blunted Body 
in the Atmosphere 


With High Supersonic Speed 


E CONSIDER the question of the influence on the dis- 

tribution of pressure over a thin blunted body (i.e., a body 
with a longitudinal dimension much greater than the dimen- 
sion of the blunted edge) of the strong increase in entropy of 
the gas behind the nose shock wave arising from the blunted 
leading edge. The question of the effect of the blunted edge 
on the distribution of pressure, without allowance for the 
foregoing phenomenon, was considered by G. G. Chernyi 
(1,2).2 


1 Let p p.V.*, pew, and wV,, be the pressure, 
density, enthalpy and absolute velocity of the gas. The sub- 
script © refers to dimensional quantities in the fore stream. 
With a constant adiabatic exponent « on the streamlines which 
intersect the shock wave at an angle ¢, and with M sin g>1, 
the flow parameters will be connected with p by means of the 
following relations? 


1/« 


w= V1 2i 


In our case the angle ¢ will range from 7/2 to small values 
of the order of the relative slenderness of the body. At a dis- 
tance of the order of the diameter of the blunted edge, down- 
stream of the nose of the body, the pressure between the 
surface of the body and the shock wave is of the same order 
of magnitude as the pressure behind the shock, and p < 1 
for slender bodies. From this we obtain from Eq. [1.1] 

% 


~ (sin g)~ 
Po 
The subscript 0 refers to the flow quantities which have 
passed through the normal portion of the shock. The rela- 
tions [1.2] show that at the surface of the blunted slender 
body there is formed a high entropy layer of gas,* in which the 
flow quantities have a different order of magnitude from those 
in the remaining region of flow, where sin ¢ < 1. Near the 
surface of the body the condition w ~ 1, which holds in hyper- 


Translated from Izvestiia Akademii Nauk SSSR, Otdelenie 
Tekhnicheskikh Nauk, Mekhanika i Mashinostroenie (Bull. USSR 

cad. Sct., Div. Tech. Sci., Mechanics and Machine Building), 
no. 4, 1959, pp. 131-133. 

* Numbers in parentheses indicate References at end of paper. 

2 The limitation M sin g >>1 is used here only for simplicity; 
the phenomenon discussed in the following takes place for slender 
bodies even for M $ 4. 

* A detailed examination of the flow in a high entropy layer for 
constant « and variable flight conditions is given in (3). 
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sonic flow for sharp slender bodies (law of plane sections), is 
not satisfied. 

For flight in the atmosphere with M > 10, the air in the 
high entropy layer is in a dissociated state. For qualitative 
deductions we can use here the visual constant adiabatic 
exponent formulas 


Ko 


etc. (p ~ 1 at the critical point) with a smaller adiabatic ex- 
ponent 


Ko = 1.15 + 1.3 


On a blunted cone of half-angle a = 10 deg at M = —, 
at a large distance from the nose and with x = 1.4, 1.3 and 
1.2, the values of % and wo are 0.186, 0.224, 0.28 and 0.79, 
0.745 and 0.665 respectively. Thus, as xo is decreased, an 
ever-increasing part of the work done by the motion of the 
blunted nose will remain in the high entropy layer in the form 
of irreversible thermal energy and in the form of kinetic 
energy of the translational motion of the gas in the direction 
of the flight of the body. In accordance with the principal 
premises of references (1, 2 and 4) it follows from the fore- 
going arguments that the dissociation of the air in the high 
entropy layer should decrease the effect of the blunting on the 
distribution of pressures along the lateral surface of the body. 

2 We shall confine ourselves henceforth to a consideration 
of axisymmetric bodies. Let ro be the radius from the axis 
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of the blunt nose; xr» and yr the Cartesian coordinates along 
and normal to the axis of the body, with the origin of the 
coordinates at the nose; ror, ro and reyo the shape (in the 
meridian plane) of the surface of the body, of the shock wave, 
and of the boundary of the high entropy layer respectively; 
and V ,u the velocity component in the direction of the x axis. 

At large values of M the shock wave has a large inclination 
only at points which lie opposite points of the surface with 
large slope. Therefore, by virtue of the remarks made previ- 
ously, we can assume that the mass flow of gas in the high 
entropy layer equals the mass flow of the incident stream 
passing through an area equal to the area of the blunted por- 
tion from the axis, and that outside of this layer u ~ 1, i.e., 
/ = Yo is the line of contact discontinuity. Outside the high 
entropy layer we consider the gas to be perfect with a constant 
x. As in (2) we assume that the pressure in the disturbed 
region, with the exception of a negligibly small region adjacent 
to the shock wave, is approximately equal to the pressure 
p = p(x) on the surface of the body, and that the velocity of 
the gas in the direction of the y axis equals its value behind 
the shock wave 


o(z) = 2R'(M*R’? — 1)/[(K + 


Such a rough approximate formulation of the problem is 
fully justified in estimating the influence of the foregoing ef- 
fect. Applying the momentum and energy theorems to a 
layer of unit width fixed in space (perpendicular to the direc- 
tion of motion), and taking into account the added work done 
by the pressure forces, mass transfer, momentum and energy 
transfer through this layer due to the longitudinal velocity of 
the gas, 1 — u, we obtain the following system of equations 
for the motion of the gas 


1 x 1 
+ (1 R'R I+ (» an) Rdx [2.1] 


1 
K 1 
+1 M?R"” 


R? x 
E-E, 
+ K(k — DM? +2 f, r'rpdx 


il — 


2 p U4 


E = (1/2)er) E, = 


Here E and J are the dimensionless energy and momentum 
(per unit angle) along the y axis, imparted to the gas by the 
blunt nose (2), and cz is the drag coefficient of resistance of 
the blunt nose. 

Equations [2.1 and 2.2] differ from the corresponding equa- 
tions obtained by G. G. Chernyi (2) in the inclusion of the 
term E,(p) which, as shown by analysis, is always positive. 


The value of £,(p) is comparable with the value £F if the air 
in the high entropy layer is dissociated, and is negligibly small 
if the gas is assumed to be everywhere polytropic with k = 
1.4. 

In order to make the system [2.1 and 2.2] closed at x = 0, 
we assume that the density of the gas in the regon yo < y < 
R is equal to the density behind the shock p; atz = 0. Then 


1 (x + 1)M?R” 
R Yo (R ) 2+ (K 1)M?2R” 


[2.3] 


pu 


From [2.1 to 2.3] we determine the initial values of p, R 
and R’ at z = 0. 

The figure shows the results of calculations of the pressure 
distribution p over cones with a = 10 deg, blunted in the 
form of a hemisphere (curve 1) or flatly truncated (curve 2), 
flying at M = 30 at an altitude of the order of 40 km (solid 
lines) and in a constant property gas with «x = 1.4 (dotted). 
The values cz = 0.92 and 1.8 were taken for the spherically 
blunted and flat faced nose respectively. The same figure 
shows curve 3, obtained at M = © by G. G. Chernyi by a 
somewhat different formulation (2). We note that almost 
the same results are obtained by solving the system [2.1 + 
2.2) at M = 30and at M = ~, 

In the determination J, E; and yo, the form of the shock 
wave arising from the blunted leading edge and the form of 
the blunt surface are assumed parallel. The distance between 
the shock and the body was determined from the formulas of 
(5), and the pressure was determined assuming a Newtonian 
pressure distribution. The values of the parameters for dis- 
sociated air were taken from the tables of (6). 

The results given are evidence of the fact that in the atmos- 
phere at very large flight velocities with thin blunt bodies 
(M > 20, for a = 10 deg, as shown by calculations) the in- 
crease in entropy behind the shock wave arising from the 
blunted leading edge noticeably affects the pressure distribu- 
tion at large distances downstream. 

The author expresses his sincere gratitude to T. G. Korchi- 
kova and A. G. Kuz’mina for much help in the performance 
of the calculations. 
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Reviewer’s Comment 


The problem treated in the present translation is that of 
determining the pressure distribution on a thin blunted bedy 
in a hypersonic flow. In the United States, Cheng (1), and 
in Russia, Chernyi (2), have shown that with small blunting 
the hypersonic equivalence principle [or hypersonic similitude; 
see Hayes and Probstein (3)] can be applied, provided an ad- 
ditional similitude parameter which depends on the bluntness 
is invariant. In Russia the equivalence principle is referred to 
as the “law of plane sections.” However, with blunting the 
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equivalence principle fails in the vicinity of the blunt nose 
itself, as well as in a layer of fluid of high entropy near the 
surface of the body. Previous works have used the equiv- 
alence principle [e.g., Lees and Kubota (4), and (1, 2 and 4) 
of the present translation] in order to apply one-dimensional 
nonsteady inviscid solutions to the calculation of steady 
hypersonic flows over slender bodies with small nose blunt- 
ing. In so doing the influence on the solution of the high 
entropy layer near the surface in which the equivalence prin- 
ciple fails was neglected. The importance of this entropy layer 
has previously been pointed out by Hayes and Probstein (5). 
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Lunev in his paper has investigated the effect of this entropy 
layer by assuming the shock layer to be divided into two 
regions by a discontinuous surface across which the longi- 
tudinal velocity can be discontinuous. The region between 
this contact surface and the body surface is termed a “high 
entropy layer.” This layer results from the high entropy 
which arises behind the normal portion of the shock wave 
generated by the the blunt nose. Lunev assumes that in the 
region outside of the entropy layer, between the contact sur- 
face and the shock wave, that the flow variables are constant 
and have values equal to their values just behind the shock 
wave. He further assumes that the pressure and normal 
velocity are constant across the entire shock layer. The only 
question which arises here is that if the normal velocity is 
constant across the shock layer and depends only on the 
longitudinal coordinate, then if there is a discontinuity in the 
longitudinal velocity across the edge of the entropy layer one 
is led to the conclusion that there should be also a discon- 
tinuity in flow direction at the edge of the entropy layer. 
It would therefore seem that one should also take into con- 
sideration the normal velocity change across the shock layer. 


However, it is possible that this effect may turn out to be 
small. 

Under Lunev’s assumptions he derives the appropriate 
momentum and energy integrals for the flow taking into ac- 
count the energy contained in the entropy layer. Using rough 
estimates for the shock shape, for the detachment distance at 
the nose, for the pressure distribution on the nose, and for 
the mass flow in the entropy layer he calculates a modified 
pressure distribution along the body surface which differs 
considerably from that obtained without consideration of the 
entropy layer. This reviewer believes these results to be 
very interesting and important. 

F. Prospstern 
Brown University 
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Investigation of a Stationary Discon- 
tinuity Surface in an Electromagnetic 


Field With a Gas Conductivity Jump 


HE PROPERTIES of discontinuity surfaces in ordinary 

gasdynamics have by now been thoroughly investigated. 
The conditions on such surfaces have been derived and their 
structure and stability have been investigated. Likewise, the 
discontinuity surfaces in an infinitely conducting medium 
have been investigated in magnetogasdynamics. This paper 
treats a stationary discontinuity surface, on which, along with 
discontinuities existing in the thermodynamic parameters 
(flow and velocity), the conductivity of the gas also ex- 
periences a discontinuity. The analysis discloses certain new 
properties of the discontinuity surfaces, which have no 
analogies in ordinary gasdynamics. 


1 We consider a plane, stationary discontinuity surface, 
on which is incident a stream of nonconducting gas (o; = 0). 
The external electric and magnetic fields are considered 
parallel to the discontinuity surface and perpendicular to each 
other. We assume the parameters of the incident stream 
and the intensities of the electric and magnetic fields in front 
of the discontinuity surface to be specified. The conductivity 


Translated from Izvestiia Akademii Nauk SSSR, Otdelenie 
Tekhnicheskikh Nauk, Mekhanika i Mashinostroenie (Bull. USSR 
Acad. Sci., Div. Tech. Sci., Mechanics and Machine Building), no. 
5, 1959, pp. 9-15. 
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of the gas behind the discontinuity surface is considered in- 
finite (2 = ©). The laws of conservation of mass, momen- 
tum, energy and tangential component of the electric field 
are written in this case in the following form (5)! 


H,? 
= 2 —_ = 2 
Pit; P22 Pi + Piri +3 Pe + Pode 


2 Cc Ve? c 
Pits (2 i) E,H, = (2 + i) a ior 
E, = FE, = —(1/c)vH 


The same relations connect the parameters of the stream 
and the field on the boundaries of a certain transition region, 
which separates two moving streams, if 7; = 0, 02 = 09 ¥ ~, 
and the motion in the transition region is one-dimensional. If 
the width of this region is much less than the characteristic 
length of the problem, it can be considered as a discontinuity 
surface with a finite jump in conductivity. Upon going 
through such a discontinuity surface, the parameters of the 
stream and of the field are related by Equation [1.1]. 


1 Numbers in parentheses indicate References at end of paper. 
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If the stream and field parameters are specified in front of 
the discontinuity surface, relations [1.1] determine the values 
of all the parameters of the stream and the field behind it. 
Using the density ratio € = p:/p2, we obtain for the parameters 
behind the discontinuity surface 


ve = = —(cE,/e;) 

1 


We can determine ¢ from the cubic equation 


x? Ay — 1) 


where 
n = pi/pw;? 
= 


Here, as everywhere else, the gas is considered perfect with 
a constant ratio of specific heats, y = c,/c,. 

Eliminating the velocities from [1.1], we arrive at a rela- 
tion analogous to the Hugoniot adiabat in ordinary gasdy- 
namics 

P2 — 


(Vi + V2) (H2* — H,*) X 


(Vi + V2) — =0 [1.3] 


and when H, = H,; [1.3] goes into the ordinary Hugoniot 
adiabat. 

Let us investigate the curves obtained by plotting [1.3] in 
the p-V plane (V = 1/p). 

a Case EF, =0. Relation [1.3] assumes the following form 


— + (Vit = 0 (11.4) 


We note that unlike ordinary gasdynamics, the adiabat 
[1.4] does not pass through the point of the original state p:Vi, 
but passes below it 


(y-1) for Vi = Vs; 


This is caused by the fact that in this formulation we can- 
not have surfaces of discontinuity of infinitesimally small in- 


tensity, since by assuming an infinite conductivity jump we 
imply a finite discontinuity in the intensity of the magnetic 
field. 

The adiabat [1.4] has the following asymptotes in the p-V 
plane 


V = ——YV,=Vd = 


Let us investigate the character of curve [1.4] and its loca- 
tion relative to the Poisson adiabat 


Pp2V2" = [1.5] 


Since we have at the points of the adiabat [1.4] 


dV, 
E V2 [1.6] 
dv 2 
then 
dp» 
—— <0 f H 1 
or 1 y¥—1 [1.8] 
[1.9] 
dV, 


From this, and also from the condition p2 > 0, it follows 
that condition [1.8] gives rise to a branch of [1.4] located to 
the right of the line V = Vi(y — 1)/(y + 1), while condition 
[1.9] produces a branch of [1.4] located to the left of line V = 
Vily — 1)/(y + 1). Relation [1.7] shows that in both cases 
the curve [1.4] is concave upwards. It is easy to show that 
adiabats [1.4] and [1.5] have only one point of intersection. 

The relative placement of the Poisson adiabat and of [1.4] 
subject to condition [1.8] is shown in Fig. 1, and that subject 
to condition [1.9] is shown in Fig. 2 (1 is the Poisson adiabat; 
2 is the adiabat [1.4]). 

According to Figs. 1, the only shocks possible when condi- 
tion [1.8] and the condition of increasing entropy are satisfied 
are those for which 


p2 > p* V2 < V* 


Thus, unlike in ordinary gasdynamics, not only rarefaction 
jumps are forbidden, but also condensation jumps with small 
increase in pressure p2 — p) < p* — pr. 
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When conditions [1.9] are satisfied, the possible discon- 
tinuity surfaces are shocks with properties that are not en- 
countered in ordinary gasdynamics of nonconducting gases. 
These shocks are characterized by a larger density increase than 
the maximum possible ones in ordinary gasdynamics, (y + 1) 
p: < (Y — 1)ps, and are also characterized by the fact that 
their density rise decreases with increasing pressure jump. 

b Case Z, #0. Relation [1.3] can be converted, with the 
aid of the expression for Hz, from [1.1] to the following form 


1 
(V2 Vi) Xx 


1 2 
—— (pV: — + or 


y-1 2 


kV, 
V2 


Ve 


Henceforth we shall retain the name adiabat for the curve 
[1.10] in p-V plane even though the parameter & includes the 
velocity of the incident stream. 


pt 


x 
113 
| 


Fig. 5 


Let us write down the derivative 


dp. _ (V2) | [1.11] 


dV2 Ay — 1) 2 
where 
Y Y 
¥+1 + 1 x 
3k°V\3V_. 
+ — = + BV.? D [1.12 
+ lor + BV.? + CV2 + [1.12] 


By specifying an initial point p.,Vi, the intensity of the 
magnetic field in front of the discontinuity surface, and also 
the parameter k, we define the unique adiabat [1.10] in the 
p-V plane. 

This curve has the following asymptotes 


Hi? y¥-1 


V deg V =0 


It is seen from [1.12] that for any values of the parameters 
ahead of the discontinuity surface we have D > 0 and ¥(0) > 
0. The form of the curve [1.10], for specific values of the 
parameters in front of the discontinuity surface, depends on 
the signs of the coefficient A and of the discriminant M of the 
equation ¥(V2) = 0. Different cases are possible here, some 
of which are shown in Figs. 3-5. On the same figures are 
plotted the curves for the Poisson adiabat and those bounding 
the region of discontinuities possible from the point of view 
of conditions of increasing entropy. 

Figs. 3-5 show that along the discontinuity surfaces, which 
have the ordinary properties of ordinary gasdynamics (Fig. 3), 
discontinuity surfaces can occur with properties not en- 
countered in the absence of an electromagnetic field. Fig. 4 
indicates, for example, that discontinuities are possible in 
which the density rise decreases with increasing pressure 
jump. It is particularly important to note the possibility of 
existence of rarefaction jumps, permissible from the point of 
view of the increase in entropy (Fig. 5). 

In addition to the jumps noted, relations are also possible 
between the parameters in front of the discontinuity surface, 
under which the discontinuities become multiple-valued. To 
verify this let us plot on the p-V plane, in addition to the 
adiabats, also the curves 7 = piv; = const. These curves are 
given by the equation 


(>. Pi V2 Sr V2) {i 13] 

It is easy to check that the line V = 0 is an asymptote for 
the curves [1.13]. Furthermore, curves [1.10] and [1.13] 
have in general three points of intersection. This gives rise 
to a situation wherein the parameters behind the discon- 
tinuity surface are not determined uniquely from the parame- 
ters specified ahead of it (p1, pi, Hi, k, 7). An example of 
such an arrangement of the curves is shown in Fig. 5. In this 
case two discontinuities B and C are possible, satisfying all 
the conditions, for specified stream and field parameters. 

In conclusion let us note that if the values of pi, 71, Hi 
and & are specified in front of the discontinuity surface, the 
adiabat [1.10] becomes specified in the p-V plane. If we now 
specify one of the quantities on the discontinuity surface, for 
example, the pressure, we find the point of intersection be- 
tween the line p2 = const and the adiabat [1.10]. The curve 
j = const passing through this point determines the mass flow 
j corresponding to the discontinuity surface with a given pres- 
sure drop, and consequently the velocity and the intensity of 
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the electric field in the incident stream. If, on the 
other hand, the quantities specified ahead of the discontin- 
uity surface are pi, p:1, H; and Ej, then it is possible to plot 
on the p-V plane a family of adiabats [1.10] corresponding to 
different values of the parameter k, as well as a family of 
curves 7 = const, which depends on the two parameters k 
and 7. If we now specify the pressure jump, we find on the 
line pe = const the point of intersection of the curves [1.10 
and 1.13], which correspond to the same value of k. This 
point determines the value of 7 and all the parameters of the 
stream behind the discontinuity surface. Naturally, it is 
necessary to verify in each case whether the conditions of in- 
creasing entropy are satisfied. 

2 To determine the structure of a discontinuity surface 
with a conductivity jump, we proceed in analogy with (3, 6 
and 7). 

We consider one-dimensional, stationary motions of non- 
vicous and non-heat-conducting gas with a conductivity that 
varies along the z axis. The intensities of the electric and 
magnetic fields are considered directly perpendicular to each 
other and lying in a plane perpendicular to the z axis. In this 
case the system of equations of magnetohydrodynamics can 
be reduced to four integrals and one ordinary differential 
equation 


We shall assume that up to a certain finite point z = x the 
conductivity of the gas in zero. At x < 2% the flow is con- 
sidered moving with specified parameters, i.e., the con- 
stants ¢), C2, ¢3 and cy are determined on the stream for ¢ = 0. 
For x > 2, the flow is described by Equation [2.1] with the 
right halves specified, and the conductivity in this region 
may either be constant or a certain function of x. In order to 
determine the flow, we obtain the dependence of v on H from 
the first four equations 


1(,-2m) 
81 


81 7 + 1 4rrc, 
[2.2] 


Inserting this expression in the latter equation and integrat- 
ing, we obtain the function H(x) and then use all the final 
formulas to obtain all the parameters of motion. 

Leaving aside the exact solution of the problem, let us 
study the behavior of the curves v = v(/), specified by rela- 
tion [2.2] in the »-H plane. The expression under the square 
root in [2.2] vanishes at four (or two) values of H, which are 
roots of the equation 


(y Sr y¥+1 \a 

At these values of H the straight line, H = const, is tangent 


to the curve v = (1). The geometric locus of such points is 
a parabola 
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Since the momentum and the energy flux for the gas should 
be positive in order to realize this kind of flow, the branches of 
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\ 
(2.2) 


Fig. 8 


the curves » = v(H) for which co — H?/8r < 0 and c3 + 
cal] /4a < 0 must be considered as nonrealizable. 

The form of the curves v = v(/1) will depend on the behav- 
ior of the function 


2 2 9 ae 
(y¥ + la 87 \a 

Let us consider the case cy < 0 (for cy > 0 the situation will 
be analogous if the sign of H is reversed). The qualitative 
behavior of the function f(H) for f(0) > 0 is shown in Fig. 6, 
and the realization of positions / and 2 will depend on the re- 
lation between the constants contained in f(H). The charac- 
teristic placement of the curves v = v(/) for cases 7 and 2 is 
indicated in Figs. 7 and 8. The same figures show a plot of 
the hyperbola 


E = —(1/c)vH [2.4] 


The points of this hyperbola correspond to the points of the 
stream, where the conductivity is ¢ = © or else to those 
points, after which there is no change in the flow parameters 
with z. The latter of the equations [2.1] establishes a 
definite direction on the curve v = v(H), in which the repre- 
sentative point should move. In the case of a finite conduc- 
tivity of the stream (8) in all points — © < x < ~&, the initial 
and final points should be located on the curve [2.4]. In this 
case any discontinuity, defined by the laws of conservation 
[1.1], can be put in correspondence with a certain continuous 
flow (AB in Fig. 7). In our ease the initial point can be 
located at any point C of the upper half plane. It is easy to see 
here that continuous motions, which can be continued to in- 
finity, are possible for certain initial conditions; these corre- 
spond to the points located above [2.4], and the points in 
the segment DB lying below curve [2.4], if such points exist. 
On the other hand, if the initial parameters are such that the 
initial point falls in the region DCA on the curve v = v(H), 
then it is impossible to have continuous motion which can 
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be continued to infinity. We note that for these initial data 
the discontinuity surfaces with conductivity jumps, defined by 
relation [1.1], will become feasible and the condition of entropy 
increase can be satisfied for them. 

We note that if the initial parameters are such that the 
curves [2.3 and 2.4] do not intersect, then there always 
exists a region of initial data, corresponding to DB, from 
which stationary rarefaction flows are possible. These flows 
do not have any analogies in ordinary gasdynamics and 
correspond to the rarefaction jumps noted above. 

Thus, an examination of the structure of the discontinuity 
surface in the formulation used here, leads to the conclusion 
that not all the discontinuity surfaces, corresponding to the 
conservation laws [1.1], can be considered as the limits of con- 
tinuous flows with variable conductivity. The discontinuity 
surfaces with conductivity jumps, which can be considered 
as the limit of continuous flows, are either jumps of increasing 
density with increasing magnetic field, or else yumps of de- 
creasing density with decreasing magnetic field. 

Allowance for the viscosity and heat conduction of the gas 
shows that stationary flows are possible only under definite 
initial values of the parameters. This circumstance allows 
an electromagnetic wave, which changes the initial parame- 
ters of the field to be propagated in front of a discontinuity 
surface with a conductivity jump, contained in a stream of 
ideal gas. 

3 Discontinuity surfaces with conductivity jumps, with 
energy liberated on the front of the discontinuity, can be con- 
sidered in analogy with the foregoing. Relations [1.1] con- 
tinue to hold in this case, with the exception of the condition 
of conservation of energy, to which one adds an additional 
term g, connected with the amount of energy liberated on the 
front of the wave. 

The detonation adiabat can be represented in the form 


2 167 


— — (Vi + V2) (= 


This relation shows that the qualitative behavior of the de- 
tonation adiabat in the p-V plane is the same as of adiabat 
[1.10], except that in the region V2 > (y — 1)Vi/(y + 1) the 
adiabat [3.1] lies above the adiabat [1.10] (¢ < 0), and in the 
Pager 0< V2< (y — 1)Vi/(y + 1) it lies below the adiabat 

1.10]. 

It is easy to show that when V2 > 2(y — 1)V:/3(y + 1) 
the detonation adiabat is steeper than the adiabat [1.10], and 
that for 0 < V2 < 2(y — 1)V:/3(y + 1) it slopes more gently 


(Af, — H2)+q =0 [8.1] 


than adiabat [1.10]. 

These considerations enable us to plot in the p-V plane the 
detonation adiabats, the forms of which will be analogous to 
the curves shown in Figs. 3-5. It is also possible to plot in 
the p-V plane the curve j = const, corresponding to the sonic 
velocity of the gas behind the detonation wave 


k* 
(cmv. in V2? (Vi — Vs) [8.2] 


If the curve [8.2] intersects the adiabat [3.1] below the 
Poisson adiabat, then, after Jouguet, a detonation at these 
values of the initial parameters will be impossible. There will 
also be realized an ‘‘overcompressed”’ detonation with the gas 
moving behind the wave at a velocity less than the velocity of 
sound. In ordinary gasdynamics the realization of detona- 
tion modes under which the velocity of gas behind the wave is 
greater than the velocity of sound is impossible without energy 
being taken away at the detonation front. If the detonation 
occurs in an electromagnetic field and the conductivity of the 
medium experiences a jump behind the detonation wave, such 
modes become possible at certain values of the initial param- 
eters. The parameters of such detonation waveswere given by 
usin (8), but the cases possible here are naturally not confined 
to these examples. 

Thus, an examination of the stationary discontinuity sur- 
faces on which, in addition to the flow parameters, there is 
also a discontinuity in the conductivity of the medium, dis- 
plays certain new properties of discontinuity surfaces, which 
have no analogies in ordinary gasdynamics. 
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Reviewer’s Comment 


The previous studies of the magnetohydrodynamic jump 
conditions (see, for example (1) ), and the possible paths of 
transition from one state to another ((2), for example), have 
dealt with the case in which the end points are both in a region 
where the product of the gas conductivity and the length over 
which the flow parameters change is very large. In this case 
the “infinite conductivity” relation E = — vxBholds. The 
present paper differs from those studies in that this relation is 
taken to hold only in back of the shock, while in front the con- 
ductivity is taken to be zero so the corresponding relation is 
that E there is arbitrary. This difference leads the author to 
some interesting phenomena which are different from those 
found in the usual “infinite conductivity” case. Only the spe- 
cial case where the magnetic field is parallel to the shock front 
and the fluid velocity is perpendicular to it, is considered. 

It should be noted that not all these phenomena may occur 
in physical cases because, for example, some of them may be 
unstable. Such is the case with the “infinite conductivity” 
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shocks, where switch-on and switch-off shocks (reference 3 
below), as well as shocks which cause the tangential compo- 
nent of magnetic field to change sign (references 4 and 5 below) 
have been shown to be unstable to small disturbances. 

The reviewer would like to express his appreciation to Dr. 


Harry Petschek for his kind assistance. 
—Ne tson H. Kemp 


Avco-Everett Research Laboratory 
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the transferring body in the vicinity of the moon: 


Tdes 244,000 miles 
r 5000 miles 


= 48.8 


r 5000 miles 
ro 4000 miles 


= 1.25 


Vescres = 11.2 km/sec = 7 miles/sec 


Marth is the primary gravitational source. From Fig. 1, 
‘orres} onding to the values of \ and y 


k = 0.99 e= 0.9 


Therefore 


= € Vescres = (0.99)(0.9)(7 miles/sec) = 6.23 miles /sec 


\ssuming the initial orbit to be circular, the initial velocity is 


| lg 4 
Vr 5 x 103 


The necessary increase in velocity is thus 


AV = Vreq — 


Veirete = 1.81 miles/sec 


Successive transfers from a 200-mile Earth orbit to a 22,400- 
mile “stationary” (24-hr period) orbit: 


Ides 26,400 miles 
A= = —— = 6,29 
r 4200 miles 


r 4200 miles 
ro 4000 miles 


= 1.05 


From Fig. 1, & = 0.933 and « = 0.98. Therefore 


View = = (0.933)(0.98)(7 miles/sec) = 6.4 miles /sec 


From Equation [7] 


Vieq 6.4 miles /see 
6.29 


= 1.02 miles/see 


To remain in the 22,400-mile orbit, it is necessary to modify V 
so as to effect a circular orbit at rac. The necessary change 
in V is also given by Fig. 1, as follows (subscript 2 is used to 
designate conditions at rae.) 


r 26,400 miles 


4000 miles 


= 6.6 


Therefore 


ke = 0.707 = 0.39 


Vieqs = ke €2 Vesey = (0.707 )(0.39)(7 miles/sec) = 1.925 miles /sec 


The necessary velocity change at rae, is then 


AV = Vreqz — V = 0.905 mile/sec 


Interplanetary transfer (considering only effect of sun) 
from Earth to Mars: 


93 106 miles 
4.35 X 10° miles 


= 214 


Vescr: = 386 miles/sec 
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Transfer from a 1000-mile Earth orbit, to an orbit placing 


The sun is the primary gravitational source. Therefore 
k = 0.77 € = 0.07 
Vieqg = € Vesey = (0.77)(0.07 (386 miles/sec) = 20.8 miles/sec 


Assuming the transferring body to be initially traveling at 
Karth’s orbital speed of 18.5 miles per sec, the necessary 
change in velocity is 


AV = 20.8 — 18.5 = 2.3 miles/sec 


To effect a “soft” landing on Mars (assuming Mars is prop- 
erly oriented at takeoff to insure interception by the trans- 
ferring body) or to transfer to the Martian orbit, the velocity 
V at rae. (see Eq. [7]) must be modified to the Martian orbital 
velocity of 15.1 mile per see. The necessary change in velocity 
at raes May be computed, using Fig. 1, as was done pre- 
viously. 


Interplanetary transfer (considering only effect of sun) 
from Earth to Venus: 


40 


0.75 
10A.U. 


y = 214 
as shown previously. Therefore 


k = 0.655 € = 0.07 


€ = (0.655)(0.07 )(386 miles/sec) = 17.7 miles/sec 


Note that Vieq << 18.5 miles per see, as required for transfer to 
an interior ellipse. V as given by Equation [7] must be 
modified, as before, at rae, to insure a “soft”? landing on 
Venus. 

The last two examples are based on assumptions similar to 
those made in computing the Hohmann interplanetary trans- 
fer ellipses (i.e., only the effect of the sun is considered, and 
takeoff and arrival are made tangentially to the initial and 
terminal orbits); the results are thus identical to the Hohmann 
results. Since the curves of Fig. 1 are based on orbital trans- 
fers which depart and arrive tangentially to the initial and de- 
sired radii, all transfer requirements given by Fig. 1 are, in 
fact, minimum energy (Hohmann) transfer requirements. 


Nomenclature 


Vescrer = reference escape velocity = V 2goro, escape velocity 
at surface of primary gravitational source (Vescrer = 
11.2 km/see = 7 miles/sec for Earth, Vescn; = 618 
km/sec = 386 miles/sec for sun) 


a = semimajor axis of elliptical orbit 

Te = apogee radius of elliptical orbit 

hy = perigee radius ot elliptical orbit 

u = ro, gravitational constant associated with primary 
gravitational source, u = 9.76 X 104 miles*/sec? for 
Earth 

Jo = gravitational acceleration at surface of primary gravi- 


tational source 
= desired radius to initial radius ratio, raes/r 
= required velocity to escape velocity ratio, Vreq/Vese 
+ = initial radius to gravitational source radius ratio, r/7r» 
= escape velocity to reference escape velocity ratio, 
V ese/ Veseret 
initial radial distance from center of primary gravita- 
tional source 
desired radial distance from center of primary gravita- 
tional source after transfer is effected (perigee or 
apogee radius of transfer ellipse, depending, respec- 
tively, on whether racs < 7 Or raes > 1) 
ro = radius of primary gravitational source (7) = 4000 miles 
for Earth, ro = 4.35 X 10° miles for sun) 
Vreq = required total velocity, at initial point of transfer, to 
effect desired new orbit 
Vese = V 2u/r, escape velocity associated with initial radius r 
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Nonlinear Pressure Oscillations in a 
Combustion Field! 


GERALD ROSEN? 


Guggenheim Jet Propulsion Center, 
Princeton, N. J. 


The unsteady combustion of liquid propellant droplets 
entrained in a burned gas flow is studied in terms of an 
idealized mathematical model. We use the model to 
derive an associated equation for nonlinear pressure oscil- 
lations. Two necessary and sufficient conditions for un- 
stable oscillations follow from a discussion of this equa- 
tion. 


AJOR practical importance has been attached to the 
problem of high frequency combustion instability, in 
particular, the unstable pressure oscillations which may affect 
the performance of a rocket motor. A thorough mathematical 
treatment of the linearized aspects of this problem has been 
given by Crocco and Cheng (1-3). Basing their theory on 
the sensitive time lag postulate, these authors have derived 
certain necessary conditions for the unstable operation of a 
rocket motor. However, the work of Crocco and Cheng does 
not elucidate the nonlinear dynamics of instability. There is a 
question of whether or not an acoustical perturbation will 
grow to become a pressure oscillation of large amplitude and 
destructive capabilities. An answer to this question is not 
within the scope of linearized theory. 

Existing mathematical methods do not permit a complete 
theoretical analysis of nonlinear combustion instability. 
Nevertheless, it is possible to understand the qualitative dy- 
namics of finite amplitude pressure oscillations if the combus- 
tion field is represented by a relatively simple model. The 
present paper reports a number of theoretical conclusions 
which are derived from such a model. Of primary interest 
is the fact that a mechanism for instability follows from an 
ordinary equation for the burning rate, without evoking 
the postulate of a time lag. The stability of a pressure wave is 
essentially determined by the frequency of the oscillation and 
the pressure sensitivity of the burning rate function, in quali- 
tative agreement with the original ideas advanced by Crocco 


(4). 
Mathematical Model 


Let us consider a special model of a one-dimensional com- 
bustion field in which physical and chemical processes trans- 
form liquid propellant droplets into burned gas. If the 
average Reynolds number of the droplets is sufficiently small, 
the droplets will remain entrained in the burned gas flow. 
This we take as a basic feature of the model. Let us put 


p, gr/em = density of the mixture 
u, em/sec = mixture velocity in the z-direction 
P, gr/em-sec? = pressure 


The equations of mass and momentum conservation take the 
form 


op 
+ (pu) = [1] 
ou ou oP 
2 
p ot Or + ox [2] 


Received Oct. 5, 1959. 
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A term due to longitudinal viscous effects has been neglected 
in writing Equation [2]. 

The entropy of the completely burned gas is assumed to be 
constant. Letting w denote the weight fraction of the burned 
gas, it follows that the density of the burned gas (wp) is given 
by the isentropic relation 


wp = [3] 


in which a and ¥ are appropriate constants. 
The rate of production of burned gas is described by thi 
equation 


» +4 ~ = f(w, P; fluid element) [4 
ot ox 


The function f embodies the rate controlling physical and/o: 
chemical processes which govern the combustion of the drop 
lets. In general, f is not simple and depends parametrically o1 
the fluid element (through the local fuel/oxidizer ratio, th« 
local spectrum of the droplets, etc.). However, the burning 
rate function must be required to vanish for all fluid elements 
at the initial and final states of combustion 


f(wo, P; fluid element) = f(1, P; fluid element) = 0 [5 


where the initial weight fraction of the burned gas is denoted 
by Wo. 

Equations [1 to 4] suffice to determine the dynamical be- 
havior of the four dependent variables p, u, Pandw. Suitabl 
boundary and initial conditions are required for a unique 
solution to Equations [1 to 4]. Here we shall be interested in 
the class of solutions characterized by the boundary conditions 


lim @ = {6a} 

limw = {6b} 


Jn view of Equations [5 and 6] the region of combustion is 
confined to those finite values of « which make f greater than 
zero. Within this region the liquid droplets are consumed, and 
the weight fraction of the burned gas increases from a to 1. 


Equation for Pressure Oscillations 

By eliminating the density and the velocity (5,6), the sys- 
tem [1 to 4] can be reduced to a pair of coupled equations. 
First we evoke Equation [1] in order to define a mass coordi- 
nate which satisfies 


dy = pdx — pudt (7] 


Observe that y is a Lagrange coordinate, in the sense that a 
fixed value of w refers to a particular fluid element. Express- 
ing the dependent variables in Equations [1, 2 and 4] in terms 
of and we have 


ou 

ou oP 
oy = 0 [9] 
ow 
=f(w, P; p) [10] 


In the latter, as well as in all subsequent equations, 0/Ot re- 
fers to a time differentiation with wy held fixed; y has been 
used to represent the parametric dependence of the burning 
rate function in Equation [4]. 

The combination of Equations [8, 9 and 3] produces 


(sar) toy 7° [11] 


Although the coupled Equations [10, 11] do not involve p and 
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u, these two equations express the full content of the system 
[1 to 4]. 

Let us now work out the time differentiations in Equation 
[11]. Making use of Equation [10], we eventually obtain 


af \ oP 
— yaP't1/7 — 2f —yP —)— = 
(y¥ + 1)w foP\? of 
{| — —P [12 
yP ol [12] 


a result which displays the essential mechanism for nonlinear 
pressure oscillations in the combustion field. 


Dynamical Interpretation 


The first and second terms on the left side of Equation [12] 
letermine the characteristics. These terms propagate a 
vressure pulse with the local mass coordinate velocity 


dy/dt = + = + VW yPp 


tecalling Equation [7], we can convert this to the classical 
vropagation velocity 


dx/dt =u + WyP/p 


Thus, the inclusion of combustion does not modify the local 
phase velocity of waves. 

The third term on the left side of Equation [12] is a “fric- 
tion term” which comes from the interaction of a pressure 
oscillation with the combustion processes. If the “friction co- 
efficient” [2f — yP(Of/OP)] is positive, the amplitude of a 
pressure wave suffers exponential damping as the wave passes 
through the combustion region. On the other hand, if the 
friction coefficient is negative, the amplitude of a pressure 
wave increases in an exponential fashion as the wave traverses 
the combustion region. Hence, a condition for the unstable 
growth of a wave is 


[13] 


The right side of Equation [12] gives intrinsic “source 


terms” for the generation of pressure waves. The first term 
on the right side is positive-definite in character and stems 
from the convective nonlinearity in Equations [1 and 2]. This 
term creates and amplifies waves of a higher frequency, 
thereby steepening a pressure pulse of finite amplitude. Low 
frequency pressure waves associated with the combustion 
process are generated by the second term on the right side of 
Equation [12]. For the case of steady deflagration this term 
gives rise to a weak rarefaction wave. 

Finally, let us consider the order of magnitude of the friction 
term on the left side of Equation [12] in comparison to the 
order of magnitude of the source terms on the right side. 
The friction term is of the order f\OP/dt|, whereas the source 
terms are of the order P~'|0P/ot|? and f?P, respectively. 
Hence, the friction term will be negligible in comparison to 
one or the other source terms if (fP)~!|OP/dt! is not of the 
order unity. The condition 


—|~|~1 [14] 


defines a range of frequencies and amplitudes for which in- 
stability is possible. Notice that in the case of small ampli- 
tude acoustical perturbations, [14] requires the period of an 
unstable oscillation to be much less than the characteristic 
time required for the droplet combustion. 

The relations stated in [13 and 14] seem to be necessary and 
sufficient conditions for nonlinear combustion instability. 
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Forced Convection Heat Transfer to 
Gaseous Hydrogen at High Heat 
Flux and High Pressure in a 
Smooth, Round, Electrically 
Heated Tube' 


J. R. MeCARTHY? and H. WOLF? 
Rocketdyne Division, North American Aviation, Inc., 
Canoga Park, Calif. 


NFORMATION concerning the heat transfer characteris- 
tics of hydrogen under conditions of high heat flux is of 
prime interest in the design considerations of high perform- 
ance rocket engines and of reactors for nuclear powered rocket 
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engines. The results presented here summarize experiments 
conducted over a sufficiently wide range of pressure and heat 
flux to be of immediate interest, although the gas temperature 
at the inlet to the test section was greater than —100 F. Heat 
fluxes as high as 14.8 B/sec in.? and pressures to 1350 psia were 
achieved. The data presented herein are the initial results 
of a program directed toward the determination of the heat 
transfer characteristics of gases under conditions of very high 
surface to fluid bulk temperature ratios 7,,/T,. It is antici- 
pated that values of the ratio 7,,/T, as high as 10 can be 
achieved with the present apparatus. 

The experiments were conducted in two electrically heated 
stainless steel test sections 13-in. long, having an 1D of 0.194 
and 0.305 in. (L/D = 67 and 42.6, respectively). Iso- 
thermal starting lengths 10-in. long were provided to insure 
a fully developed velocity profile and a uniform temperature 
profile at the entrance. The gas flow rate was measured with 
calibrated sharp edged orifices, and the gas bulk temperatures 
at the inlet and exit from the test section were measured with 
four calibrated Range I, iron-constantan thermocouples in 
each of two specially designed mixing chambers. Power was 
supplied to the test section by four d-c motor-generator units 
and was computed from the measured current through the 
test section and the voltage drop along it. The inner wall 
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Fig. 1 Correlation of forced convection heat transfer results for hydrogen at different pressures 


temperature was calculated from the measured outer wall 
temperature and the energy generated in the tube wall ac- 
cording to the method given in (1).4 

The gas exhausting from the test apparatus was burned at 
a vent approximately 30 ft from ground level. During the 
course of the experiments, flares were kept burning in the test 
cell. Accordingly, several leaks and a ruptured test section 
(caused by a malfunction in flow control) resulted only in 
fires; no explosions occurred. 

The heat transfer results presented are average values for 
the region of flow having fully developed velocity and temper- 
ature profiles. Reference (2) demonstrated that the local 
heat transfer coefficient approaches to within 5 per cent of its 
fully developed value in less than eight diameters for the 
Reynolds number range to 500,000. Accordingly, the heat 
transfer coefficients employed in the Nusselt number were 
computed from Newton’s law of cooling for the center 10-in. 
length of the test section. The energy flux through the heat 
transfer surface was assumed equal to the electrical energy 
dissipated in the center 10-in. length less the energy con- 
ducted through the insulation. Heat balances for the ex- 
periments were always less than +10 per cent. 

In correlating the heat transfer results the thermal conduc- 
tivity at l-atm pressure was taken from the National Bureau 
of Standards tables (6), and the effect of pressure was ac- 
counted for according to the experimental results of Keyes in 
(7). The values for viscosity and the pressure effect upon vis- 
cosity were those given by Hilsenrath and Touloukian (8). 
The specific heat, Prandtl number and the compressibility 
data are those reported in (6). 


* Numbers in parentheses indicate References at end of paper. 
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Fig. 1 presents the forced convection heat transfer results 
for hydrogen in a smooth round tube as a function of the bulk 
Reynolds number, at different pressures with the physical 
properties evaluated at the average bulk temperature of the 
gas. A correlation based upon bulk temperature has been 
selected, since its use in design calculations is more convenient 
than correlations based upon other reference temperatures. 
The results encompassed the following range of variables: 

Bulk Reynolds number 7.8 (10%)-1.50 (10°) 


Bulk inlet temperature 360-560 R 
Average wall to average bulk 

temperature ratio 1.5-2.8 
Average wall temperatures 860-1640 R 
Average bulk temperature 559-729 R 


Inlet pressure 64-1354 psia 

The figure demonstrates that pressure, per se, has no influence 
on the convective heat transfer characteristics of the gas. 

References (3 and 4) show that the results of heating experi- 
ments having a uniform heat flux boundary condition cannot 
be correlated with the McAdams equation without including 
the surface to bulk fluid temperature ratio T,,/T,. The in- 
fluence of the ratio T,,/T, was obtained by crossplotting the 
values of Ny./Np,?-4 at constant wall temperature and Reyn- 
olds number as a function of 7,/T>. The resulting slopes 
varied from 0.29 to 0.31. An average of 0.30 was selected for 
the correlation; the equation for the line shown in Fig. 1 is 


hD VD 0.8 es —0.3 
— = 0,023 ) (1) 
ke My 


where 


h = heat transfer coefficient, B/hr ft?R 
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Fig. 2 Comparison of forced convection heat transfer results for 
hydrogen with those for helium from (5). Properties evaluated 
at the wall temperature 


D = test section ID, ft 

k, = thermal conductivity evaluated at 7',, B/hr ft R 

Pr = fluid density at the average bulk temperature, lb/ft* 
V = average bulk velocity of the fluid, fps 

My = viscosity of the fluid evaluated at 7), lb/ft sec. 
Npryp = Prandtl number evaluated at 7; 

T., T) = average inner wall and bulk temperature, R 


The correlation at Reynolds numbers less than 10° can be im- 
proved by reducing the exponent on the Reynolds number to 
0.76 and increasing the value of the constant in Equation [1] 
to 0.038. 

It is of interest to compare the present results with the heat 
transfer results for helium gas reported by Taylor and Kirch- 
gessner (5). Fig. 2 illustrates the hydrogen data (presented 
in Fig. 1) correlated with the modified Reynolds number, 
properties evaluated at the wall temperature, as recom- 
mended in (5). The solid line represents the equation 


hD 
= 0.023 (N [2] 


where the subscript w signifies evaluation at the wall tempera- 
ture. A very satisfactory correlation of the hydrogen and the 
helium data in Fig. 2 results if the constant in Equation [2] 
is increased to 0.027. 

The authors would like to express their appreciation to the 
management of Rocketdyne for permission to publish these 
results, and to the personnel at the Combustion and Heat 
Transfer Laboratory for their assistance in conducting the 
experiments. 
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Lateral Speed Indicator 
ROBERT L. SOHN? 


Space Technology Laboratories, Inc., Los Angeles, Calif. 


HIS note describes a lateral speed sensor which can 

measure the transverse, or drift speed, of an object with 
respect to the surface of a satellite or planet as the object 
approaches the body from a remote distance, e.g., the lateral 
speed of a missile approaching the moon from several thousand 
miles distance. 

For soft landings on the surface, the speed of the vehicle 
with respect to the body must be reduced by retrorockets. 
In order to cancel the total velocity, both the vertical and 
lateral components of speed must be known. Determination 
of vertical speed is accomplished relatively easily by a radar 
altimeter. But measurement of the lateral speed component 
is much more difficult, since measurements must be made 
from altitudes of over several miles, where optical devices 
are not too practical. 
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Some lateral speed measurement schemes operate on the 
principle of scanning the horizon with three or more sensors. 
To illustrate: As the vehicle approaches the surface, the 
angles subtended by the horizon sensors vary in proportion 
to both vertical and lateral speed. By noting the rates of 
change in each of the subtended angles, lateral speed can be 
computed. But the great difficulty with this system lies 
in the fact that the angular rates are extremely small, so that 
good accuracy cannot be obtained. This note describes a 
lateral speed indicator that does not depend upon rate informa- 
tion, but only upon position information. 


Analytical Background 

Consider, for example, the case in which a vehicle is ap- 
proaching the surface of the moon. Nominally, the vehicle 
would be designed to impact vertically, with no lateral speed 
component. In actuality, the vehicle probably will be off 
course by a certain distance. As a result, the vehicle will not 
approach along a radius and impact vertically, but will ap- 
proach along a slightly curved path and impact with some 
residual lateral speed. This is illustrated in Fig. 1. 

It can be shown that the lateral speed is uniquely related to 
the angle through which the trajectory curves during ap- 
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Fig. 1 Trajectory path and nomenclature 


proach. In other words, lateral speed can be predicted by 
combining certain known quantities, such as initial speed 
and altitude, with a measurement of the angle through which 
the trajectory curves during approach. For example, a 
horizon scanner can be used to establish a local vertical when 
the vehicle begins its approach to the body (position 1). A 
second measurement of the local vertical is made as the 
vehicle draws near to the surface (position 2). The differ- 
ence between these angles (the local vertical at position 1 is 
noted by a gyro) is, then, a unique indication of lateral 
speed both as to direction and magnitude. Since vehicle 
motion with respect to the gyro reference system (three- 
axis gyros) is nulled, sufficient information is known to orient 
a retrorocket so as to null the lateral speed as well as vertical 
speed. 


Analytical Development 


The motion of a body in a central force field has been 
analyzed in several texts.? Fig. 1 (from the cited reference) 
illustrates the nomenclature. It should be noted that the 
case of a hyperbolic orbit has been selected because a vehicle 
approaching the moon, for example, normally will be traveling 
faster than escape speed with respect to the moon. The 
derivation can be developed easily for the case of parabolic 
or elliptic paths. 

First of all assume that a horizon scanner will measure 
local vertical at positions 1 and 2; thus A@ will be known, 
since the angle of rotation of the local vertical is precisely 
Aé. (A gyro will be used to store the local vertical angle at 
position 1.) Also the altitudes at positions 1 and 2 will be 
measured by a C-W range and range rate radar altimeter (which 
is directed along the local vertical); thus 7; and re will be 
known. Further, measurement of range rate along the local 
vertical at position 1 will give V; to a very close approxi- 
mation. Actually the vehicle velocity vector at position 1 
will be rotated very slightly from the local vertical (by less 
than 2 deg, as will be shown later). The true velocity will 
then vary from the measured velocity by the cosine of the 
rotation angle. But since this angle is usually small, the 

2A summary of the pertinent relations is given in “A Pro- 


posed Kepler Diagram,” by R. L. Sohn, ARS Journat, vol. 29, 
no. 1, Jan. 1959, pp. 51-54. 
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measured velocity will differ only slightly from the true 


velocity. 
The pertinent equations are as follows 
= V%r/GM (1) 
r Tp 
A-2 e-1 
(e + 1)r,/r 
-1 3 
1; = (1 = 


= cos~! E — (1+ e) [4 
e r 


At position 1, V and r are measured, so that \; can be com- 
puted 


= 


As noted, A@ is measured by the lateral speed sensory instru- 
mentation. From Equation [4] 


A6é = cos“! E — (1+ e) = 
1 —(1+e)[- [5} 
e re 


Equation [2] can be rearranged to solve for rp 
rp = n(e — 1)/(Ar — 2) (6! 


Combining Equations [5 and 6] gives 


All quantities in Equation [7] are measured, except the 
eccentricity e. Hence a solution can be obtained for e. 
Actually, the computer incorporated into the lateral speed 
indicator will solve the equation by an iterative process. 

Finally, it is necessary to solve for the flight path angle 
and velocity at the surface of the body. The radius of the 
body R will be known beforehand, of course. The solution 
for the velocity is as follows. From Equation [2] 


As = (Ar — 2) (R/n) + 2 [8] 


and from Equation [1] 
V2 (9) 


Combining Equations [8 and 9] gives 


24+ (A 2)R/n (2) 2 
J R V; {10} 


and by combining Equations [3 and 6] 


Ke? — 1)/(A1 — 2))]n/R 
2+ (Ai — 2)/(n/R) 


A direct solution for the lateral speed at the surface of the 
body can now be given as 


Vu = Vo sin By = R 2) [12} 


B> = sin-! {11} 


Numerical Example 


A vehicle approaches the moon for a controlled landing. 
It is assumed to be off course by 100 nautical miles, with a 
resultant lateral speed component of 435 fps, which is com- 
puted by the sensor as previously described. 
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Description 


A vehicle is approaching the moon along a nominally radial 
path. However at a distance of 25,000 nautical miles from 
the center of the moon, the vehicle is proceeding parallel to 
its nominal path, but 100 nautical miles to one side. Speed 
at this point is 3930 fps. It is desired to predict impact 
speed as described. 

At position 1. Assume that position 1 is 1000 nautical 
miles above the surface of the moon. The lateral speed 
indicator measures the following: Velocity, 6480 fps (radar 
range rate); altitude, 1000 nautical miles (radar range); local 
vertical, by means of horizon scanner. 

At position 2. Assume that position 2 is 100 nautical 
miles above the surface of the moon. The lateral speed 
indicator measures the following: Velocity, 8380 fps; alti- 
tude, 100 nautical miles; A@, 1.38 deg (by comparing horizon 
scanner reading with position 1 horizon scanner reading). 


Computation 


From Equation [7], eccentricity e is 1.00123. 

From Equation [10], V, is 8570 fps. 

From Equation [11], 8, is 2.92 deg. 

From Equation [12], V,, is 435 fps. (Direction of motion 
is given by comparing the horizon scanner reading with the 
three-axis gyro reference.) 


Error Analysis 

The sensitivity of the lateral speed indicator to errors in the 
measured quantities can be determined as follows. 
Equation [7] can be rewritten as shown by inserting Equation 
[12] 


1 R* 


By assuming that 
62 


es 1.0 Al — — [14] 


the following relations are obtained 


2k? 

Nou 

va 

Hence 
GM 

- = 15 
Ve — [15} 


The sensitivity equations are then obtained with suitable 
differentiation 


or; 2h? (avr, 1/V 2r,3/2 
avy _ 


1 
ore QR? (1/V — 1/Vr2)? 18] 


By using the numerical example presented previously the 
following sensitivities are indicated 


h=n—R 

%OVy = 0.847%6( AA) 
= 0.598% 
= —0.152% (he) 


where the terms on the left are percentage changes in V;, with 
percentage changes in measured inputs. The percentage 
error in A@ is likely to be larger than those of the altimeter. 
For example, the time of flight between points 1 and 2 is 
approximately 0.2 hr. If a gyro drift of 1 deg per hr is 
assumed, a 0.2-deg error in A@ will be introduced. The cor- 
responding error in Vy, is 53 fps. 

It should be noted that the assumptions of Equation [14] 
are valid for the example chosen, but cannot be used in the 
general case. 


New Method for Studying 


Polymerization of Solid Propellants 
and Propellant Binders 


R. W. WARFIELD! 
U. S. Naval Ordnance Laboratory, Silver Spring, Md. 


ODAY’S great usage of solid propellants has emphasized 

the need for new experimental techniques to study the 
polymerization of these composite systems. In general, the 
conventional methods are not applicable because of the poly- 
functional nature of the propellant binder, and because of the 
presence of large amounts of inorganic oxidizer particles and 
other additives. The author has recently shown (1, 2)? that 
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electrical volume resistivity measurements, made contin- 
uously and isothermally, can be used to quantitatively define 
the bulk polymerization characteristics of thermosetting poly- 
mers. This technique has now been applied to the study of 
the polymerization of solid propellants and propellant bind- 
ers. It is the purpose of this note to indicate how resistivity 
measurements may be used to study the polymerization of 
solid propellants and to describe briefly the types of data 
which may be obtained. 

The continuous current monitoring device (CCMD), which 
is a simple application of Ohm’s Law, is used to determine the 
resistivity of the propellant during and after polymerization 
(3). This device consists of a Keithley Model 210 electrom- 
eter, a Keithley Model 2008 decade shunt and a sample cell. 
With this instrumentation it is possible to monitor the large 
change in resistivity (10° to 10'* ohm-cm) which propellants 
may undergo during polymerization. The sample cell, in 
which polymerization is conducted, consists of two small con- 
centric cylinders between which about 13 gm of the liquid 
propellant or binder are poured. A small laboratory oven 
was used to maintain the isothermal conditions within +0.5 C. 
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The temperature dependence of the resistivity of the poly- 
merized propellant was determined by monitoring the re- 
sistivity as the sample was heated at the rate of about 15 C 
per hr. The sample cell and experimental arrangement em- 
ployed to continuously record the electrical volume resistivity 
have been previously reported (3). 

From a series of these polymerizations conducted on both 
the binder and the composite propellant, it is possible to deter- 
mine the rates of polymerization at various temperatures and 
the overall activation energies for the polymerization process 
of both the binder and the composite propellant. From the 
values of the activation energies it can be determined if the 
presence of oxidizer and/or additives serves to catalyze and/or 
accelerate the polymerization. It is expected that a catalyzed 
reaction would exhibit a lower activation energy than an un- 
catalyzed reaction. A rapid and simple method for determin- 
ing catalysis in these systems is particularly significant, in view 
of the work of Mishuck and Carleton (4) who have reported 
that the presence of ammonium perchlorate will accelerate the 
polymerization of a polyester-styrene propellant. Results ob- 
tained in this laboratory (5) by the electrical resistivity method 
indicate that the presence of ammonium perchlorate will cata- 
lyze the polymerization of a polyurethane propellant. 

By studying the rates of polymerization at different tem- 
peratures, it is possible to select the optimum manner in which 
to conduct polymerization and thus obtain the optimum phys- 
ical properties. Also, providing that the polymerizations are 
conducted over a broad temperature range, Arrhenius plots 
may be used to determine if different reaction mechanisms oc- 
cur. It has been previously suggested (6) that the method 
employing electrical resistivity during an isothermal poly- 
merization may give results superior to those obtained by con- 
ventional dilatometric techniques. 


Technical Comments 


The temperature dependence of the resistivity of the poly- 
merized propellant can be used as an index of the extent of 
polymerization. It has been shown (1, 2) that during poly- 
merization the temperature dependence of the resistivity in- 
creases with the extent of polymerization. Incomplete poly- 
merization is evidenced by a resistivity vs. temperature plot 
exhibiting a reduced slope and a displacement toward lower 
temperatures when compared with a plot of the completely 
polymerized material. 

It is suggested that electrical volume resistivity measure- 
ments may be used for long-term surveillance testing of com- 
posite propellants. On the basis of past work on thermoset- 
ting polymers, it would be expected that volume resistivity 
measurements would be sensitive to degradation, decomposi- 
tion, crystallation, and possibly could be used to study deterio- 
ration at the binder-oxidizer interface (7) and other effects 
encountered in long-term aging. 
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Comment on “Effect of Thrust 
Termination Process Upon Range 
Dispersion of a Ballistic Missile” 


STEPHEN L. JOHNSTON! 


Army Rocket and Guided Missile Agency, 
Redstone Arsenal, Ala. 


HE ELLIPTICAL range equation has appeared in many 
forms. The recent article by Kelly (1)? presents the 
following range Equation [14] in (1) 


Vb? + 4ac + 4ac 


Although Kelly indicated that this equation was taken from 
the 1953 edition of (2), he did not show the derivation by 
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which he obtained this form. Equation [14] does not appear 
in the 1957 edition of (2). 

The range equation given by Baker and Hart (3) indicates 
that Kelly’s Equation [14] may be incorrect. Baker and 
Hart’s equation in Kelly’s notation is 


cos" ( b— cos" ( b — WZp ) 
P + 4ac a/b? + 4ac 


The quantity c must be redefined’ as [(ro + ho)Uo cos Oo]. 
Other quantities are as defined by Kelly. 

No attempt was made to verify Equations [16 and 17] of 
(1), but their values for the sample problem were compared 
with results obtained by using the equations of (4). These 
equations were obtained by differentiating Baker and Hart’s 
range Equation [2]. This comparison is shown in Table 1. 
The small value of 0R/06 (57.3 nautical mile/radian) quoted 


® Alternatively, quantity c could remain as defined by Kelly 
provided 6 were redefined as the angle with the vertical rather 
than the horizontal. Inspection of Table 1 indicates that Kelly 
may have inadvertently done this in his calculations. For ex- 
ample, compare his value of 480 naut mile for R; for 0 = 30 deg 
with the value of 482 naut mile for 6) = 60 deg in Table 1. 
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by Kelly indicates that 6) ~ Opt. In fact, for the ho and Up 
values shown, Oop: = 40.954 deg. It was not possible to obtain 
Kelly’s values of R, OR/OUo, and OR/O6 for any one value 
of 6. The results shown in Table 1 for (4) were obtained on 
an IBM 704 computer for various values of 6 and were 
rounded off for use herein. Values of g and ro are the same 
as those used by Kelly. 

Comparison of Kelly’s values [method of (1)] indicates 
that there may be errors in his Equations [16 and 17] in 
addition to his errors in Equation [14] and his confusion over 
A. 
While a generalized approach as was used by Kelly is very 
desirable, two possible simplifications for many practical ap- 
plications should be mentioned. The angle & used by Kelly 
is actually the angle of attack. Both load factor considera- 
tions and control system linearity problems demand that & be 
small, especially in the vicinity of thrust termination. For 
the preliminary analysis purposes of Kelly’s paper, it would be 
permissible to use £ ~ 0 and thus to simplify Equation [11] to 
= OR/OU>. 

A further justification for this simplification of Equation 

11] arises from maximum range considerations. It is well 
known that maximum free flight range is obtained for that 
value of 6) for which 0R/06) = 0. In the initial design con- 
sideration of a missile, performance aspects at maximum 
range are of principal interest. For this case 0R/06) = 0, and 
again simplification of Equation [11] and the elimination of 
;quation [17] are permitted. 

For most cases (i.e., 0R/O0%) = 0 and/or & = 0), Equation 

16] can also be eliminated. From Fig. 4 of (4), it can be seen 
that for optimum trajectories OR/OU> is very insensitive to 
changes in values of ho for a gross evaluation. Similarly, 
OR/OU> is approximately a linear function of free-flight range. 


Table 1 Comparison of results obtained by methods of 
references (4) and (1) 


90, R,, oR/ol oR ‘OOo, 

deg naut miles naut miles/fps naut miles/rad 
| Method of (4)} 

30 554.3 0.1138 400.0 

40 592.8 0.1218 35.57 

50 567 .0 0.1150 —326.2 

60 481.9 0.0961 — 637.2 
| Method of (1)! 

30 480 0.117 57.3 


Thus Equation [11] could be further simplified to the form 
Q = K, + AR, where 


K, = 0.0335 naut mile/fps K, = 1.33 X 10~‘ sec/ft 
These values of A, and Ke are valid for missiles having ranges 
of approximately 200 to 2500 nautical miles. 
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Aeromechanieal Effects in Ablation 


MARK V. MORKOVIN! 


The Martin Co., Baltimore, Md. 


2 


N HIS article on ‘Recent Advances in Ablation’ (1)? 
Adams states that “Mechanical erosion, though not dis- 
cussed here, will reduce the effective heat of ablation. Since 
surface erosion effects cannot be separated from ablation as 
such, the shear and pressure forces anticipated for application 
should be duplicated in laboratory ablation tests.”” At the 
discussion on ablation at the Conference on Physical Chemis- 
try and Space Flight, held Sept. 1-3, 1959 at the University of 
Pennsylvania in Philadelphia, Barry and Gruntfest of the 
General Electric Co. emphasized that good mechanical 
strength of the char is an important property of good ablators. 
The following is the essence of the comments at the same con- 
ference by the author on the question of laboratory simula- 
tion of the mechanical erosion aspects of ablation. 

Let us consider another possible interplay between fluid 
mechanics and the rate of ablation. The importance of 
mechanical strength of the char suggests that the dynamic 
aspects of the local air loads during the testing and during 
flight should be considered. 

From information presented at the International Sym- 
posium on Fluid Mechanics in the Ionosphere (2), it appears 


Received Jan. 15, 1960. j 
1 Principal Scientist, Research Staff, Manned Vehicles Division. 
? Numbers in parentheses indicate References at end of paper. 


Aprit 1960 


highly unlikely that there are any irregularities in that part of 
Karth’s atmosphere where ablation may be important which 
has a scale of lessthan 5m. Thus for all practical purposes the 
ablating vehicle encounters a uniform and steady (or quasi- 
steady) free stream and the aerodynamic disturbances present 
must arise in the boundary layer of the vehicle itself (or due 
to the vibrations of the body). 

The situation is quite different in the aerodynamic testing 
facilities where one should distinguish between steady and 
unsteady temperature, sound, and turbulence nonuniformi- 
ties of the free stream (3). For instance, steady nonuni- 
formity in velocity has been observed to cause considerable 
unevenness in ablation rate in the Caltech Hypersonic Tunnel 
which is otherwise a clean and quiet facility (private com- 
munication by Lees’). In hot jets with diameters less than 2 
diam of the tested body the steady nonuniformities are also 
likely to be important. However, in hot jets, plasma jets and 
rocket exhausts, the unsteady temperature, sound, and tur- 
bulence nonuniformities may well be most important and 
often critical. While nothing definite is known about the 
effect of rather intense moving hot-spots on the body, we 
know that the extremely high sound levels encountered in 
rocket exhausts can be highly destructive of structural com- 
ponents. This may occur in a number of ways depending 
upon the frequency of the sound, since in one way or another 
the critical effects are usually associated with resonance of 
local structural elements or of aerodynamic cavities. Strong 
sound fluctuations may create local spots hotter than stagna- 
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At x 


tion temperature and local instantaneous high pressure regions, 
especially inside cavities, which must be detrimental to the 
structure, see e.g., (4) on “frontal’’ cavities, and unpublished 
data on “side” cavities. (Even the mild pressure and shear 
fluctuations associated with a low speed turbulent boundary 
layer are known to accelerate substantially the mechanical 
scrubbing of paint from ship hulls.) Finally, besides the 
above direct effects of the nonuniformities in the testing 
facilities, there must be an indirect effect, namely the earliest 
possible transition to turbulent flow over the model itself with 
consequent higher heat transfer and still higher pressure 
fluctuations. 

Clearly, all these effects depend upon the amplitudes of the 
disturbances which have not been measured in the hot facili- 
ties. Whatever these amplitudes may be, the gross informa- 
tion indicates that they are one or two orders of magnitude 
larger than those we are used to in ordinary aerodynamic 


tests and those the ablating vehicle is likely to encounter in 
actual flight. The reported importance of the mechanical 
strength of the char seems to indicate that the “‘aerodynamic 
abrasion” in the “noisy’’ testing facilities may have to be con- 
sidered in evaluating the significance of the tests and that 
actual flight vehicles may ablate slower. 
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Book Reviews 


Analysis and Design of Aircraft Structures, 
Vol. I, Analysis for Stress and 
Strain, by E. F. Bruhn, assisted by A. F. 
Schmitt, Tri-State Offset Co., Cincin- 
nati, Ohio, 1958. 

Reviewed by GEORGE GERARD 
New York University 


In approaching a subject as broad as 
aircraft structures, there are certain 
courses available to an author: He may 
present a large volume of material that can 
cover most of the problem areas likely to 
be encountered in one’s _ professional 
career; or, he may proceed in a highly 
selective manner by choosing broad areas 
of importance where the student is 
thoroughly grounded in the scientific 
principles of the subjects selected. With 
this fundamental background, one can 
proceed to handle the specific problems 
which arise in his professional work, as 
well as treat the new developments 
characteristic of this rapidly evolving 
field. 

A large measure of the success of the 
earlier editions of this book, as measured 
by its routine use in industry, can probably 
be attributed to the use of the first 
approach. The current edition consists of 
two volumes, the one under review and 
Volume II, “Strength Analysis and 
Design” yet to appear. Thus, there is a 
distinct effort to cover much more 
material in detail in the present revision. 

To this reviewer, the needs of the 
modern academic world are more con- 
sistent with the second approach discussed, 
since progressive engineering education 
tends to emphasize the science rather than 
the art of the profession. The great 
amount of detailed material covered in 
“Analysis and Design of Aircraft Struc- 
tures’”’ will require each instructor to 
make his individual selection of material. 
The fierce competition for the student’s 
time in his education as an aeronautical- 
astronautical engineer necessitates con- 
siderable judgment in selection of subject 
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material for aircraft structures courses. 
Thus a major role that should be per- 
formed by the author of a textbook has 
been left to the instructor. - 

The material covered in Volume I is 
concerned with the stress analysis of 
statically determinate and indetermin- 
ate structures, axial loads and shear flow 
in thin wall shells, and contains an 
introduction to wing and fuselage analysis. 
The excellent photographs of modern 
aircraft structural components throughout 
the book are commendable. 

To the stress analyst in industry, the 
appearance of this new edition may be of 
considerable interest particularly because 
of the enlarged coverage. To those 
concerned with the educational aspects of 
aeronautical structures, the enlarged scope 
of the work may prove these volumes to 
be unattractive as a text. 


Hypersonic Aerodynamics, by Robert 
Wesley Truitt, The Ronald Press Co., 
New York, 1959, 474 pp. $10. 
Reviewed by H. Dorrance 

Convair Scientific Research Laboratory 
As stated in the Preface, “‘the purpose of 

this book is to present the fundamental 
principles of hypersonic aerodynamics and 
their application to missile design prob- 
lems.”’ In the opinion of this reviewer the 
author comes closer to the latter than to 
the former, as will be elaborated upon in 
what follows. 

This book has some merit. The 
stated objective is excellent, the organiza- 
tion by chapter headings good and the 
references at the end of each chapter 
excellent. Certainly the inclusion of 
problems with their answers, according to 
methods presented in the book, is a meri- 
torious accomplishment for teachers and 
those who learn by example. The subject 
matter suggested by the table of contents 
is excellent. It is organized into 10 
chapters which, starting with an intro- 
duction to hypersonics, ranges through 


methods applicable in analyzing the two- 
and three-dimensional flow about bodies 
the development of minimum drag bodies 
using impact (Newtonian) theory and 
impact theory with centrifugal force 
corrections, hypersonic small disturbance 
theory, hypersonic slender body theory, 
hypersonic flow over blunt bodies, com- 
plete configurations, stability and control, 
and ends with a chapter on real gas effects. 

This reviewer found Chapter 3 useful. 
It is on impact theory (Newtonian), 
modified to account for centrifugal force 
effects for bodies of revolution, and em- 
pirical methods of estimating cone-cylin- 
der normal force and pitching moment. 
Similarly, Chapter 4 on the development 
of minimum drag bodies using impact 
(Newtonian) and impact plus centrifugal 
force correction theories, Chapter 6 on 
Cole’s slender body theory, and Chapter 9 
on stability and control, are well organized 
and, probably useful, chapters. In most 
cases the author’s efforts to compare 
theory with experiment were clearly and 
pertinently presented. 

Despite the virtues cited, this reviewer is 
of the opinion that this is not a good book. 
The virtues are primarily those of the 
author’s source material; the short- 
comings, for the most part, are those 
contributed by the author. For example, 
in pages 34-36 the author sets out to 
derive the potential flow equation for 
compressible flow (although he doesn’t 
identify it that way but instead calls it the 
“Governing Differential Equation for 
Hypersonic Flow,’ certainly an unneces- 
sarily restrictive description). He uses 
what at first glance appears to be vector 
analysis but which makes no acknowl- 
edgment as to the formalities of the use of 
scalar and vector products. He persists in 
calling the Euler equation the ‘‘Navier- 
Stokes” equation. Certainly a book 
dedicated to elaborating the fundamental 
principles should avoid this error. 

The author consistently ignores the 
contribution of axial force distribution to 
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pitching moment and to the location of the 
center of pressure, both for axially 
symmetric and nonsymmetric _ bodies. 
For example, on page 84 he gives the 
center of pressure location for a cone as 
two thirds the length, whereas, because 
of the axial force distribution, it is actually 
2/;[1 + tan%é,] where 6, is the cone semi- 
vertex angle. For a cone with 6, = 25 
deg (not unreasonable for ICBM-type 
nose cones) the author’s expression is in 
error by about 20 percent. For practicing 
engineers such an omission can be damag- 
ing. 

The book varies in sophistication from 
chapter to chapter. Chapter 8 on com- 
plete configurations is so elementary as to 
suggest (however improperly) that the 
more difficult material presented earlier 
(such as hypersonic small disturbance 
theory, or hypersonic flow over blunt 
bodies) is superfluous for a_ practical 
engineer. This chapter detracts from the 
material presented fairly skillfully in 
Chapters 3 and 4. 

The important areas of hypersonic 
boundary layer theory and nonequilibrium 
hypersonic flow are merely mentioned, and 
their treatment is inadequate. 

The use of three references for Chapter 
8, which were papers “presented at the 
Engineering Section of the Virginia 
Academy of Sciences,”’ certainly raises a 
question as to their availability to the 
general reader. 


Analysis of Linear Systems, by David K. 
Cheng, Addison-Wesley Publishing Co., 
Inc., Reading, Mass., 1959, xiii + 
431 pp. 

Reviewed by J. L. Stewart 
University of Southern California 


For many years, one classic work has 
dominated the market for senior-graduate 
books for electrical engineers and others, 
in the general area of linear systems 
analysis by means of the Laplace trans- 
form method. 

The book reviewed here is typical of 
recent books, although it is above average 
in the care of its presentation and in its 
thoroughness. 

The content of the text is typical in 
that the principal purpose is the exposition 
of the Laplace transform method. The 
author is aware of the shortcomings of 
presenting the Laplace method as a mere 
definition in the form of a suddenly 
presented formula. So he has wisely 
prefaced his book (Chapters 1 and 2) 
with material on classical differential 
equations. 

The author also believes that adequate 
exposition of the Laplace method demands 
knowledge of the equations of specific 
systems. He thus gives a fairly general 
presentation of electric circuit analysis in 
Chapter 3, and of mechanical and electro- 
mechanical circuits in Chapter 4; these 
subjects come well before the introduction 
of transform methods. 

The familiar procedure is followed in 
Chapter 5; the Fourier series and the 
Fourier integral are introduced as a 
preliminary step to a derivation of the 
Laplace transform. The author does pause 
however, to discuss the invaluable concept 
of the frequency spectrum. 
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Chapters 6, 7 and 8 contain the principal 
material on the Laplace method. The 
presentation is fairly thorough (consider- 
ing the fact that contour integration is 
presumed to be beyond the level of the 
book), although to some it may be dis- 
appointingly lacking in novelty. 

Later chapters concern themselves 
with the usual supplementary material, 
such as feedback and block diagrams 
(Chapter 9), the Z transformation (Chapter 
10) and systems with distributed con- 
stants (Chapter 11). An appendix on 
numerical factoring methods is also 
provided. 

The book is a careful and highly con- 
servative presentation of a widely taught 
group of subjects; it is in fact superior to 
most books of its type. However, in its 
conservatism it has the same oversights 
and is subject to the same limitations as 
its predecessors. For example, the pole- 
zero interpretation is given almost no 
attention. The author still shows in- 
hibitions to the effect that the rudiments of 
contour integration are somehow too 
mysterious to explain. Numerous addi- 
tional items of this type can be cited. 

This reviewer believes the book to be 
adapted to the needs of a number of 
institutions; it is a good example of its 
type. However, the reason for its exist- 
ence apart from that of furnishing competi- 
tion can be questioned; it adds little to the 
literature. 


Space Flight, by Carsbie C. Adams, 
McGraw Hill-Book Co., Inc., New York, 
1959, xvi + 373 pp. $6.50. 

Reviewed by F. C. Duranrt IIT 
Avco—Research and Development 


An increasing flood of works on space 
flight have been appearing on the book 
market in the past few years. No 
longer can the ‘“aficianado’’ add to his 
library everything published on astro- 
nautics. Selection is difficult. Many books 
are narrow in scope, highly technical 
state of the art surveys, aimed at a 
limited audience. Others, intended to be 
educational on a general level, are com- 
pilations of summaries of other writings. 
Often the authors of the latter type are 
professional writers who are new to their 
subject, but are educating themselves as 
they assemble facts and figures in a 
mechanically skilled manner. 

Seldom does a general work appear 
which is not only broad in scope and 
selectively documented, but imaginative 
and stimulating as well, as_ history 
evolves into present-day projects and the 
paths into space become defined. Dr. 
Adams with his collaborators, F. I. 


Ordway III, Heyward E. Canney Jr. and ° 


Ronald C. Wakeford, have produced a 
remarkable book of wide readership 
value. Fordifferentreasons “Space Flight” 
should satisfy the needs of public and 
technical libraries, college undergraduates, 
the motivated high school senior and the 
professional rocket engineer and scientist. 
The layman needs education in the inter- 
play of widely differing disciplines of 
science and technology which are the 
components of astronautics. The pro- 
fessional often desperately needs a view of 


the ‘‘forest’’ while he is daily applying his 
calipers to the individual “‘tree.”’ 

In his Forward, Dr. Wernher von 
Braun writes that he has ‘‘no doubt that 
‘Space Flight’ will become the spaceman’s 
book on space flight,’ and that he is 
“certain that it will soon attain the 
stature on one of the few great classics on 
this fascinating and many-faceted subject.” 
This is high praise by a leader in the field. 

The technique of presentation follows an 
excellent organizational pattern in which 
historical review in each case leads the 
reader to the present. Dr. Adams uses 
skill and good judgment in his treatment 
of the recent past with regard to the 
development of satellite plans, programs, 
controversies and the personalities con- 
cerned. In particular, the positions of the 
USSR and the United States is finely 
drawn with respect to capabilities, drives 
and accomplishments. Much of this 
recent history has not appeared before as 
completely or as well. 

Particularly refreshing are the footnote 
references and extensive bibliographies 
listed at the end of all but one of the 15 
chapters. A perusal of these selected 
references is like opening a door to a 
whole new library of material. Many of 
these are the classic, technical papers 
written for meetings of the American 
Rocket Society, British Interplanetary 
Society and the International Astro- 
nautical Federation. Librarians of uni- 
versities and industrial research organiza- 
tions can usually obtain these papers 
without much effort. British, French, 
German and Italian works are also 
referenced. The writings of persons who 
have played key roles in the identification 
of astronautical problems and approaches 
to their solutions are listed, including not 
only the well-known names of von Braun, 
Van Allen, Singer, Whipple, Ehricke, 
Stehling, Stuhlinger, Bonney and Merrill, 
but foreign writers, such as Clarke, 
Cleaver, Moore, Saenger, Gartmann, Par- 
tel, Burgess, Lawden, Gatland, Carter, 
Ananoff, and dozens more. 

The writing style is light in touch—a 
difficult and seldom achieved objective. 
There is humor and whimsy, as ap- 
propriate, in treatment of history and 
philosophy. Fundamental parameters of 
rocketry and orbital and celestial mechan- 
ics are developed with refreshingly little 
mathematics. 

No critical review could fail to find 
faults in every book. The treatment of 
atmospheric re-entry could and should be 
stronger. The section on electric pro- 
pulsion seems too short, but the reviewer 
recognizes his own personal bias in this 
regard. Offsetting these shortcomings is 
the remarkably fine discussion of human 
factors and space medicine of which the 
author is professionally well aware. A few 
typographical errors, such as misspellings 
of: “fuze,” and ‘Allegany’ Ballistics 
Laboratory, and at least one incorrect 
middle initial in names. But this is 
carping, and a new edition will undoubt- 
edly make these corrections. 

In short, “Space Flight’ is a valuable 
and exciting book—informative, thought 
provoking and readable. It is definitely 
recommended for the layman and the 
professional. 
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M. H. Smith, Associate Editor 
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Propulsion and Power 
(Combustion Systems) 


An Investigation of Ejectors Without 
Induced Flow—Phase I, by D. L. Barton 
and D. Taylor, Arnold Engng: Dev. 
Center, AEDC TN 59-145, Dec. 1959, 42 
pp. 

Dimensioning Tables for Ramjet En™ 
gines. Part 2. One-dimensional inter~ 
nal Flow With Chemical Heating at 
Flight Mach Numbers | to 6 and Altitudes 
from 0 to 42 Km, by Irene Singer-Bredt 
and Hermann Stiimke, Forschungsinst. 
fiir Physik der Strahlantriebe, Stuttgart, 
Mitteilungen no. 18, Aug. 1959, 829 pp. 
(In German.) 


Propulsion and Power 
(Non-Combustion) 


Experimental Investigations of Plasma 
Accelerators for Space Vehicle Guidance 
and Propulsion, by B. Gorowitz, K. Moses 
and P. Gloersen, General Electric Co., 
Missile & Space Vehicle Dept., TIS 
R598D466, Nov. 1959, 25 pp. 

Magnetohydrodynamic Channel Flow 
of a Perfect Gas for the Generation of 
Electrical Power, by G. W. Sutton, General 
Electric Co., Missile & Space Vehicle Dept., 
TIS R598D472, Dec. 1959, 29 pp. 

Ion Wave Instabilities, by Ira Bernstein, 
Edward Frieman, Russel Kulsrud and 
Marshall N. Rosenbluth, Princeton Univ., 
Proj. Matterhorn, Matt-23, Nov. 1959, 
7 pp. 

Stability of Large Amplitude Waves in 
the One-dimensional Plasma, by David 
Montgomery, Princeton Univ., Proj. Mat- 
terhorn, Matt-24, Nov. 1959, 9 pp. 

Electromagnetic Radiation from an 
Ionized Hydrogen Plasma, by S. M. Ber- 
man, Space Tech. Labs., Inc., Physical 
Res. Lab., PRL 9-29, Sept. 1959, 25 pp. 

Space Charge Limited Current Flow 
from a Plane, Plasma ‘‘Cathode,’’ by 
B. D. Fried and L. O. Heflinger, Space 
Tech. Labs., Inc., Physical Res. Lab., Tech. 
Rep. 59-0000-00870, Oct. 1959, 12 pp. 

Akademiia Nauk, SSSR, Institute 
Atomnoi Energii, Physics of Plasmas and 
Problems of Controlled Thermonuclear 
Reactions, Moscow, Akademiia Nauk, 
1958, 4 vol. (in Russian). Vol. 1 (300 pp.): 

Theory of the Magnetic Thermonuclear 

Reactor, by I. E. Tamm, pp. 3-19 

(Chap. 1) (1951). 

Theory of the Magnetic Thermonuclear 

Reactor, by A. D. Sakharov, pp. 20-30 

(Chap. 2) (1951). 

Theory of the Magnetic Thermonuclear 

Reactor, by I. E. Tamm, pp. 31-41 

(Chap. 3) (1951). 

Questions Associated with the Drift of 

Particles in a Toroidal Magnetic Ther- 

monuclear Reaction Chamber, by G. I. 

Budker, pp. 66-67 (1951). 

Stationary Conditions of a Magnetic 


Ep1tor’s Note: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged. 
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Thermonuclear Reactor, by D. N. 
Zubarev and N. V. Klimov, pp. 249-288 
(1952). 

Vol. 3 (363 pp.): 
Thermonuclear Reactions in a System 
with Magnetic Mirrors, by G. I. Budker, 
pp. 4-31 (1954). 


Propellants and Combustion 


Electrical Initiation of Insensitive Ex- 
plosives, by D. B. Moore and G. M. 
Muller, Stanford Res. Inst., Poulter Labs., 
Tech. Rep. 016-59, Nov. 1959, 36 pp. 

Some Safety and Handling Characteris- 
tics of a Mixture of 70°, Propyne and 30% 
Allene, by Lloyd E. Line Jr., Experiment 
Inc., T M-1104, July 1959, 17 pp. 

Experimental Investigation of Flame 
Stabilization in a Deflected Jet, by Gunnar 
Erik Broman, Calif. Inst. Tech., Jet Prop. 
Lab., Rep. 30-5, Sept. 1959, 60 pp. 

Mathematical Approach to Solid Pro- 
pellant Grain Design, by H. K. Lumpkin, 
Army Rocket & Guided Missile Agency, 
Ord. Missile Labs. Div.. ARGMA TN 
1G5N, Dee. 21, 1959, 30 pp. 

Special Report; Detonation Charac- 
teristics of PBAA Propellant, by R. C. 
McCauley, Thiokol Chem. Corp., Redstone, 
Rep. 43-59, 15 pp. 

Absorption and Dispersion of Micro- 
waves in Flames, by Jiirgen Schneider 
and F. W. Hofman, Phys. Rev., vol. 116, 
Oct. 15, 1959, pp. 244-249. 

The Effect of an Electromagnetic Field 
on the Detonation Regime, by G. 
Lyubimov, Soviet Physics-Doklady, vol. 4 
(Russian vol. 126, nos. 1-6), no. 3, Dee. 
1959, pp. 526-528. 

Exact Solutions for the Propagation of 
Explosion Waves in a Gravitating Gas with 
Zero Temperature Gradient, by EK. V. 
Ryazanov, Soviet Physic s-Doklady, vol. 4 
(Russian vol. 126, nos. 1-6), no. 3, Dee. 
1959, pp. 537-540. 

On Particularities Near the Detonation 
Center and on the of Two 
Shock Waves, by N. Kochina, Soviet 
Physics-Doklady, vol. 4 ge Mean vol. 126, 
nos. 1-6), no. 3, Dee. 1959, pp. 544-548. 

Reverse In Situ Combustion, by R. P 
Gilbert, Zeitschrift fiir Angewandte Math. 
und Phys., vol. 10, no. 6, 1959, pp. 544—- 
551. 


Materials and Structures 


Nuclear Radiation Effects on Structural 
Plastics and Adhesives, PartI: Literature 
Survey, by Edwin M. Kinderman and 
Shirley B. Radding, Wright Air Dev. 
Center, Tech. Rep. 56-534, Part I, June 
1957, 71 pp. 

Nuclear Radiation Effects on Structural 
Plastics and Adhesives, Part II: Exten- 
sion of Literature Survey, by R. Y. Mixer 
and Shirley B. Radding, Wright Air Dev. 
Center, Tech. Rep. 56-534, Part II, June 
1959, 33 pp. 

Nuclear Radiation Effects on Structural 
Plastics and Adhesives, Part III: Experi- 
mental Research, by R. Y. Mixer and D. B. 
Parkinson, Wright Air Dev. Center, Tech. 
Rep. 56-534, Part III, Aug. 1957, 44 pp. 
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—— by E, J. Opik and S. F. Singer, 

Phys. Fluids, vol. 2, Nov.-Dec. 1959, pp. 
653-655. 
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Champion and R. A. Minzner, Planetary 
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264. 

Characteristics of the Radiation Emerg- 
ing from the a of a Rayleigh Atmos- 
phere, I, by K. L. Coulson, Planetary & 
Space Sci., vol. 1, Sept. 1959, pp. 265-271. 

Characteristics of the Radiation Emerg- 
ing from the Top of a Rayleigh Atmos- 
phere, II, by K. L. Coulson, Planetary & 
Space Sci. ., vol. 1, Sept. 1959, pp. 272-284. 
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Considerations, by F. F. Marmo, J. 
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291-305. 

Artificial Electron Clouds, III: Release 
of Atomic Potassium at 121 km, by F. F. 
Marmo, L. M. Aschenbrand and J. Press- 
man, Planetary & Space Sci., vol. 1, Sept. 
1959, pp. 306-318. 

The Association of Solar Radio Bursts 
with Auroral Streams, by A. Maxwell, 
A. R. Thompson and G. Garmire, Plane- 
‘ary & Space Sci., vol. 1, Sept. 1959, pp. 
325-332. 

The Photolysis of Nitrous Oxide in the 
Far Ultraviolet, by W. IX. Groth and H. 
Schierholz, Planetary & Space Sci., vol. 1, 
Sept. 1959, pp. 333-342. 

Ionizing Radiation Detected by Pioneer 
II, by A. Rosen, P. J. Coleman Jr. and 
C. P. Sonett, Planetary & Space Sci. vol. B 
Sept. 1959, pp. 343-349. 

Diffusion of Ions in a Static F, Region, 
by J. E. C. Gliddon, Quart. J., Mech. & 
Appl. Math., vol. 12, Aug. 1959, pp. 340- 
346. 

Use of Green’s Function in the Solution 
of Ionospheric Diffusion Problems, by 
J. E. C. Gliddon, Quart. J. Mech. & Appl. 
Math., vol. 12, Aug. 1959, pp. 347-353. 

Some Properties of Primary Electro- 
magnetic Disturbances Caused by Light- 
ning Discharges, by V. E. Kashprovskii, 
Telecommunications (trans. of Electrosviasz’), 
no. 7, Oct. 1959, pp. 717-731. 

Ultraviolet Reflectivities of Mars and 
Jupiter, by A. Boggess III and L. Dunkel- 
man, Astrophys. J., vol. 129, nos. 1 & 2, 
Jan. and Feb. 1959, pp. 236-241. 


Some Aspects of the Electrical Prop- 
erties of the Upper Atmosphere, by F. A. 
Goldberg, J. Aero/Space Sci., vol. 26, 
Dec. 1959, pp. 63-68. 
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A. H. Klotz, Nuovo Cimento, vol. 14, no. 1, 
Oct. 1959, pp. 135-141. 

Measurement of Phase and Amplitude 
Fluctuations of Radio Waves Which Have 
Traversed the Ionosphere, by V. V. Vitke- 
vich and Iu. L. Kokurin, Radio Engng. & 
Electronics (trans. of Radiotekhnika 7% 
Elektronika), vol. 3, no. 11, 1958, pp. 70- 
78. 


Combined Photographic and Radio 
Echo Observations of Meteors, by J. 
Davis, J. S. Greenhow and J. E. Hall, 
Proc., Roy. Soc., London, vol. A 253, Nov. 
17, 1959, pp. 121-129. 

The Effect of Attachment of Radio Echo 
Observations of Meteors, by J. Davis, 
J. S. Greenhow and J. E. Hall, Proc., 
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1959, pp. 130-139. 
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RESEARCH INGENUITY —AND NITRASOL 


That wonderful American ability to fix 
almost anything with a piece of baling wire 
often seems lost among the complexities of . 
an age in which almost any problem requires 
a nine-figure budget. 

That this “baling-wire” ingenuity is an 
unalterable American heritage was re-estab- 
lished recently by Dr. Leon Foreman, for- 
mer college chemistry professor, now in 
Grand Central Rocket Co.’s Research 
Department. 

A very interesting new solid propellant, 
Nitrasol, was being developed at GCR. It 
had great promise. But there was no eco- 
nomical way to schedule a test batch in one 
of Grand Central’s large complex mixers. 

Dr. Leon Foreman was given this prob- 
lem to study. Before anyone realized the 
significance of what he was doing, he had 
bought a polyethylene waste basket at the 
dime store, attached an air-driven shaft and 
propeller, and, with this equipment costing 
$150, mixed the first test-size batch of Nitra- 


Explorer 


sol. It was cast and fired—with complete _ 
success. al 
So efficient was Dr. Foreman’s waste- 4“ J 


basket mixing technique that GCR scaled 
up the batch size to a one-ton-a-day capacity 
in constructing the nation’s first commercial 
pilot plant for Nitrasol. With typical team 
aggressiveness, GCR personnel built the 
new plant in eight weeks. During the first 
three weeks of operation, thirty-five Nitrasol 
rocket motors were mixed, cast, and success- 
fully tested. 

Today the Nitrasol mixed in this new 
plant—revolutionary in its simplicity and 
low cost—offers a new promise to America’s 
future in propulsion. 


Positions open for chemists, engineers 
and solid rocket production specialists. 


Grand Central 
KK ockel. Co. 


P.O. Box 111 Telephone : PYramid 3-221! 
DLANOS, CALIiFORNIA 
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